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Multiple Scattering of Electrons in the Gas of a Wilson Cloud Chamber 


LyLe W. SMITH AND P. GERALD KRUGER 
University of Illinois, Urbana, Illinois 
(Received April 30, 1947) 


The deviation from linearity of straight electron tracks appearing in a cloud chamber has 
been measured and found to be in agreement with that predicted by the theory of multiple 
scattering insofar as the energy of the electrons is known. The method described could be 
easily modified to give a more valuable and exact check on the theory of multiple scattering. 





1. INTRODUCTION 


OME earlier reports’? have suggested the 
presence of a new penetrating radiation, pre- 
sumably a low mass neutral particle which decays 
at rest. There was some cloud-chamber evidence* 
that this assumption was correct. In order to 
evaluate the cloud-chamber data it was necessary 
to determine to what extent multiple scattering 
in the gas of the cloud chamber and turbulent 
motion of the gas affected the determination of 
electron energy. The effect of multiple scattering 
on curvature measurements has been discussed 
by Williams‘ and Bethe. 


2. EXPERIMENTAL ARRANGEMENT 


In order to obtain data on multiple scattering 
in a gas, the cloud chamber was operated under 
conditions which were the same as those which 
prevailed when data on the supposed decay 
electrons were taken, except that there was no 


'G. K. Groetzinger, P. G. Kruger, and Lloyd Smith, 
Phys. Rev. 67, 52 (1945). 
194s) K. Groetzinger and Lloyd Smith, Phys. Rev. 68, 55 
ws Kruger and Lyle W. Smith, Phys. Rev. 68, 104 

‘E. J. Williams, Phys. Rev. 59, 460 (1940). 

°H. A. Bethe, Phys. Rev. 70, 821 (1946). 
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magnetic field. The cloud chamber contained air 
at atmospheric pressure and three aluminum 
foils 0.80 mm thick, spaced 5 cm apart. The 
chamber was surrounded by a ring, 19 cm thick 
and 21 cm high made of lead blocks, and was 
placed about four meters from the target of the 
cyclotron and almost in a direct line with the 
deuteron beam striking the target. The cyclotron 
was completely enclosed by water tanks of 120- 
cm thickness, and the tanks on the side between 
the target and the cloud chamber contained a 
saturated solution of boric acid. 

In order to get an approximation to the energy 
of the electrons which produce the straight tracks 
which were measured, only those tracks were 
used which penetrated one or more of the 
aluminum foils. Since the minimum energy at 
which this can occur is known roughly, and the 
energy spectrum of electrons in the chamber 
under these conditions is known from data® taken 
with a magnetic field, one can assign an average 
energy to these tracks. The average kinetic 
energy will be assumed to be 1.5 Mev. 

The actual quantity which was measured on 
the straight tracks was their deviation from 


* P. Gerald Kruger and Lyle W. Smith, Phys. Rev. 72, 
360 (1947). 
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Fic. 1. Exaggerated shape of scattered track. S and ¢ 
represent the quantities measured. 


linearity. To do this the same measuring engine 
was used which was used to measure track 
curvature in the case when a magnetic field was 
present. The measured quantity was the perpen- 
dicular distance between the midpoint of the 
track and a straight line connecting the two 
ends. This is similar to the measurement of a 
sagitta of an arc in the case of curvature measure- 
ments. The following further criteria were used in 
selecting tracks for measurement. 


(1) All sections of the track must lie in the same plane, 
and this plane must lie within 10° of the median plane of 
the cloud chamber. 

(2) All sections of the track must be of the same age 
and density of ionization. 

(3) There can be no obvious single scattering over the 
region of the track which is measured. The minimum angle 
at which single scattering can be detected is estimated to 


be 0.1 radian. 
(4) Adjacent sections of the track, when extended, must 


meet in the foil. 


Pictures of the tracks were taken with a 
stereoscopic mirror system and _ reprojected 
through the same camera and mirror system onto 
a ground glass screen for measurement. Because 
of the way in which the measurements were made 
the results will give the projection of the scat- 
tering angle on a plane, not necessarily the 
median plane of the cloud chamber. 

In general, the track may have any sort of 
shape between its end points as is shown greatly 
exaggerated in Fig. 1. The quantities which are 
measured experimentally are ¢ and S. 

Rossi and Greisen’ have given a formula, 
attributable to Fermi, which expresses the distri- 
bution of S as a function of the length of the 
track and the multiple scattering angle @. By 6 we 


7B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 268 
(1941). 
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shall always mean the projection of the scatte 
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is the relation which we shall use to compare our 
results with the theory. 

The mean square average angle of scattering in 
a length of track ¢ is*® 


4nre*Z*Nt (W/mc?)? = 
= In 
(mc?) . C( W/mc?)? as 1 F Amin 


Ay 








» (3) 


where Z is the atomic number, N the number of 
atoms per cm* of the scattering substance, and W 
is the total energy of the particle. 6,,;, is de- 
termined by screening and is given by 


eZ: (W/mc?)? 
he (W/mc?)?—1 





min > (4) 
We shall take @max as the minimum angle (as- 
sumed to be 0.1 radian) which can be recognized 
as single scattering in the gas of the cloud cham- 
ber. Since from (3) we have (6?),4,~?, and we know 
from (1) that (S)*/f~ (6), we shall expect the 
relation (S)4~#! to obtain for the various lengths 
of tracks which were measured. Table I gives the 
results of the measurements where C and (S),, are 
the length of track and average deviation from 
linearity, respectively. Columns 4 and 5 give the 
experimental and theoretical values of the mean 
square scattering angle, respectively. The experi- 
mental value of (24—(.S*),)? is included in the 


TABLE I. Results of measurements on straight tracks. 
The first column gives the length of the straight track in cm 
over which the average deviation (S)s4, has been measured. 
(24-S*)y) gives the experimental value for this quantity. 
(® ye? and (6),4*%ee give the experimental and theoretical 
values of the mean square scattering angle. The last 
column gives the number of tracks on which measurements 
were made over the lengths given in column one. 











No. of 

g(cm) (S)ay(cm) (2471(S?)ay)4 (62) _yexP (62)qytheo tracks 
3 0.0155 0.0156 20.2 X10~* 19.4 x10 103 
4+ 0.0231 0.0234 25.2 X10-* 25.4 X10-* 107 
5 0.0311 0.0288 30.5 X10~* 32.4X10~* 44 
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table as an indication of the consistency of the 
measurements, since if the distribution of S is 
gaussian this quantity should be equal to (S),. 

Because of the fact that we know only an ap- 
proximate value for W, i.e., 2 Mev, and that there 
is certainly some range of energies represented by 
the measured tracks, exact agreement with the 
theory cannot be expected. For instance, if W is 
taken to be 3 Mev, the theoretical results will be 
lowered by approximately a factor of two. 

Figure 2 shows the experimental relationship 
between (5S) and ¢!. Since this plot is a straight 
line we conclude that all of the measured devia- 
tions from linearity are due to multiple scattering 
and that any turbulence of the gas in the cloud 
chamber is negligible under these conditions. 
This is in agreement with the fact that no 
turbulence is observed visually in the cloud 
chamber. 


3. CONCLUSIONS 


The distortion of straight electron tracks in a 
Wilson cloud chamber has been measured and 
found to be due to multiple scattering, within 
experimental error, the effect of turbulent motion 
of the gas not being observed. The magnitude of 
the multiple scattering is in agreement with the 
theory insofar as the energy of the electron is 
known. 
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tou? 


Fic. 2. Experimental relation between (S), and #!. 


It is obvious that the above method may be 
modified by injecting electrons of a known energy 
into a cloud chamber, e.g., from a Van de Graaff 
generator. This modification could give a more 
valuable and exact check on the theory of mul- 
tiple scattering which at present has its only 
experimental confirmation in the scattering of 
electrons by thin foils. When a foil is used as the 
scattering substances certain rather arbitrary 
criteria* must be applied to insure that no large 
angle scattering masks the multiple scattering. 
Contributions due to large angle scattering can 
be completely avoided if the path of the electron 
can be observed over the whole region in which 
the scattering is measured. 


*G. Wentzel, Ann. d. Physik 69, 333 (1922). 
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Cloud-Chamber Search for a Low Mass Neutral Particle 


P. GERALD KRUGER AND LYLE W. SMITH 
University of Illinois, Urbana, Illinois 
(Received April 30, 1947) 


The possibility that certain anomalous electron tracks in a cloud chamber may be due to 
decay of a low mass neutral particle has been examined, and it is concluded that such tracks 
are not caused by decay electrons. This conclusion is necessary because of the failure to observe 
such tracks in 5000 pictures when hydrogen gas was used in the cloud chamber instead of air. 





N order to check the results of two previous 
investigations? on the possible existence of a 
new low mass neutral particle, a Wilson cloud 
chamber containing three aluminum foils has 
been used as a detector in a search for decay 
electrons which might be emitted by such a 
particle. The arrangement of the cloud chamber 
is described elsewhere.* By observing the loss of 
an electron passing through a foil it is possible to 
determine the direction of motion of the electron. 
Thus, if an electron and a positron start from a 
common point and proceed collinearly in opposite 
directions, as would be the case if a neutral 
particle at rest decayed into a positive and 
negative electron, a decrease in energy in both 
directions should be observed provided they 
penetrate one or more foils each. In any case the 
common point of origin is not determinable ex- 
actly but should lie somewhere between or in one 
of the two foils which have been penetrated. If 
any such events occur, one must carefully evaluate 
the extent to which multiple scattering and 
chance affect their identification. Several such 
events have been observed, and it is the purpose 
of this note to describe how the validity of these 
tracks was checked. 

Tracks chosen for measurement (1) penetrated 
two or more foils and (2) satisfied other criteria.’ 
The majority of such tracks selected are electrons 
which apparently lose energy at each successive 
foil penetration. However, some anomalous tracks 
occur whose central portion exhibits a smaller 
curvature than either end portion. Also, some 
tracks were observed whose central portion was 


1G. K. Groetzinger, P. G. Kruger, and Lloyd Smith, 
Phys. Rev. 67, 52 (1945). 

$ £5) K. Groetzinger and Lloyd Smith, Phys. Rev. 68, 55 
(1945). 

3 Lyle W. Smith and P. Gerald Kruger, Phys. Rev. 72, 
357 (1947). 
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of larger curvature than the end portions, Jp 
addition to tracks which were selected according 
to the above criteria, all + electron pairs pro. 
duced by y-rays occurring in the chamber were 
measured. As a check on the chance coincidence 
of independent tracks, tracks were selected which 
penetrated one foil and ended on another, and 
which penetrated no foil and ended on one. 

Several independent determinations of the 
radius of curvature of each track were made with 
a measuring engine which determined the chord 
and sagitta of the arc of track. 

Several experimental arrangements were used. 


(1) Cloud chamber filled with one atmosphere of air and 
shielded by lead.* Tungsten target in the cyclotron (1800 
pictures). 

(2) Cloud chamber filled with one atmosphere of air and 
shielded by lead. Copper target in the cyclotron (2600 
pictures). 

(3) Cloud chamber filled with one atmosphere of air 
with reduced lead shielding. Copper target in the cyclotron 
(1940 pictures). 

(4) Cloud chamber filled with one atmosphere of hydro- 
gen and shielded by lead. Copper target in the cyclotron 
(5000 pictures). 


In all cases the cyclotron target was bom- 
barded with 10-Mev deuterons. When reduced 
lead shielding was used around the cloud cham- 
ber, the cyclotron beam was held to approxi- 
mately } of its value when maximum lead 
shielding was used. With arrangement (4) the 
cyclotron beam was approximately 20 percent 
larger than for arrangement (2). 

Table I gives the results of measurement on 
electron tracks which have penetrated two or 
more foils. Column 1 gives the track number, 
columns 2 and 4 the measured energies of the end 
portions in Mev, and column 3 the measured 
energy of the center portion in Mev. In the case 
of anomalous tracks the calculated* root mean 
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square deviations due to multiple scattering are 
given by the bracketed quantities. 

The greater number of the anomalous tracks 
can be explained by supposing them to be 
electrons, one or more portions of which suffer 
multiple scattering several times larger than the 
root mean square scattering. This in some cases is 
quite improbable but nevertheless possible. The 
number of such tracks which could be caused by 
chance can be easily estimated. The total area of 
the foils is 9800 mm?, and we assume the “‘area of 
confusion” on the foils to be 3 mm’; then, knowing 
the number of fractional parts of tracks which 
appear in the cloud chamber, one can arrive at 
the probabilities for chance coincidence for the 
arrangements (1), (2), and (3) to be 1/11,200, 
1/7050, and 1/1030 per picture, respectively. The 
calculation for arrangement (3) was checked by 
looking for ‘‘S-shaped” tracks which satisfied all 
the criteria. One such track was observed in 940 
pictures. 

The possibility that the anomalous tracks are 
due to decay electrons is indicated by the internal 
consistency of the measured energy of their 
central portions and by the fact that when lead 
shielding was removed from around the cloud 
chamber the number of electrons per picture 
increased by a factor of 3.8 while the anomalous 
tracks with larger energy center portions re- 
mained approximately constant. Since under 
these conditions the deuteron beam current was 
about } of that in arrangements (1) and (2), this 
should have resulted in a 19-fold increase in 
anomalous tracks if they had been produced by 
gamma-rays. However, when hydrogen was put 
in the cloud chamber no anomalous tracks were 
observed. The probability for this occurring is 
approximately 0.05 if the average ratio of 
anomalous tracks to electrons is the same here as 
for arrangements (2) and (3). The decrease in the 
number of electron tracks observed in this case is 
in agreement with the measured decrease in 

sensitive time of the cloud chamber when filled 
with hydrogen. 
Pair production by gamma-rays in a foil could 
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TABLE I. Results of measurements on tracks from electrons 
which have penetrated two or more foils. 








Electron energy 


Electron energy 
from measure- 


Electron energy 
from measure- 


from measure- 





Track ments on one ments on center ments on last 
num- end portion of portion of end portion of 
ber track (Mev) track (Mev) track (Mev) 
Arrangement 1 

10-1 2.3 +0.2(0.2) 3.4+0.1(0.3) 2.3 +0.1(0.2) 
13-1 1.8+0.2(0.1) 3.0+0.1(0.1) 2.4+0.2(0.2) 
17-2 3.0 +0.2(0.2) 3.4+0.1(0.1) 2.0 +0.2(0.2) 
13-2 2.740.3 2.220.3 1930.4 
15-2 2.9+0.3 2.5+0.2 2.3240.1 
170 1.6+0.1 1.1+0.1 0.7+0.1 
19-1 6.8+0.6 4.920.5 3.1+0.7 
23-4 1.9+0.1 1620.1 1020.1 
18-1 44+0.2 3.9204 3.740.2 
Arrangement 2 
34-3 1.19 +0.04(0.1) 2.19 +0.33(0.2) 0.47 pelea) 
34-4 0.08 +0.02(0.005) 2.32 +0.17(0.2) 2.10 +0.10(0.2) 
35-2 2.29 +0.09(0.1) 2.39 rey 1.91 +0.09(0.1) 
35-3 0.47 +0.07(0.1) 2.03 +0.08(0.2) 0.80 +0.07 (0.01) 
47-5 2.03 +0.12(0.2) 2.32 +0.17(0.2) 1.95 +0.16(0.1) 
50-1 2.69 +0.22(0.2) 3.13 Oni 2.52 +0.2(0.1) 
24-1 3.0 +0.2(0.1) 1.75 +0.05(0.1) 2.1+40.2(0.1) 
37-2 2.0 +0.2(0.2) 1.40 +0.05(0.1) 1.50 +0.02(0.1) 
26-3 4.3+0.5 2.9+0,2 2820.1 
27-1 3.340.3 2.340.2 2.140.2 
30-1 3.0+0.1 1.8+0.1 0.46 +0.05 
31-3 10.4+1.0 7420.5 3.240.3 
34-4 2.3+0.2 2.1240.2 0.40 +0.05 
35-1 1.7+0.2 14+0.2 1.2+0.1 
39-1 1.7+0.1 1.10+0.05 0.88 +0.04 
48-2 2.4240.2 1.8+0.2 0.50 +0.05 
Arrangement 3 
53-6 {S020 Ds0O1) 2.15 +0.13(0.2) 1.67 +0.03(0.2) 
55-15 1.50 +0.05(0.1) 1.70 +0.04(0.1) 1.34+0.05(0.1) 
55-7 3.7 ery} 2.0 +0.1(0.2) 2.2 a9 305} 
55-8 5.5 +0.5(0.6) 3.3+0.2(0.3) 3.4+0.2(0.3) 
53-3 2020.2 1.5+0.1 0.99 +0.05 
53-5 3.9+0.4 3.0+0.3 2.0+0.1 
54-3 1.7+0.1 1.5+0.1 0.72 +0.05 
55-3 2.82+0.2 1.9+0.2 1620.1 
55-9 2.1+40.2 1.5+0.1 1.1020.03 
55-10 3.5+0.2 2.6+0.2 2.1+0.1 
56-7 2.8+40.3 2.0+0.3 1.90.1 
56-10 2.9+0.1 2.0+0.1 1.30 +0.05 
56-12 2.3+0.1 2.1240.2 0.87 +0.03 
58-6 3.340.2 2.2+0.2 1.7+0.2 
58-13 4.3+0.2 3.520.4 2.920.3 
Arrangement 4 
86-3 2.0+0.1 1.55 +0.05 0.3 +0.1 
91-1 1.6+0.1 1.0+0.1 0.3+0.1 
91-2 1.7+0.1 1.2+0.1 0.20 +0.05 
102-1 §.5+1.5 2.2+0.2 1.3+0.05 
107-1 3.5+0.5 2.420.2 1820.1 








give rise to the anomalous tracks in their energy 
range only if one of the pair were scattered 
through a large angle in the foil. Also, if the 
anomalous tracks were caused by such pair 
production, their frequency of occurrence should 
have increased markedly when the lead shielding 
around the cloud chamber was decreased. 

It must be concluded that the evidence is con- 
trary to the assumption that the anomalous 
tracks are caused by the decay of a neutral 
particle. 
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On the Slow Neutron Absorption Resonances in Iodine 


Ws. B. Jones, Jr. 
Cornell University, Ithaca, New York 


(Received May 3, 1947) 


The 12-channel slow neutron velocity spectrometer at 
Cornell University is described. Reported here is a measure- 
ment with this spectrometer of the neutron cross section of 
iodine over an energy range extending from 0.0026 to 
1000 ev. Below 3 ev, the iodine cross section was found to 
obey the 1/v absorption law. The total cross section at 
thermal energy (k7=0.025 ev) was found to be 10.3 X 10-* 
cm?/atom. The scattering cross section, found by extrapo- 
lating to zero time of flight the straight line which best 
fitted the data on a cross-section vs. time-of-flight graph, 
is 3.6X10-* cm?/atom. A single resonance was observed 
at 20.3 ev. A region containing at least three resonances 
was observed between 25 and 50 ev. An indication of a 
resonance or group of resonances was found near 85 ev. 
An attempt was made to match a symmetrical type Breit- 
Wigner formula to the observed data at the 20.3 ev 
resonance. The constants chosen in this manner for the 
resonance are: E,=20.3+0.5 ev, oo=80X 10-* cm?/atom 
with o9>40X 10-** cm?/atom, and l'=0.45 ev with [<0.8 
ev. Consideration of the activities for three thicknesses 
of the absorber for this resonance showed that the absorbers 
were thick and that oof?= 15 x 10-* ev? cm*/atom, neglect- 
ing the Doppler effect. In the neutron energy region be- 


tween 25 and 50 ev, the relation between the activity and 
the absorber thicknesses used showed that the absorbers 
were thick for the resonances in this region. From the 
curve of activity as a function of thickness, it was found 
that 2- (ooI'*)t = 32.5 10-" ev-cm/atom! for these reso- 
nances, neglecting the Doppler effect. An attempt at fitting 
two resonances of about equal strength to the data proved 
to be impossible, and it was concluded that there are 
three or more resonances in this region. Consideration of 
the maximum widths of these resonances, that account for 
the absorption, indicated that the widths were less than 1 
or 2 ev. This conclusion is reached regardless of the relation 
of the level widths to the Doppler width. Therefore, it js 
believed that there are no resonances for neutron energies 
between 0 and 70 ev which are so wide as to be incom- 
patible with the present theory of nuclei. Since iodine has 
only one stable isotope, the resonances observed are for 
one type nucleus. It would be valuable to know exactly 
the number of resonances for a wide range of neutron 
energies; however, the apparatus used in the present 
experiment does not have sufficient resolution to resolve 
resonances with a spacing as small as 7 ev at 40-ev energy. 





1. INTRODUCTION 


HE slow neutron cross section of iodine has 
been investigated by a number of workers. 

The earliest measurements were made by in- 
direct methods. Because the resonance absorp- 
tion in iodine occurs for rather high neutron 
energies, it has not been possible to determine 
very accurately by indirect methods the charac- 
ter of the cross section as a function of energy. 
However, the results of experiments by indirect 
methods! indicated that the resonance-level 
widths in most nuclei were a few tenths of an 
electron volt but that the iodine resonance-ab- 
sorption width was about 15 ev. It was concluded 
from these experiments that either there is an 
unusually wide absorption resonance in iodine or 
there are perhaps several absorption resonances 
of about equal strengths and with small spacings. 
It is much more ¢ompatible with theoretical 
ideas and with much of the experimental data 
on resonances in various nuclei to attribute the 


1H. A. Bethe, Rev. Mod. Phys. 9, 134 (1937). 
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resonance absorption of neutrons by iodine to 
several closely spaced resonances rather than to 
one wide resonance. Recently, some experiments 
of a survey nature with a slow neutron velocity 
spectrometer® indicated that there is one weak 
absorption resonance in iodine of normal width 
at 20.6 ev and a region of strong absorption 
extending from 30 to 45 ev. The wide region of 
absorption was interpreted as resulting from two 
absorption resonances, because the measured 
absorption as a function of energy was not of the 
shape usually observed for a single resonance. 
However, even under the assumption of two 
resonances for this region, one cannot come to 
the conclusion that the levels are of widths ex- 
pected from the theory. 

In view of the doubt as to the nature of the 
slow neutron absorption properties of iodine, a 
more detailed investigation of iodine with a slow 
neutron velocity spectrometer was indicated. In 
addition, iodine has only one stable isotope, and 


2C. S. Wu, L. J. Rainwater, and W. W. Havens, Phys. 
Rev. 71, 175 (1947). 
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if several resonances were found, information 
would be obtained on the energy levels in a single 
nucleus. The experiment to be described here is 
a study of the neutron cross section of iodine 
with a slow neutron velocity spectrometer. 


2. APPARATUS 


The time-of-flight method and an apparatus 
used for time-of-flight measurements have been 
described earlier.* The apparatus of Baker and 
Bacher* has been improved and expanded by 
Bacher, Baker, and McDaniel‘ and by Mc- 
Daniel. The timing equipment used in the 
present experiment is, with minor improvements 
the same as that used, but not described, by 
McDaniel in the measurement of the slow 
neutron cross section of indium. 

The timing equipment of the apparatus is 
composed of electrical circuits which will turn 
on the ion source of the cyclotron for a chosen, 
definite period of time and accurately modulate 
pulse amplifiers. The cyclotron ion source is an 
arc type source and will be referred to as the arc. 
The modulated amplifiers are connected so that 
certain pulses from the linear amplifier used 
with the BF; ionization chamber are selected. 
Each modulated pulse amplifier will amplify only 
those pulses originating from neutrons detected 
by the BF; chamber within a time-of-flight 


tisai) P. Baker and R. F. Bacher, Phys. Rev. 59, 332 
1). 

*R. F. Bacher, C. P. Baker, and B. D. McDaniel, Phys. 
Rev. 69, 443 (1946). 
*B. D. McDaniel, Phys. Rev. 70, 832 (1946). 





interval determined essentially by the on-time of 
the arc, the on-time of the pulse amplifier, and 
the time separation of these on-times. There are 
12 of these modulated pulse amplifiers so that 
the intensity of the neutron beam in 12 time-of- 
flight intervals may be measured simultaneously. 
The cycle of turning on the arc and turning on 
the modulated amplifiers is repeated at a fre- 
quency which is as high as possible without 
having neutrons produced in one cycle detected 
in the succeeding cycles. 

Figure 1 is a block diagram of the timing 
equipment. The arrows on the lines connecting 
the blocks indicate the directions in which the 
signals proceed between the components of the 
equipment. In order to simplify the diagram, the 
cathode-ray oscilloscope and its associated linear 
sweep circuit, which comprise an important 
component of the equipment, are not shown. 
The oscilloscope sweep circuit is so arranged that 
the sweep may be triggered accurately in any 
phase with respect to the timing wave which 
determines the repetition frequency for the 
equipment. The sweep speed may be adjusted 
continuously from a speed of 5 us for a sweep 
across the screen of the 5-inch cathode-ray tube 
to a speed of 0.02 second for a sweep across the 
screen. The sweep speed was calibrated fre- 
quently by comparison with the time standard 
of the equipment. Tap switches are connected 
so that the vertical deflection plate of the 
cathode-ray tube may be conveniently connected 
to any one of the 35 points of the equipment 
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Fic. 2. Timing diagram for timing apparatus. 


where there are wave forms which should be 
viewed in the adjustment and checking of the 
operation of the equipment. 

The time standard is a 100 kc/sec. (+0.02 
percent) crystal oscillator. A six-stage frequency 
divider is employed to obtain the timing waves 
of the desired frequencies for determining the 
on-time of the modulated amplifiers, the time 
separation of the arc on-time and the on-time 
of the modulated amplifiers, and the repetition 
frequency. The timing wave which determines 
the on-time of the modulated amplifiers has a 
period twice as long as the on-time. This timing 
wave will be referred to as the fast timing wave. 
The detector on-times used in this experiment 
were 5, 10, 50, and 250 us corresponding to fast 
timing wave frequencies of 100, 50, 10, and 
2 kc/sec., respectively. The uncertainty in the 
detector on-times is less than 0.15 us for the 
5 us on-time and less than 1.5 percent of the on- 
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time for the longer on-times. The repetition 
periods used were 300, 320, 400, and 5000 us. 

Figure 2 shows a timing diagram for the Opera- 
tion of the timing equipment during a portion 
of a repetition period. Reference is made to 
Figs. 1 and 2 in outlining the operation of the 
equipment. 

The BF; ionization chamber gives an electrical] 
pulse upon the detection of a neutron. The 
pulses from the chamber are amplified by the 
linear amplifier, and an amplitude discriminator 
is used to pass only the neutron pulses which are 
larger than the thermal noise pulses from the 
amplifier. The rise time of the amplifier is slightly 
less than 2 ws; this is somewhat greater than 
the collection time for the electrons in the 
ionization chamber. The neutron pulses from the 
discriminator are sent through a pulse shaper 
which makes all pulses passed by the discrimi- 
nator of the same size and duration. From the 
pulse shaper the pulses are sent to a scale-of-64 
circuit driving a mechanical counter so that all 
detected neutrons are counted. From the pulse 
shaper, the pulses are also sent into the modu- 
lated amplifier J which separates the pulses into 
two sets of timed groups. One set of timed groups 
is composed of those pulses resulting from neu- 
trons detected in time intervals equal in length 
to one-half the period of the fast timing wave, 
and these intervals are separated by a time equal 
to one-half the period of the fast timing wave. 
The other set of timed groups of pulses are those 
not included in the first set of time groups. 
Figure 2a shows the fast timing wave and Fig. 2b 
the modulation envelope determining one set of 
the timed pulse groups. 

The low frequency timing wave is sent into 
phase shifters and from the phase shifters are 
obtained three pulses which may be adjusted so 
as to have any desired duration and to start at 
any desired phase with respect to the low fre- 
quency timing wave. Coincidences between these 
pulses occurring at the repetition frequency and 
pulses derived from the fast timing wave (shown 
in Fig. 2c) are used to trigger the arc modulator 
control and the two sets of 6 modulated amplifiers 
of modulated amplifier J7. Thus the accuracy of 
all the timing is determined by the fast timing 
wave. 

A selected timing pulse which triggers the arc 
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modulator control is shown in Fig. 2e. This 
pulse initiates a linear time sweep in the arc 
modulator control so the arc may be turned on 
at a phase which is continuously adjustable over 
a time interval of about twice the length of a 
fast timing wave period. This linear time sweep 
is also used to determine the on-time of the arc. 
The arc on-time is always made equal to the 
detector on-time. Figure 2f shows a typical 
phasing for the arc modulator pulse. 

The grouped neutron pulses from modulated 
amplifier J are sent into modulated amplifier JJ 
where certain groups of the pulses are selected 
to be recorded by the mechanical counters. The 
selector pulse (Fig. 2g), which selects the fast 
timing pulse (Fig. 2h) to start one of the chain 
of 6 amplifiers of modulated amplifier J7, may 
be phased to start the modulation chain with any 
one of the timing pulses. A chain of the modu- 
lated amplifier JJ selects six successive groups 
of pulses of one set of the grouped neutron pulses 
from modulated amplifier J. Thus the neutron 
pulses appearing in these six time intervals may 
be recorded on mechanical counters. The other 
chain of modulated amplifier JJ selects any six 
successive groups of pulses of the other set of 
grouped neutron pulses from modulated ampli- 
fier J. Figure 2i shows the modulation envelopes 
of the amplifiers of one chain. Thus, there are 
twelve time intervals in which detected neutrons 
are counted. These twelve time intervals are 
referred to as the detector channels. The double 
modulation scheme insures stability and uni- 
formity in the detector channel on-times. 

Since the resolving time of the mechanical 
counters and their driving circuits is relatively 
large (about 1/30 sec.), it is necessary to use 
scaling circuits so that no appreciable number of 
neutron counts will be missed by the recording 
devices. A scale-of-64 circuit is used with the 
recorder of the total neutrons detected. A scale- 
of-4 is used with each of the recorders of the 
twelve timed intervals. 

The neutron source used was the Cornell 
cyclotron. Deuterons produced by the arc are 
accelerated in the cyclotron, are brought out 
of the accelerating region of the cyclotron, and 
allowed to fall upon a Be target, thus producing 
neutrons by the (d,m) reaction. The peak deu- 
teron current falling upon the target during the 
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experimental runs was between the limits of 80 
and 150 micro amperes. The target is almost com- 
pletely surrounded with water, contained in a 
copper tank, which serves as the moderator for 
the neutrons. A diagram of the moderator tank 
is shown in Fig. 3a. 

A diagram of the neutron beam collimator and 
the neutron detector is shown in Fig. 4. The 
collimator is made up of concentric hollow 
cylinders of paraffin and BO; with an especially 
thick hollow cylinder of B:O; surrounding the 
BF; ionization chamber. This collimation is 
necessary so that the neutron source will be well 
defined; i.e., so that neutrons scattered by 
materials in the cyclotron room will not be 
scattered into the BF; chamber. The collimator 
is sufficiently effective so that a fast neutron 
background is detectable only at energies above 
100 ev. The BF; ionization chamber used as the 
neutron detector is very similar in construction 
to the chamber described by Baker and Bacher.* 
The ionization chamber has an effective length 
of 15 cm and is filled with BF; of high purity to 
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Fic. 3. (a) Moderator tank arrangement, (b) room plan. 




















































a pressure of 2 atmospheres. The chamber col- 
lecting potential used is of the order of 3 kv. 
Figure 3b is a diagram showing the relative 
positions of neutron source, collimator, and de- 
tector. 
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Fic. 4. The neutron beam collimator and neutron detector. 


A monitor for the neutron source was used. 
The monitor neutron detector is a parallel plate 
ionization chamber in which one plate is coated 
with uranium. The chamber is filled to a pressure 
of one atmosphere with argon. This chamber 
depends on the fission properties of uranium for 
the detection of slow neutrons. This type cham- 
ber is convenient and useful as a slow neutron 
monitor because of the flat plateaus obtained in 
counting rate as a function of discriminator bias 
and as a function of the chamber collecting 
voltage. During the experiment the monitor 
chamber was placed directly beneath the moder- 
ator tank. 


3. EXPERIMENTAL PROCEDURE 


The total cross section of iodine as a function 
of neutron energy was determined by measuring 
the transmission of samples of iodine of uniform 
thicknesses as a function of time of flight of 
neutrons for chosen source-detector distances. 
In most respects, the procedure for measuring 
the transmission was similar to that described in 
the references.* * 

For measurements below 12-ev neutron energy, 
the source-to-detector distance was chosen to be 
3 meters. For measurements above 12 ev, the 
source-to-detector distance was about 6 meters. 
The greatest distance used was 6.84 meters, and 
the exact value of the distance in each of the 
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high energy runs was determined by the neutron 
filter which is described below. The distance was 
measured from the effective position of the source 
to the effective position of the detector. 

The repetition period of the timing equipment 
was chosen so that neutrons produced in one 
cycle would not be detected in succeeding cycles 
For each repetition period used, it was shown 
that no significant fraction of the neutrons pro- 
duced in one cycle was detected in succeeding 
cycles. This was done by testing with a repetition 
period much longer than that to be used in the 
experimental run. With this long repetition 
period, the detector channels were set to measure 
the neutron intensities at times of flight equal 
to the times of flight for recycled neutrons for 
the repetition period and times of flight of in. 
terest for the experimental run. For measure- 
ments at neutron energies below 0.6 ev, a repeti- 
tion period of 5000 us was used with the source- 
detector distance of 3 meters. In the energy 
range 0.6 ev to 10 ev, a thick Cd (0.9 g/cm?) 
filter was used which absorbs low energy neu- 
trons, and this allowed a repetition period of 
400 us to be used at the 3-meter source-detector 
distance without neutron recycling difficulty, 
Since the average source intensity is propor- 
tional to the repetition frequency, ceteris paribus, 
the Cd filter is used to advantage in increasing 
the channel counting rates. At energies higher 
than 10 ev, a filter of very thick Cd (3 to 5 g/cm?) 
and 1.0 g/cm? of In was used. An In absorber of 
1.0 g/cm? thickness absorbs more than 99 percent 
of the neutrons with time of flight between 56 
and 66 us/m. If the source-detector distance, d, 
the repetition period, T,, and 72, a time-of-flight 
between 56 and 66 us/m, are chosen for a time- 
of-flight measurement, 7; so that r1=72—T,/d, 
the recycled neutrons of the first preceding cycle 
can be removed by a 1.0 g/cm? In filter. If 
recycled neutrons for the other preceding cycles 
can be removed with another material, such as 
Cd, an In filter can be used to advantage in 
making the repetition period short. Such a filter 
was found feasible for the high energy meas- 
urements. 

Corrections must be made for the various time 
lags in the apparatus in order to obtain the 
exact timing of the neutron time-of-flight in- 
tervals from the timing of the various modulation 
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envelopes observed on the oscilloscope. One cor- 
rection is the time lag obtained from the fast 
neutron time distribution function measured 
with the apparatus. This function will be referred 
to as the FN-function, and the method of ob- 
taining it is explained below. The FN-function 
gives the sum of the time lags of the apparatus 
involved in neutron production and detection. 
Another correction is the mean time of the 
neutrons in the moderator. A third type of 
timing correction is the effect of the variation 
of the neutron source intensity with energy. 
This correction is small if the intensity variation 
is small within the resolution width of the 
apparatus. The resolution was always made 
sufficiently good so that this correction was 
negligible. 

To obtain the FN-function, the BF; chamber is 
replaced by a similarly constructed butane filled 
chamber. Such a chamber detects fast neutrons 
by detecting the recoil protons. The FN-function 
is obtained by measuring the fraction of the total 
neutrons detected which are counted in a de- 
tector channel as a function of the time delay 
of the channel, with respect to the arc modulator 
pulse. For this measurement, the water is re- 
moved from the moderator tank. Not only does 
the FN-function indicate the time lag, but also 
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the sharpness of the arc modulation and detection 
of the neutrons by the equipment. The ideal 
shape of this function is that of an isosceles 
triangle with a base twice as long as the on-time. 
Figure 5 shows a typical FN-function obtained 
with an on-time of 5 us. The dotted curve shows 
the ideal shape of the function. The “‘tailing”’ of 
this function on either side is caused, for the most 
part, by the finite rise time of the linear amplifier 
and the straggling in the acceleration of the 
deuterons in the cyclotron. Thus, for all on-times, 
the “tailing” is about of the same duration and 
the FN-function is, therefore, more nearly tri- 
angularly shaped for longer on-times. The time 
lag is taken as the delay of the centroid of this 
function with respect to the center of the arc 
modulator pulse. The FN-function was taken at 
some time during each set of observations, and 
the delay was always found to be within 0.15 us 
of 7.4 us for the 5 us on-time and very nearly 
7.5 us for the other on-times used. 

For energies above 0.7 ev, the mean time of the 
neutrons in the moderator is taken as the mean 
deceleration time, assuming the protons in the 
water are not bound. For example, this time is 
1.5 ys for 1-ev neutrons. At energies below 
0.04 ev, the mean time of the neutrons in the 
moderator is taken as the decay time constant 
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Fic. 5. FN-function or observed time distribution of fast neutrons for 5 ws on-time. 
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TABLE I. Important features of the experimental conditions. 














Repetition Absorber 
Run On-time Distance period Source modifier thickness Time of flight Ener 
No = 2 m | hae IbiL cf 20 a ter g/cm? us/m ad 
250 3 5000 Tray 16.7 490-1420 0.0026-0.001¢ 
2 50 2.97 5000 Tray 31.0 290-477 0.023-0.063 
3 50 3 5000 Tray 31.0 98-274 0.0695~.549 
4 10 3 400 \ 0.9 g/cm? Cd 31.0 64-101 0.51-1.28 
5 10 3 400 0.9 g/cm? Cd 31.0 23.8-60.6 1.42-9,2 
6 5 6.84 300 5 g/cm? Cd 16.2 12.8-20.85 12-31.9 
1 g/cm? In 
7 5 6.32 300 5 g/cm? Cd 5.96 9.09-17.8 16.5-63.3 
1 g/cm? In 
8 5 6.32 300 5 g/cm? Cd 2.5 9.09-17.8 16.5-63.3 
1 g/cm? In 
9 5 6.6 320 5 g/cm? Cd 9.78 5.66-14.0 26.6-164 
1 g/cm? In 
10 5 6.04 320 3.4 g/cm? Cd 16.2 2.07-11.2 41.6-1200 
1 g/cm? In 








of the time distribution of neutrons emerging 
from the moderator. Using water alone as the 
moderator gave a value for the mean time of 
85 us. Therefore, for measurements in the low 
energy region, the moderator was modified by 
placing a sheet of thick Cd directly on the face 
of the moderator tank to remove the low energy 
neutrons emerging from the tank and placing 
adjacent to this a sheet of paraffin 2 cmi thick 
to regenerate low energy neutrons. This paraffin 
and Cd sheet combination will be referred to as 
the “tray.’’ The mean time of neutrons emerging 
from the tray is 30 us. In the energy region be- 
tween 0.04 ev and 0.7 ev, it was assumed that 
the mean time of the neutrons in the tray varied 
linearly with time of flight from 30 us at 0.04 ev 
to 1.8 us at 0.7 ev. 


4. RESULTS 


The transmission by several thicknesses of 
iodine absorbers was measured over the energy 
interval from 0.0026 ev to 1000 ev, corresponding 
to a time-of-flight interval from 2.2 us/m to 
1400 us/m. A resonance was found at 20.3 ev. 
A region of high cross section was found between 
29 and 46 ev, and an analysis shows that this 
region must contain at least 3 resonances of 


about equal strength. Indications of resonances 
at higher energies were found. 

The iodine used for the absorbers was re- 
sublimed iodine crystals with a minimum purity 
of 99.985 percent. Brass containers were used 
for the iodine and a “blank” of brass of the 
same thickness as the sum of the thicknesses of 
the two faces of the containers was interposed in 
the neutron beam when counts were being taken 
with no absorber. The containers were packed 
tightly with the iodine in such a way that each 
absorber was believed to be of uniform thickness. 

Table I gives the important features of the 
experimental conditions under which the data 
were taken. All the data were corrected for back- 
ground and time lags. Each separate run is 
assigned a run number for the purpose of refer- 
ence. On all graphs showing experimental data, 
a legend is given which indicates the run number, 
on-time, source-detector distance, and absorber 
thickness. Thus the legend ‘1-250-3m 31 g/cm*” 
refers to run number one in which the on-time 
was 250 us, the source-detector distance was 3 
meters, and the absorber thickness used was 
31 g/cm’. 

In the energy range below 10 ev, measure- 
ments were made with a resolution which was 
sufficiently high to detect possible resonances. 
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Fic. 6. Low energy neutron cross section of iodine. 


Sufficient accuracy was obtained so that the 
l/v law would be tested. Figure 6 shows a cross- 
section plot for times of flight greater than 
90 us/m. The statistical probable errors in the 
points are indicated by the lengths of the vertical 


lines through the points, or for those points with 
no vertical line the probable error is no greater 
than the diameter of the point. No correction 
for the effect of the apparatus resolution function 
on the transmission was made. No such correc- 
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Fic. 7. Neutron cross section of I in the 0.5 to 15-ev region. The dotted line is the best 


straight line fit for the low energy data. 
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Fic. 9. Observed transmission of I absorbers at high energy. 


tion is necessary, at least for the most interesting Breaks in the curve are noted at times of flight 
part of the curve, because the transmission varies of 970 us/m and 1200 us/m, corresponding to 
nearly linearly with time of flight within a neutron wave-lengths of 3.9 and 4.8A, respec- 
resolution width, and if the transmission is_ tively, which were interpreted as the effect of 
linear with time of flight, the measured trans-__ the crystalline structure of iodine on the scatter- 
mission is equal to the actual transmission. ing cross section. Except for the crystalline 
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effects, the points lie very nearly on a straight 
line, thus verifying the 1/v absorption law. The 
equation for the best straight line through the 
ints, excluding the crystalline effects, is given 
by «= (3.6+0.0146t) X 10-*4 cm?/atom, where ¢ is 
the time of flight in us/m and o is the cross 
section. This is the same as o=(3.6+1.06E-) 
10-4 cm*/atom, where E is the energy in ev. 
The error in the slope of this curve is believed 
to be no more than 4 percent. The intercept of 
this curve is taken as the scattering cross section 
at the energies above 0.5 ev. At energy kT 
(0.025 ev), c= 10.3 X 10-*4 cm?/atom. 

Figure 7 shows the time-of-flight region be- 
tween the low energy resonance and 100 us/m. 
Again, no correction has been made for the 
effect of the resolution function on the cross 
section. The statistical errors are shown in the 
usual manner. The dotted straight line has the 
same slope and intercept as the straight line 
which best fits the data shown in Fig. 6. 

Figures 8 and 9 show the measured trans- 
mission as a function of time of flight for the 
data taken above 14 ev. Data were taken with 
several thicknesses of absorber in order to deter- 
mine the character of the resonances indicated 
in these figures. The curves through the points 
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serve to indicate the probable shape of the 
transmission curve as seen by the spectrometer. 
The curves indicate that there is a single absorp- 
tion resonance near 16 us/m. The flat region of 
low transmission of the curves between 11 and 
12 us/m contains a number of unresolved reso- 
nances. At 8yus/m there is an indication of 
another resonance or group of resonances. At 
higher energies, there are indications of other 
resonances. 

In a region where the transmission varies non- 
linearly by a large amount within a resolution 
width of the apparatus, it is expected that large 
corrections have to be made in the data to obtain 
the true transmission of the absorber. With the 
best resolution used with the present apparatus 
there are large variations in the cross section of 
iodine within a resolution width, in the region of 
the resonances. The variation in the transmission 
can be reduced by making the absorber thin, but 
it is very difficult to obtain a good indication of 
the cross section from the measured transmission 
of a very thin absorber. One can obtain some 
information on the character of the cross section 
by measuring the transmission for several thick- 
nesses of the absorber as was done in the region of 
the resonances. 
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Fic. 10. Cross section of I in the region of the main resonances calculated from the formula 


T=exp(—on), with no correction for apparatus resolution. 
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FiG. 11. Cross section of I as a function of neutron energy for the energy region studied. No 
correction has been made for apparatus resolution. 


Figure 10 demonstrates the effect of changing 
the absorber thickness in transmission measure- 
ments. The cross section was obtained from 
the observed transmission in the formula 
T =exp(—on), where T is the transmission and n 
is the absorber thickness in atoms/cm?. The cross 
section thus obtained is plotted as a function of 
the logarithm of the energy. It is noted that the 
apparent cross section increases with decreasing 
absorber thickness in the region of the resonances, 
as is expected. Thus the best indication of the 
cross section in the region of the resonances is 
given by the thinnest absorber. One fact demon- 
strated by this graph is that there are large 
variations in the cross section within a resolution 
width in the energy region between 30 and 45 ev, 
where the spectrometer indicates little variation. 
Of course, in regions where there are no large 
variations of cross section within a resolution 
width, the cross section is given best by the thick 
absorber data. 

Figure 11 shows a plot of the cross section as a 
function of energy for the whole energy range 
investigated, the abscissa having a logarithmic 
scale. Only those data are plotted which give the 
best representation of the cross section in the 
regions where data were taken for several thick- 





nesses of absorber. However, it must be re- 
membered that in the region of the resonances the 
best measured values of the cross section are sub- 
ject to large corrections because of the finite 
resolution width. Using a logarithmic energy 
scale shows the data in such a way that the reso- 
lution width is the same order of magnitude over 
the whole energy range. 


5. DISCUSSION 


For the most part, this discussion is concerned 
with the treatment of the data in the region of the 
resonances. In the region below 15-ev neutron 
energy the variations in the cross section are 
sufficiently slow (except for the crystalline effects) 
so there is no significant difference between the 
actual transmission and the observed transmis- 
sion of an absorber. In the region of the resonances 
the corrections required in the observed trans- 
mission to obtain the actual transmission are 
very large. 

Baker and Bacher® have given a treatment of 
the problem of obtaining the transmission curve 
for an absorber from the transmission observed 
with a slow neutron velocity spectrometer. They 
demonstrate that it is rather difficult to proceed 
directly from the observed transmission to the 
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actual transmission of an absorber. A much 
simpler procedure is to assume a certain trans- 
mission, find the effect of the resolution of the 
apparatus on the assumed transmission, and com- 

this result with the observed transmission. 
This process is repeated until a good comparison 
with the experimentally determined transmission 
is found. In the region of a resonance, the trans- 
mission assumed is that obtained from the Breit- 
Wigner single-level formula® with a correction for 
the Doppler effect, and the parameters in the 
formula are adjusted for the best correspondence 
to the observed data. 

The resolution function is the sensitivity of the 
apparatus, for any one channel, as a function of 
the neutron time of flight. If the resolution func- 
tion were determined solely by the on-times, its 
shape would be that of an isosceles triangle with 
the base of the triangle twice the on-time. How- 
ever, several other factors influence the shape of 
the resolution function. The combined effect of 
the on-time of the arc and detector channels, the 
straggling in the arrival of deuterons at the Be 
target, and the uncertainty in the time of detec- 
tion caused by the finite rise time of the linear 
amplifier is shown by the FN-function. A typical 
FN-function is shown in Fig. 5 for an on-time of 
5 us. It is seen that by far the largest contribution 
to the spread of the FN-function is the on-time. 

Contributions to the spreading of the resolution 
not measured by the FN-function are the uncer- 
tainty of the time spent by the neutrons in the 
moderator, the uncertainty in the source position, 
and the uncertainty in the position of the detec- 
tion of the neutrons in the BF; chamber. Under 
the conditions of interest here, i.e., for neutron 
energies greater than 15 ev and distances of flight 
greater than six meters, only the last of these 
three contributions is of significance. The neutrons 
are almost equally likely to be captured at any 
point along the 15-cm length of the BF; chamber. 

The time required for a 20-ev neutron to traverse 
15 cm is 2.4 ws. The uncertainty in the source 
position is less than one centimeter, and this is 
small compared with 15 cm. The uncertainty in 
the time of the neutrons in the moderator is a few 
tenths of a microsecond, and this is small com- 


*H. Feshbach, D. C. Peaslee, and V. F. Weisskopf, Phys. 
Rev. 71, 145 (1947). 
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pared with the time for a neutron to traverse the 
chamber length. Therefore, the resolution func- 
tion is determined essentially by the FN-function 
and the simple step-function representing the 
sensitivity of the apparatus as a function of 
distance. 

Figure 12a shows the apparatus resolution 
function for an on-time of 5 us, a source-detector 
distance of 6.5 meters, and for 20-ev neutrons. In 
making comparisons between the Breit-Wigner 
formula and the data, it is convenient to use a 
step-function approximation for the resolution 
function. Figure 12b shows the step-function ap- 
proximation used for the resolution function 
shown in Fig. 12a. The step-function is made 
symmetrical with respect to the line corre- 
sponding to zero time. Also, the second moment 
of the step-function taken at zero time is made 
equal to the second moment of the resolution 
function taken with respect to the vertical line 
through its centroid. 

A symmetrical form of the Breit-Wigner single- 
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Fic. 12(a). Apparatus resolution function for 54s on- 
time, 20-ev neutrons, and 6.5-m distance of flight. (b). Step- 


function approximation for resolution function above. 
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Fic. 13. Comparison of the Breit-Wigner formula with experimental points for three 
thicknesses of absorber used for the 20.3-ev resonance. 


level formula was used of the following form: 
ool? 


oa(E) = ’ 
lr?+4(E—E,)? 





where o,=neutron cross section, E=neutron 
energy, E,=resonance energy, oo=cross section 
at resonance, and [' = level width. Thecross section 
given by this formula is modified by the Doppler 
effect resulting from the thermal motion of the 
absorber atoms. Thus the effective cross section 
is given by 


o(E) =(1/x!A) f oa(E’) 


Xexp[ —(E—-E’)?/A* dE’, 


where A = 2(E,Eym/M)'! = Doppler width, 
m = neutron mass, Ey, = thermal energy, M = mass 
of target nucleus. Since the Debye temperature 
for iodine is 106°K, the thermal energy is taken 
as kT, in accordance with the calculation by 
Lamb,’ to take into account the lattice binding. 
The value of A is then 0.13 ev for E,=20 ev. 

The low energy resonance is now considered. In 


7W. E. Lamb, Phys. Rev. 55, 190 (1939), 





order to make a comparison between a theoretical 
curve and an observed transmission, the ob- 
served transmission curve must be corrected for 
the scattering cross section. The value of the 
cross section obtained by extrapolating the curve 
which gives the best fit in the 1/v region, to zero 
time of flight, is used as the scattering cross 
section. Thus the observed transmission in the 
region of the resonance was divided by the trans- 
mission corresponding to a scattering cross sec- 
tion of 3.6X10-* cm?/atom. Values are chosen 
for the parameters of the Breit-Wigner formula 
to give o,. From og, o is then determined. For 
each absorber thickness, « was used to determine 
the transmission as a function of time of flight for 
the source-detector distance used in taking the 
data under consideration. The step-function ap- 
proximation to the resolution function is then 
used to obtain the transmission curve due to 
absorption that would be observed for a cross 
section given by the Breit-Wigner formula with 
the chosen parameters. Such transmission curves 
were compared with the observed curves and 
parameters adjusted until the best correspond- 
ence was obtained. 

Figure 13a, b, and c shows the transmission 
curves obtained from the theory with a particular 
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set of parameters which gives the best corre- 
spondence to the observed data for the three 
absorber thicknesses used on the low energy 
resonance. The solid curves are the transmissions 
given by the theory, the dotted curves are these 
transmissions after being put through the resolu- 
tion function, and the crosses show the experi- 
mental data corrected for the scattering cross 
section. The values of the parameters giving this 
correspondence are E,=20.3 ev, oo=80X10-* 
cm?/atom, '=0.45 ev. The error in £, is about 
0.5 ev. Because of the magnitude of the experi- 
mental errors, the values of oo and I may be 
varied over quite a wide range, with correspond- 
ence between the theory and experiment within 
the experimental error. Also, because of the 
Doppler effect, no minimum value of I or maxi- 
mum value of oo can be given. The maximum 
value of I and minimum value of 9, which give a 


’ theoretical curve which could be said to corre- 


spond to the experimental points, are 0.8 ev and 
40X10-** cm*/atom, respectively. 

Another method of analysis which is useful 
yields information on the values of the parameters 
in the Breit-Wigner formula, and is largely inde- 
pendent of the resolution of the apparatus. The 
method has to do with the behavior of the 
activity integral for a resonance as a function of 
absorber thickness. The activity, I, for a reso- 
nance and a particular absorber thickness is 


given by 
I= fa —T)dy, 


where T is the transmission of the absorber as a 
function of time of flight ¢, and dy=dt/t,, where 
t, is the time of flight corresponding to the £, for 
the resonance. The limits of the integral are taken 
as values of ¢ on either side of the resonance where 
T ~1. For a given absorber thickness the value of 
| is the same for any resolution that is used to 
measure T if the resolution is sufficiently good to 
separate the resonance in question from other 
resonances. 

If the symmetrical form of the Breit-Wigner 
formula is used, then 


nool? 
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If noo is sufficiently small so that the exponential 
term in the integrand may be approximated by 
the first two terms of a power series, then 


] = rnT'00/4E,. 


No error is made in this equation by neglecting 
the Doppler effect. Thus the activity is pro- 
portional to the absorber thickness for thin 
absorbers: 

If moo>1, little error is made by assuming 
I’«(E—E,)* in computing the activity using the 
Breit-Wigner formula. Making this approxima- 
tion, one obtains 


[= (rnI a9) 4/2E,. 


Thus for thick absorbers, the activity is pro- 
portional to the square root of the absorber 
thickness. 

If the apparatus resolution is so wide that 
levels cannot be resolved, if noo>1 for each level, 
and if the separation of the levels is large com- 
pared with the level widths (as must be the case 
for the Breit-Wigner formula to hold). 


I=(n)! >, (oo1',?/4E,*)}, 


where r denotes parameters characterizing the 
rth level, and the sum is taken over the unresolved 
group. Again, for thick absorbers, the activity is 
proportional to the square root of the absorber 
thickness. The Doppler effect introduces an error 
in all these formulas for I for thick absorbers. 
However, if AI, for the resonances or if the 
absorber is sufficiently thick, the error is neg- 
ligible. 

In Fig. 14a is shown a plot of activity calcu- 
lated from the experimental curves (corrected for 
scattering cross section) as a function of the 
square root of the absorber thickness for the low 
energy resonance. With the possible exception of 
the point for the thinnest absorber, these points 
fall very nearly on a straight line. From the slope 
of this curve, one finds that oI?=15x<10-* 
ev? cm?/atom, assuming no error to be introduced 
by the Doppler effect. This is in agreement with 
the value of 17X10~-** ev? cm?/atom, by use of 
the values of o») and I obtained from the curve- 
fitting method described above. The Doppler 
effect is neglected because an analysis taking into 
account this effect is difficult, and the data are 
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Fic. 14. Activity vs. square root of absorber thickness. 


not sufficiently accurate to warrant this refine- 
ment of the analysis. 

The energy region 25 to 50 ev is now con- 
sidered. If the cross section were constant be- 
tween 30 and 40 ev, the apparatus resolution is 
sufficiently good so that the measured trans- 
mission would be an exponential function of the 
thickness. That this is not the case is demon- 
strated in Fig. 10. Therefore, it may be concluded 
that the 25- to 50-ev region must contain at least 
two resonances. By making use of the curve- 
matching method it was not found possible to fit 
two resonance curves of about equal aol to the 
data for all thicknesses of absorber used. One sees 
that this is evidently the case from two con- 
siderations. The width of the dips in the trans- 
mission curves are about the same for all thick- 
nesses of absorber. This indicates that there must 
be narrow resonances near the edges of the dip in 
transmission. However, two narrow resonances 


JONES, 











JR. 


near the edges of the dip would not give a flat 
bottom to the dip for the thick as well as the thin 
absorbers. Therefore it is concluded that in this 
energy region there are at least 3 resonances, No 
attempt was made at curve matching to de. 
termine the characteristics of 3 resonances that 
might fit the data, because there are insufficient 
data to justify such an attempt. There js no 
method of deducing from the data the maximym 
number of resonance levels in this energy region, 

Shown in Fig. 14b is a plot of activity for the 
25- to 50-ev region as a function of the square root 
of the absorber thickness. The points fall very 
nearly on a straight line, indicating that the 
absorbers are thick for all or most of the reso. 
nances in this region. One finds from this curye 
that >. (col)! for the resonances is 32.5 x 10-1 
ev-cm/atom}!. This value is given under the as. 
sumption that 38 ev is near the resonance energy 
for each resonance and that no error is introduced 
by the Doppler effect. 

The data on the 25- to 50-ev region are now 
reviewed to determine what evidence exists for 
believing that the resonances in this region have 
widths of the order expected from the theory. It 
is assumed that the Breit-Wigner formula holds 
for these resonances. This means that the widths 
of the resonances must be small compared with 
the spacing. If, for example, it is supposed that 
there are 3 levels with about 10-ev spacing 
of about equal strengths, ooI?=120X10™ 
ev? cm?/atom for each level. The data indicate 
that oo for such resonances is large compared with 
40X10-* cm?/atom. Thus, for this possible 
distribution of levels, the widths are no greater 
than 1 ev. If it is supposed that there are 3 levels 
unequally spaced, then the levels nearest each 
other must be weak and narrow because of their 
small spacing, and the third level could have 
gol? as large as 200 10-* ev? cm?/atom corre- 
sponding to a width perhaps as great as 2 ev. 
Thus, to correspond to the data, no one of the 
three supposed levels could have a width of more 
than 1 or 2 ev. If it is supposed that there are four 
or more resonances within this region, it is evi- 
dent that the level widths can be no more than | 
or 2 ev. The conclusions drawn above are based 
on the assumption that >A. On the other hand, 
the level widths may be about the same as, or 
smaller than, the Doppler width, A=0.19 ev, at 
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this energy. However, if this is so, then one is led 
to the same conclusion, namely, that abnormally 
wide levels are not indicated by the data. 





With regard to the energy region above 60 ev, 
little can be said about the constants of the 
resonances. It seems evident that there is at least 
one resonance in the neighborhood of 85 ev. 
However, the resolution of the apparatus is not 
sufficiently good at this and higher energies to 
make it possible to obtain more information than 
this from the data. 

Because of the fact that the iodine resonances 
occur at rather high energies, it is not expected 
that the data obtained by indirect methods 
should compare well with the present data. The 
recently reported results of Wu, Rainwater, and 
Havens,” using the Columbia neutron spectrome- 
ter on iodine, agree with the results obtained here 
with regard to the energies for the resonance ab- 
sorption and with regard to the slope and _inter- 
cept for the best straight-line fit for the data in 
the low energy region. However, their result for 
the value for ool of about 4 10- ev? cm*/atom 
for the 20.3-ev resonance is only } the value ob- 
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tained here. They assume that the 38-ev group of 
resonances is composed of two levels of equal 
strengths with ool = 135 X 10- ev? cm?/atom. In 
the present experiment it is found that there 
must be at least three resonances in the 38-ev 
region. If it is assumed that there are only three 
resonances of equal strength, the data presented 
here indicate that ool’ = 115 x 10-* ev? cm?/atom 
for each resonance level. Using their data and 
postulating three levels of equal strength, one 
obtains ool’ =60X10-* ev? cm?/atom. There is 
no apparent reason for these discrepancies unless 
an error has been made or unless the resolution of 
their apparatus is not sufficiently good to sepa- 
rate the resonance groups in the particular way 
they are separated in the present experiment. 
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Professor P. Morrison, and many other members 
of the Laboratory of Nuclear Studies at Cornell 
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carrying out this investigation. Support under 
contract N6-ORI-91 of the Office of Navy Re- 
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Energy Levels in the Nucleus Mn‘*** 


A. B. Martint 
Yale University,** New Haven, Connecticut 


(Received May 15, 1947) 


The energy distribution of the protons emitted in the reaction Mn**(d,p)Mn* has been 
studied by an absorption method. For excitations up to 4.38 Mev above the ground state of 
the Mn* nucleus, the protons are found to fall into six distinct energy groups, corresponding 
to excited states at 1.07, 1.77, 2.48, 3.61, and 4.38 Mev above ground. Results are also presented 
which indicate that the yields of the individual proton groups vary slowly with deuteron 
energy, showing no sharp resonance effects. From the Q-value for the most energetic proton 
group the mass of the nucleus Mn* has been computed to be 54.9634+.0022 mass units. 





INTRODUCTION 


NUMBER of investigators! have observed 
the energy levels in nuclei of low mass 
number, at excitations up to about 5 Mev, which 
are exhibited in the process of (d,p), (d,m), 
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Fic. 1. Sketch of the bombardment chamber, mechanical 
foil changers and proton counters used in measuring the 
range of protons emitted in a direction perpendicular to 
the direction of the deuteron beam. The target is manga- 
nese metal (Mn**) evaporated on gold foil. 
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1 For example, J. D. Cockcroft and W. B. Lewis, Proc. 
Roy. Soc. 154, 246 (1936); W. L. Davidson, Jr., Phys. 
Rev. 56, 1062 (1939); E. B. M. Murrell and L. C. Smith, 
Proc. Roy. Soc. 173, 410 (1939); Ernest Pollard and 
W. W. Watson, Phys. Rev. 58, 12 (1940); T. W. Bonner, 
Proc. Roy. Soc. 174, 339 (1940); M. G. Holloway and 
B. L. Moore, Phys. Rev. 58, 847 (1940). 
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(a,p) and other particle emission reactions, The 
results do not appear to be in agreement with 
generally accepted theories,? which predict a 
greater density of levels than is observed and, 
also, that the level density should increase 
rapidly with excitation and with atomic weight. 
This experimental method has been applied 
mainly to the study of elements of low atomic 
weight, and the energy level spacing observed in 
these cases is on the average about 1 Mev or 
somewhat greater; moreover, no uniform trend 
of level spacing, either with excitation in a 
particular nucleus or, at a given excitation, with 
atomic weight, is apparent. It is of interest to 
extend studies of this type to elements of higher 
atomic weights to determine if better agreement 
with theory can be obtained for particle emission 
reactions of this type, particularly with regard 
to the theoretical prediction of smaller level 
spacings in heavier nuclei. 

In the present work, the energy levels in the 
aucleus Mn** have been determined by measur- 
ing the energies of the protons that are emitted 
when the element Mn* is bombarded with 
deuterons. These nuclei contain a sufficient num- 
ber of particles to make the statistical approach 
of theory at least reasonably valid. The assign- 
ment of levels to Mn* is unique since manganese 
contains only a single isotope. 

In addition to the main problem of measuring 
proton energies, a study has been made of the 
yields of the individual proton groups as 4 
function of deuteron energy. 


2H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 173 
(1936); H. A. Bethe, ibid. 9, 86 (1937); John Bardeen, Phys. 
Rev. 51, (1937). 
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EXPERIMENTAL ARRANGEMENT 


The experimental arrangement used to meas- 
ure proton range is shown in Fig. 1. The target 
is mounted in the outer bombardment chamber 
of the cyclotron and is inclined at an angle of 45° 
to both the incident beam of deuterons and to 
the direction of the protons that are to be 
counted. The geometry of the system is such 
that only protons emitted at angles of 90+4° to 
the direction of the deuteron beam will be 
detected. 

The target material was 99.986 percent pure 
manganese obtained from Electro Manganese 
Corporation, Knoxville, Tennessee. Principal 
impurities were 0.012 percent sulfur as sulfide, 
0.001 percent sulfur as sulfate, and less than 
0.001 percent iron. Thin targets of uniform 
thickness were prepared by evaporating the 
manganese metal on gold foils. The target used 
to obtain most of the data to be presented here 
had a thickness of 1.0 cm air equivalent. A 
second target, having a thickness of 0.2 cm air 
equivalent also was used to obtain data on the 
high intensity groups and to study the effect of 
target thickness on the resolution of the method. 

Various thicknesses of aluminum foils are held 
in position between the target and the proton 
counters by the mechanical foil changers, shown 
in Fig. 1, which are operated remotely from the 
cyclotron control room. Also, a considerable 
portion of the data was obtained using a gas 
absorption cell in place of the foil changers. 
This cell was designed so that one of its end 
plates could be bolted directly to the exit port of 
the bombardment chamber and the other end 
plate could be fastened directly to the face of 
the first counter, thus eliminating all air spaces 
in the path of the protons and providing a 
continuous control over the total absorption. 
A pressure line made up of one-half inch copper 
tubing was run from the cyclotron to the control 
room so that the pressure of air in the cell could 
be adjusted from the operator’s position. 

Two proton counters, operating into a coinci- 
dence circuit, are used in order to reduce the 
background count due to the high flux of neu- 
trons in the vicinity of the cyclotron. The 
counters are operated in the proportional region, 
and the counting level of the second one is 
normally set to a comparatively high level in 
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Fic. 2. Detailed drawing of the proton counter. The 
preamplifier unit is attached directly to the base plate of 
the brass cylinder enclosing the counter. Aluminum foils, 
held in place by rubber gasket seals, normally cover the 
counter ports. 


order to peak up the proton groups. The gas 
used in the counters was tank argon; counter 
pressures of 10, 20, 40, and 65 cm Hg were 
employed in the course of this work, the counting 
level being changed in each case to achieve the 
highest possible resolution. 

A detailed drawing of the proton counter is 
shown in Fig. 2. The counter cylinder is of brass 
and is 3 cm long and 1.9 cm in diameter. The 
central wire is 10-mil tungsten; a small glass 
bead encloses the free end. High voltage is 
applied to the counter cylinder, and a negative 
pulse is taken from the central wire. A small 
preamplifier unit is attached to the base plate of 
each counter case so that the lead from the 
central wire to the grid of the first amplifier tube 
is made as short as possible, thus keeping the 
capacity of the wire to a minimum. 

A circuit diagram of the preamplifier unit is 
shown in Fig. 3. This unit is essentially an 
impedance transformer, which takes pulses from 
the counter wire at a high impedance level and 
has a cathode follower output stage matched to 
a 75-ohm line. The pulses are then piped over 
approximately 60 feet of shielded cable to the 
main amplifier and the remainder of the counting 
system, which are located in the cyclotron control 
room. The feedback from the plate of the second 
tube, through the variable condenser, to the grid 
of the first, is a capacity neutralization scheme 
which effectively reduces the capacity of the grid 
input circuit, enabling a sharp pulse (about 1- 
microsecond half-width), with a fast rate of rise 
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Fic. 3. Circuit diagram of the preamplifier unit. The feedback from tube 2 to tube 1 


pooonsiny | reduces the capacity of wed og input circuit and sharpens the pulses from the 


counter. Negative feedback is employ 
fluctuations and tube characteristics. 


(0.1 microsecond or less), to be obtained from 
the counter. 

The portion of the counting system at the 
operator’s position consists of two main video 
amplifiers, a coincidence circuit, a slightly modi- 
fied Higinbotham scale-of-64 circuit, and the 
usual mechanical recorders. The main amplifiers 
have a band pass of about 4 megacycles and a 
gain of 80 db. The response characteristics favor 
high frequencies; the low frequency response is 
made poor intentionally to reduce the effect of 
microphonics, hum, and similar sources of 
interference. 

The coincidence circuit employs a pentode 
tube. The output of one of the counters is 
applied to the control grid of the pentode; the 
output of the other goes to the screen grid. Both 
grids are biased beyond cut-off so that the tube 
will conduct only when positive pulses are applied 
to the two grids simultaneously. The resolving 
time of the cgjncidence circuit is of the order of 
one to five microseconds. This figure was deter- 
mined by measuring the number of random 
coincidences, Nr, occurring in the interval during 
which JN, counts were recorded in counter 1 and 
Nz counts registered in counter 2, the counters 
being activated by independent particles. The 
resolving time 7 is then given by Ne=2rN1N2. 

As mentioned above, a number of different gas 
pressures were used in the counters in the course 
of these observations, in an effort to attain higher 
resolution, and at each pressure it was necessary 
to readjust the counting level in order to obtain 


in the cathode circuits to stabilize against voltage 


the best results; that is, the maximum peaking 
of the proton groups consistent with a reasonable 
counting rate. The change in counter “depth” 
resulting from the use of different counting levels 
was taken into account by a method similar to 
that described by Holloway and Moore? A 
number of calibrations using ThC’ alpha-particles 
and an analytical study based on proton counting 
observations were both used to determine this 
correction term. 

The range-energy relation for protons was 
taken from the curves of Livingston and Bethe.‘ 


ENERGY AND HOMOGENEITY OF THE 
DEUTERON BEAM 


It was necessary to determine as accurately as 
possible the energy of the deuterons from the 
cyclotron, since the energy of the bombarding 
particle enters into the calculation of the disinte- 
gration energy (the Q-value) of a nuclear reac- 
tion. The homogeneity of the deuteron beam is of 
interest in connection with the resolution of this 
method of inquiry. Measurements were made on 
the range of the deuterons by causing the beam 
to pass through an absorption cell, placed at the 
target position, and then into a detection cham- 
ber. Two types of detectors were used in this 
study. One was a simple Faraday chamber, in 
which the amount of beam collected on an in- 
sulated probe in vacuo was measured by means 

3M. G. Holloway and B. L. Moore, Phys. Rev. 58, 847 
(1940). 


4M. Stanley Livingston and H. A. Bethe, Rev. Mod. 
Phys. 9, 245 (1937). 
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ENERGY 


of a sensitive galvanometer. The other was a 
standard counter of the type shown in Fig. 2. In 
order to use such a counter in this application, 
‘t was necessary to reduce the beam current to 
a low value, employ a very small aperture 
(0.025-cm diameter), and measure the output of 
the counter by means of a counting rate meter 
circuit which would record up to 50,000 pulses 
per second. The results obtained from the use of 
these two detectors were in good agreement ; the 
averaged data from six separate runs give 
3.68+0.05 Mev for the mean energy of the beam. 
When correction is made for the straggling of 
the deuterons resulting from their passage 
through the absorbing material, the half-width 
of the differential energy curve turns out to be 
0.18+0.04 Mev. 

A study was made to determine the effect on 
beam homogeneity of varying certain of the 
cyclotron parameters, namely, the RF voltage, 
deflector voltage, gas pressure, and cross section 
of beam used as determined by the position and 
separation of the collimator jaws at the entrance 
to the bombardment chamber. Although a some- 
what more homogeneous beam was indicated at 
low deflector plate voltages, no significant change 
in beam homogeneity was found for reasonable 
variations in these parameters. 

The number of deuterons striking the target 
during each proton counting run was measured 
by a current integrator of the neon discharge 
type. The unit beam which was used for each 
counting interval was 3X 10'° deuterons, which 
corresponds to a beam current of 0.25 micro- 
amperes integrated over a period of approxi- 
mately 50 seconds. 


PROTON GROUPS FROM Mn**(d,p)Mn** 


Figure 4 shows the results obtained over the 
range of proton energies from 3.0 to 9.0 Mev. 
The lower limit to the measurements is set by 
the range of the deuterons scattered from the 
target. The proton group ending at about 8 Mev 
is the end group for the reaction in question. On 
several occasions attempts were made to find a 
group of higher energy. The results were in all 
cases negative; it can be stated that if a group 
having an energy greater than 8 Mev exists, 





*B. E. Watt, Rev. Sci. Inst. 12, 362 (1941). 
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the yield is less than 5 percent of that of the end 
group indicated in Fig. 4. 

The less prominent groups having energies 
greater than 4.7 Mev are shown more clearly in 
Fig. 5, where the scales have been expanded. A 
total of 21 separate runs were made over this 
interval of proton range, ard all of these have 
been used for analytic purposes, although only 
11 were used to obtain the averaged points shown 
on this graph. In Fig. 5, the number of counts 
taken to determine a particular point is given 
by the ordinate in each case. 

The numerical results obtained from these 
data are shown in Table I. The proton range has 
in each case been corrected for the effects of the 
straggling of the protons in the absorbing ma- 
terial, the inhomogeneity of the deuteron energy, 
and the effective counter depth. The error in 
the measured proton range caused by finite 
geometry of the detection system was calculated 
to be 0.3 percent for the most unfavorable case. 
The error caused by the finite thickness of the 
target is also small: 0.1 cm in range for the 30-cm 
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Fic. 4. The yield of protons per unit deuteron beam is 
plotted against proton energy. Eight runs were made over 
the energy range from 3.0 to 4.7 Mev; 21 runs were made 
over the interval from 4.7 to 8.5 Mev. 
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Fic. 5. A plot of the proton yield as a function of proton 
range, showing the less prominent groups having energies 
greater than 4.7 Mev. The number of counts taken to 
determine a particular point is given by the ordinate in 
each case. 


proton group, for example. These last two sources 
of error have been neglected in computing the 
disintegration energies for these reactions. The 
reproducibility of the results obtained is indi- 
cated in Table II, where several values for the 
various disintegration energies, obtained from 
different runs and different sets of averaged data, 
are listed. An analysis of the experimental errors 
which enter into this method yields the figure 
+0.11 Mev for the reliability of the Q-values 
obtained, and a study of the half-widths of the 
proton groups observed indicates that the resolu- 
tion of this method was 0.30 Mev in the present 
case. 

In Fig. 5, a number of small deviations from a 

TABLE I. Numerical results for the proton groups from 


the reaction Mn**(d,p)Mn**. The mean proton range, after 
corrections have been applied, is shown here. 











Proton range Proton energy Q Levels Spacing 
(cm air) (Mev) (Mev) (Mev) Mev) 
22.9 3.98 0.38 4.38 
0.77 
30.8 4.73 1.15 3.61 
1.13 
44.7 5.85 2.28 2.48 
0.71 
54.4 6.55 2.99 1.77 
0.70 
64.7 7.23 3.69 1.07 
1.07 
82.3 8.29 4.76 0 


Average spacing 0.85 
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smooth curve, suggesting partially resolved 
groups, are noticeable, for example at ranges of 
about 34 cm, 41 cm, 44 cm, and 70 cm. None of 
these small irregularities survived a statistical 
analysis, however. Moreover, from a series of 
experiments that were performed to determine 
the effect of small amounts of the most probable 
impurities that might be present in the manga- 
nese metal or deposited on the surface of the 
target, it was found that approximately 49 
percent of the small deviations were to be 
associated with proton groups from (d,p) reac. 
tions in sulfur and in deuterium. It is concluded, 
therefore, that these small irregularities are not 
to be taken as evidence of partially resolved 
groups from the reaction Mn**(d,p) Mn*, 

The results of the present work are in general 
agreement with those of other investigators who 
have studied heavy particle emission reactions 
of this general type. The persistence of level 
spacings of the order of 0.85 Mev in the nucleus 
Mn* supports the view that the nuclear energy 
levels exhibited in the process of heavy particle 
emission reactions are restricted by some selec- 
tion principle, not yet fully understood, to a 
small fraction of the number that are theoreti- 
cally possible. The fact that the level spacing 
observed for Mn** is somewhat smaller than the 
spacing found in lighter elements might be 
regarded as at least qualitative agreement with 
the prediction of theory that the density of 
levels should increase at higher atomic weights. 


YIELD OF INDIVIDUAL GROUP AS A FUNCTION 
OF DEUTERON ENERGY 


The energy of the deuterons striking the target 
was changed by introducing aluminum foils in 
the path of the beam at the entrance to the 
bombardment chamber. The deuteron energies 
attained in this way were 3.68 Mev, 3.33 Mev, 
and 3.15 Mev. Measurements were then made 


TaBLE II. Several values for the disintegration energies, 
obtained from different runs and sets of averaged data. 











Excitation Q-values (Mev) Best value 
Ground 4.75, 4.74, 4.78, 4.79 4.76 
First 3.64, 3.72, 3.67, 3.74 3.69 
Second 2.96, 2.98, 3.02, 3.03 2.99 
Third 2.Bty Bchty Sede 250 2.28 
Fourth 1.14, 1.15, 1.16, 1.17 1.15 
Fifth 0.31, 0.38, 0.38, 0.43 0.38 
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to determine the relative yields from the four 
most prominent proton groups of the Mn* 
spectrum, namely, the first, second, fourth, and 
sixth, in order of increasing proton range. Data 
similar to those shown in Fig. 4 were obtained 
at each beam energy, and care was taken to 
define the peak of each group as accurately as 
possible. b ; ; 

The results, shown graphically in Fig. 6, 
suggest that broad resonance phenomena are in 
operation in the case of groups A and B, the 
maximum yield for these groups occurring for a 
beam energy of about 3.4 Mev. Groups C and 
D, however, show no evidence of a resonance 
process for beam energies between 3.15 and 3.68 
Mev. These results are in general agreement with 
the findings of Bennett, et al, who observed that 
the yield of a single proton group from the 
reaction C"(d,p)C™ was a slowly varying func- 
tion of the deuteron energy, with no indication 
of sharp resonance effects. 

The dashed curves in Fig. 6 are approximately 
the envelopes of the proton energy distributions 
for the two deuteron energies indicated, and 
correspond to the distribution of proton energies 
suggested by Volkoff’s’ treatment of the Oppen- 
heimer-Phillips process. It is of interest to note 
that the present results are not in accord with 
Volkoff’s development: (1) the amplitude of the 
distribution falls off more rapidly than £,~+ for 
the proton energies obtained in the present work, 
and (2) the distribution does not exhibit a 
maximum in the vicinity of 


, Yl Ss 
‘Ti 2(14+11.1/Z) 





|- 4.9 Mev. 


The possibility of a maximum in the distribution 
at a somewhat lower proton energy (<3.0 Mev) 
could not be investigated in the present work 
because it was not possible to count protons 
effectively in the presence of the comparatively 
large flux of deuterons scattered from the target. 


MASS OF THE NUCLEUS Mn 


If it is assumed that the emission of the proton 
group of longest range leaves the Mn* nucleus in 


the ground state, a knowledge of the Q-value for 


_*W. E. Bennett, T. W. Bonner, H. T. Richards, and B. 
E. Watt, Phys. Rev. 59, 781 (1941). 
7G. M. Volkoff, Phys. Rev. 57, 866 (1940). 
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Fic. 6. The relative yields of the four most prominent 
proton groups from the reaction Mn**(d,p)Mn* for three 
different deuteron energies. 


this group enables the calculation of the mass of 
the nucleus Mn*5, 

From Dempster’s® value for the packing frac- 
tion of Fe®*, the data of Elliot and Deutsch® for 
the beta- and gamma-rays emitted by Mn*, and 
the maximum Q-value obtained in the present 
work, the mass of Mn® is found to be 54.9634 
+.0022 mass units. This result agrees, within 
the limits of error, with the value found by 
Davidson" by the use of the same general method. 
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4, 
capture has been produced! by proton bom- 
bardment of Mo. Gamma-rays of energy 0.05 
Mev and 0.5 Mev and beta-rays of energy 0.6 
Mev were found to be associated with the ac- 
tivity. A Tc activity with a period of approxi- 


Radioactive Isotopes of Mo and Tc 
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SEPTEMBER 1, 19, 


Evidence is presented which permits the assignment of a 4.3-day Tc activity to Tc. It 
decays by K-capture. X-ray photographs of long-period Tc activities show Mo x-rays resulting 
from K-capture decay of a 52-day Tc activity. Tc x-rays resulting from an isomeric transition 


involving a 95-day Tc*’* activity were also photographed. 





I. 67-HOUR Mo” AND 4.3-DAY Tc 


mately 2 days was also reported.? It was produced 


6-DAY Tc activity which decays by K- by deuteron bombardment of Mo and strong 
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x-rays characteristic of Mo were found. 

A 67-hour activity assigned to Mo* has been 
produced by the reactions Mo(d, p)*, Mo(n, 38, 
Mo(n, 2n)**, Zr(a, m)*® and is found in fission 
products from U and Th.** The isotope, with 


Fic. 1. Decay curves 
from Mo-+d showing pro- 
cedure of x-ray photo 


graphy. 
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Fic. 2. Photographs of the combination of 4.3-day Tc® and 
67-hour Mo®., 


characteristic radiations of 1.5 and 1.2-Mev 
8-particles and gamma-rays of 0.4, 0.24, and 0.75 
Mev, decays to a 6-hour daughter product Tc®*. 
The 6-hour activity with the emission of 0.136 
and 0.18-Mev gamma-rays and conversion elec- 
trons drops to the ground state,® 4 x 10° year Tc®®. 
The K and LZ conversion electrons® and x-rays"® 
characteristic of Tc have been identified pho- 
tographically. 

43-day Tc and 67-hour Mo activities are 
both produced by deuteron bombardment of Mo. 
A study of the x-rays associated with the decay of 
these isotopes and long-lived Tc isotopes was 
undertaken using curved mica crystal spectro- 
graphs." A deuteron activated Mo sample was 


*D. C. Kalbfell, Phys. Rev. 54, 543 (1938). 

” P. H. Abelson, Phys. Rev. 56, 753 (1939). 

"J. E. Edwards. M. L. Pool, and F. C. Blake, Phys. 
Rev. 67, 150 (1945). 
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Fic. 3. Exposures of 11 days; chemical fractions from 
Mo-+d were used. Exposures started 10 hours after 
bombardment. 


photographed for a total of 422 days. Five suc- 
cessive photographs were made. The photographic 
procedure for the first 4 photographs is shown in 
Fig. 1. Photographs A and B (Fig. 2) were ob- 
tained while the 67-hour Mo" and the 4.3-day Tc 
were predominant in the activity. Very strong 
Mo K lines and a weak Tc K line are found in both 
photographs. The absence of any marked differ- 
ence in the relative intensities of the Mo Ka and 
Tc Ka lines indicates that the sources of the Mo 
and Tc x-rays have half-lives of the same order of 
magnitude. 

In order to identify the atomic number of the 
x-ray sources in photographs A and B, photo- 
graphs of the Mo and Tc chemical fractions from 
a deuteron-activated Mo sample were obtained 
(Fig. 3). The exposures were made in similar 
spectrographs for 11 days, starting 10 hours after 
bombardment. The Tc fraction showed strong 
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Fic. 4. Decay characteri ti 
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Mo lines while weak Mo and Tc x-rays were 
found in the Mo fraction. A small amount of Tc 
in the Mo fraction was responsible for the weak 
Mo Ka line. Additional exposures of the two 


MASS NUMBER 
93] 94] 95} 96] 97} 98 


fractions for 14 days gave no lines from the Mo 
fraction and a very weak Mo Ka line from the Tc 
fraction. These experiments demonstrate that 
Mo x-rays are associated with a Tc activity of a 


99} 100} (Ol 























Fic. 5. Section of the 


atomic chart. 


























67H | 9.8% 14.6M 
ate. 49 


rT 3) 


“KS 
































i} —_ 


3, 
4x10*Y ity 

















teristics 
om 4,3. 


1e Mo 
the Tc 
> that 
y ofa 





ISOTOPES OF 


$0,000}— 











10,000 
5.000 


$3 


_ 6. Decay characteristics 
™” of 67-hour Mo”. 


INTENSITY 








MO AND TC 387 


67-HOUR Mo% 
DECAY CURVES 


Mo+d (Mo Fraction) 
(B+y)- Decay 


Mo+d (Mo Fraction) 
(y only)- Decay 








100 
$0 
10} 
se 
! ! | | Lo 
0 Te) 20 30 40 50 


DAYS AFTER 


few days half-life and Tc x-rays are associated 
with an Mo activity also of a few days half-life. 

Decay characteristics (Fig. 4) and radiation 
energies of a 4.3-day Tc activity were obtained. 
The decay characteristics from Mo+p, Mo+d 
(Tc fraction), and Cb+e (Tc fraction) all show 
half-lives of nearly 4.3 days. This 4.3-day activity 
can thus be definitely assigned to Tc**. Long- 
period activities were also present in the Mo+p 
and the Mo+d (Tc fraction) samples but no 
activity longer than 4.3-days half-life was ob- 
served in the Cb+a (Tc fraction) sample. By Al 
and Pb absorption measurements 4.3-day Tc*® 
was found to emit x-rays characteristic of the Mo 
region, 0.64-Mev charged particles in agreement! 
with other work, and gamma-rays of energy 0.92 
Mev. The photographic evidence of Mo x-rays 
associated with this activity leaves no doubt that 
it decays mainly by K-electron capture (Fig. 5). 

The Tc x-rays (Fig. 3) associated with 67-hour 
Mo* are due to internal conversion of a gamma- 
ray as the 6-hour daughter product, Tc®*, drops 
to the ground state, 4X10° year Tc®*. Sixty- 
seven-hour Mo®’ was produced by deuteron bom- 
bardment of Mo and by alpha-bombardment of 
Zr (Fig. 6). A rise of 6-hour period was observed 
in the gamma-radiation of the Mo fraction im- 
mediately after chemical separation from Tc. 
Half-life values and radiation energy measure- 
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ments on Mo** and Tc** are in agreement with 
other work.*'" No activities of half-life longer 
than 67 hours were observed in Mo. 


II. LONG-PERIOD ACTIVITIES IN Tc 


Long-period Tc activities (90 days and 62 
days) have been produced by deuteron bom- 
bardment® "-"* of Mo. The 62-day activity emits 
x-rays characteristic?“ of Mo while radiations 
from the 90-day activity are thought to be Tc 
x-rays“ and conversion electrons.'* Two con- 
verted gamma-rays in the long half-life activities 
have energies of 0.087 Mev and 0.184 Mev. 
Evidence of the growth of the 90-day Tc activity 
from 2.8-day Ru*” has been reported,!*!’ per- 
mitting its assignment to Tc®’. 

X-rays associated with long half-life activities 
after deuteron activation of Mo were photo- 
graphed according to the procedure shown in 
Fig. 1. Exposures C and D (Fig. 7) were obtained 
with a mica crystal curved to a radius of 8 inches. 
Mo Ka and Tc Ka lines appear in both the 78- 
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on ees Te, Ka Mo KB day exposure and the 131-day exposure with ap. 
— fre proximately the same intensity. However, it is 
Photograph C quite evident from the original films that the 

Mo Ka line in exposure C is slightly stronger than 

78-day exposure starting the Tc Ka line, while in the later photograph, D 
ment the Tc Ka line is slightly stronger. A third mon 

posure of 195 days following exposure D showed a 

very faint Tc Ka line, but no Mo lines. From 

these three photographs it is obvious that the 

source of the Mo x-rays was decaying faster than 

Mo —" the source of the Tc x-rays. To explain these two 
x-rays arising from Tc activities one must assume 

two activities in Tc of slightly different period, 

The shorter period decays by K-capture giving 

Mo x-rays, while the longer period Tc x-rays 

must result from converted gamma-radiation in 

Photograph D an isomeric transition in Tc. The two activities 
cannot be genetically related since the Mo x-ray 

131-day exposure starting intensity dropped below the Tc x-ray intensity in 













148 days after bombard- ( 
_———" the later photographs and no positrons were § 
observed. ; ‘ 
Decay characteristics followed for 900 days \ 
show two Tc activities fairly well resolved having ‘ 
half-lives of 9545 days and 52+3 days (Fig. 8). | 
—_ The activities were produced by deuteron bom- t 
Mo calibration 
K lines bardment of Mo and by proton bombardment of é 
Mo. Radiation energies from these two activities ‘ 
obtained by Pb and Al absorption measurements " 
Fic. 7. Photographs of long-period Tc activities. show two groups of soft electrons (Fig. 9), x-rays 
] 
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Fic. 8. Decay characteristics of 
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characteristic of the Mo-Tc region, and two 
gamma-rays of energy 0.84 Mev and 0.25 Mev. 
The soft electron groups are no doubt the con- 
version electrons previously reported.!® Since the 
strong Tc x-rays following the longer 95-day half- 
life arise from an internally converted gamma-ray, 
the more intense 87-kv electrons should be 
associated with this activity. The evidence pre- 
sented here requires the assignment of the 95-day 
Tc activity to a metastable state in Tc (Fig. 5). 
In the isomeric transition which occurs a 0.108- 
Mev gamma-ray is strongly internally converted 
giving rise to the Tc x-rays. This interpretation 
requires the existence of another stable or very 


20 40 60 : 80 100 
MILLIGRAMS PER CM* Al 


long-lived isotope in Tc. Such a long-lived isotope 
could logically be placed at Tc*’ permitting the 
assignment of the 95-day isotope to Tc*’* in 
agreement with previous work."” 

Observations on the gamma-ray decay charac- 
teristics indicate that the hard gamma-rays 
follow the 52-day half-life. This activity decays 
by K-capture with associated gamma-rays of 0.84 
Mev and 0.25 Mev and is tentatively assigned to 
Tc® (Fig. 5). 

The authors are grateful for the support re- 
ceived from Mr. Julius F. Stone, the Ohio State 
University Development Fund, and the Ohio 
State University Research Foundation. 











































PHYSICAL REVIEW 


VOLUME 72, 





NUMBER 5 SEPTEMBER 1 


1947 


A Class of Exact Solutions of Einstein’s Field Equations 


SupHANsU Datta MAJUMDAR 
Palit Laboratory of Physics, University College of Science, Calcutta, India 
(Received March 14, 1947) 


The work of Weyl on the gravitational field occasioned by an axially symmetric distribution 
of matter and charge is generalized to the case in which gq, and ¢ for an electrostatic field are 
functionally related, with or without spatial symmetry. It is shown that the most general 
electrostatic field in which gq and ¢ are related by an equation of the form ua =4(6+c) can 
be represented by a line element of the form (ds)*=—e[(dx')?+ (dx*)?+ (dx*)*]+e(dt)?. 
Certain of the field equations are then identically satisfied while the remaining ones reduce to 
a single equation for w. The substitution w= —2 log(1+v) transforms this into Laplace's 
equation for v, so that the solution can be expressed in terms of harmonic function. 





1, INTRODUCTION 


_ a well-known paper H. Weyl! showed that 
the gravitational field in empty space caused 
by any static axially symmetric distribution of 
matter and charge can be represented by the line 
element 


(ds)? = —e*[(dz)?+ (dr)? ] 
—r'e—”(d60)?+e(dt)? (1) 


with a suitable choice of coordinates. Using these 
coordinates, which he calls “canonical coor- 
dinates,” Weyl completely solved the problem 
of the pure gravitational field with axial sym- 
metry and also obtained a particular class of 
solutions for an axially symmetric electrostatic 
field; namely, those involving a functional rela- 
tion between the electrostatic potential, ¢, and 
the component g«, of the metric tensor. 

In this paper the more general problem of any 
electrostatic field (with or without spatial sym- 
metry) in which ga and ¢ are functionally re- 
lated, has been considered. It is found that the 
only type of functional relation that can exist 
between gas and ¢ is of the form 


gu=A+Bo+}¢’, (2) 


where A and B are arbitrary constants and 
¢/(8x)* is the electrostatic potential in rela- 
tivistic Lorentz units. If now the constants, A, 
B, be so chosen that the right-hand side of (2) 
becomes a perfect square, then the relation (2) 
reduces to 


gu =3(O+C)*. (3) 


1H. Weyl, Ann. d. Physik [4] 54, 117 (1917). 
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It is then shown that the most general electro- 
static field in which ga and ¢ are connected by 
the relation (3) can be represented by the line 
element 


(ds)? = —e~*[ (dx)? + (dx*)?+ (dx*)*]+e(dt)? (4) 


with a suitable choice of coordinates. Some of 
the field equations are then identically satisfied, 
and the remaining ones reduce to a single equa- 
tion for determining w. The substitution w= 
—2 log(i+v) transforms this differential equa- 
tion into Laplace’s equation for the function ». 
The solution can, therefore, be expressed in 
terms of harmonic functions. 

The solution thus obtained is, in a sense, more 
general than that of Weyl, since in this case there 
need not be spatial symmetry of the field. In 
another sense, however, it is less general. For, the 
functional relation (3) adopted here is a special 
case of the most general relation (2) adopted by 
Wey] in the case of axial symmetry. 

The solution outlined above is worked out in 
details in Sections 4 to 6. In Section 5 it has 
been pointed out that the nature of the solution 
is such that if the field be due to several charged 
bodies separated by empty space then to the 
Newtonian approximation they remain in equi- 
librium because of their mutual gravitational 
attraction and electrostatic repulsion. 

In Section 7 Weyl’s axially symmetric solu- 
tions have been examined, and it has been shown 
that in certain cases no solution of the field 
equations for empty space exists. A necessary 
and sufficient condition for the existence of static 
axially symmetric solutions has also been ob- 
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tained. It has been proved that static field is 

ible only when the masses are in equilibrium 
under the influence of one another. This brings 
to relief a distinct superiority of the relativistic 
theory over the older Newtonian theory. 

In the last section an attempt has been made 
to bring out the strong analogy between the old 
and the new theories by establishing a few 
simple theorems on the uniqueness of solutions 
and dependence of solutions on boundary values. 
The general boundary value problem for any 
static field is, however, far from being solved. 


2. THE FUNDAMENTAL EQUATIONS 


As usual we represent the static field by the 
line element 


(ds)? = gaxdx*dx? + g44(dt)?, 


the coordinate system so chosen that gu=go 
=gs,=0, and the metric tensor gj is independent 
of t. Since no magnetic field is present, the four- 
potential has the components (0, 0, 0, ¢/ (82)#), 
where ¢/(87)! is the electrostatic potential in 
relativistic Lorentz units. The only surviving 
components of the antisymmetrical tensor, Fi, 
and its associate F*/, representing the electro- 
magnetic field-strengths, are, therefore, 


(a, b=1,2,3) (5) 


Op 

(8x) 'Fys=—=¢,, 
Ox* 

(8x) Fti=e—’g*g,, (i=1, 2, 3) 


where e” has been written for gas. From these the 
components of the electromagnetic energy tensor 
Ej are calculated by means of the formula 


Eji=— FF ig+igiF* Fo. 


The gravitational tensor G,; is calculated by the 
formula 


0 
Gij= ——{1j, a} + {ta, b} { 7d, a} 
ax* 


oe? 





a 
+ log(—g)!— {tj, a} —— log(—g)!. 
axe 


Ox'dx 


The condition for the absence of electric 
charge and current then comes out to be 


Fi) 
—(e~"g*g»( —g)*) =0, (6) 
axe 
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which corresponds to the classical equation 
V*o =0. 
In addition we have the seven field equations 


Gij= —8rE;; (i, j= 1, 2, 3) (7) 
a= —8rEu, (8) 
the remaining three equations G4 = —8rE, for 


i=1, 2, 3 having been identically satisfied. The 
eight equations, (6)-(8) form the set of funda- 
mental equations of the electrostatic field. 


3. SOLUTIONS INVOLVING A FUNCTIONAL 
RELATION BETWEEN gu AND ¢ 


In this paper we shall be mainly concerned 
with those solutions of the field equations in 
which ga and ¢ are functionally related. First 
we shall show that the only type of functional 
relation that can exist between them is of the 
form 


£4=A+Bo+}¢’, 


where A and B are arbitrary constants. In order 
to prove this we shall require only the two 
Eqs. (6) and (8). 


re] 
Gu= ——{44, a} +2{4a, 4} {44, a} 
Ox* 
0 
— {44, a}— log(—g)! 
Ox 
1 @ 
=-—(g fr) —3-f-' fag fr 
2 dx 


1 a 
+—g" fr: — log(—g)!, (9) 
2 ax 


where f is written for ga. Also 
— 8rE = }2" bad. (10) 


Therefore, by Eq. (8) 
re) 
—(g" fr) —f-* a” fafe 
Ox* 
] 
+2” fr-— log(—g)*=g*badr. 
Ox* 


Since f is a function of ¢, this equation can be 
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written as 

0 
—(g*( —g) *¢s) 
ox* 


=(f)-M1+ SF? —f") -8°°(— 8) ado, (11) 


where f’=df/dd and f’ =d*f/d¢*?. Equation (6) 
can be written as 


0 
—(g(—g)*- os) = f-'f’ -g(—g) babs. (12) 
axe 


Now Eggs. (11) and (12) can hold simultaneously 
only if 


(i+f-f?—f)/f =f /f, 
that is, if 1—f’’=0. From this it follows that 
f=gu=A+Bot+i¢ 


is the only functional relation possible between f 
and ¢. 

Since any metric approaches that of the special 
relativity theory at great distances from matter 
and charge, gis must tend to unity at infinity 
with the proper choice of the unit of time. And 
if we choose the arbitrary constant in ¢ so as to 
make ¢ vanish at infinity, then the above relation 
may be written in the form 


£u=1+Bo+}¢". 


From general considerations, solutions of this 
type are expected to be a bit easier to obtain, 
since the relation (13) makes the Eqs. (6) and 
(8) identical, so that there is one equation less 
to be satisfied. The utmost simplification is 
expected to result from taking the right-hand 
side of (13) to be a perfect square. The relation 
(13) then reduces to 


gs =3(6+V2)?. 


(13) 


(14) 


4. A CLASS OF EXACT SOLUTIONS 


It has been possible to find the most general 
solution of (6)—(8) involving the relation (14). 
The relation (14) itself, however, restricts the 
generality of the solutions to a great extent. 

We proceed to obtain the expressions for the 
various tensors in terms of g;;, w, and @ with the 
help of the formulas given in Section 2. 
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— 82E;;=e~"$.b;— 3-8 8" bade, 


(i, j=1, 2, 3), (15) 


—_ 8rEuy = 32 °badr, (16) 


1 a 
Gas =— —(g%e"wr) — Fe’g™ wows 
2 ax? 


f) 
+ 3e"g"ws— log(—g)!. (17) 
axe 


On account of the relation (14), Eqs. (15), (16), 
and (6) can be written as 


—8rE;;= 4wwj— 42 ig Was, 


(7, j=1, 2,3) (18) 


—8rE 4, = te"g wows, (19) 


ts] 
—(e~**g*wy( —g)!) =0. (20) 
Ox* 

From (8), (17), (19), and (20) we find that 
the field equations, (6) and (8), become identical 
on account of the relation (14). 

Next we shall show that the relation (14) 
makes it possible to reduce the general static 
line element (5) to the isotropic form 


ds? = —e~[ (dx)? + (dx*)?+ (dx*)? ]+-e"(dt)? (21) 


by suitable transformation of coordinates. In 
order to do this, we introduce two sets of func- 
tions *g;; and 9g;; (t, 7=1, 2, 3) defined by the 
equations 


2i5= —*g5= —E-"9i;, (4, J=1, 2, 3) 


whence 


*9:;is the metric tensor of the three-space in (5) 
and @,; that of an associated three-space con- 
formal to it. The contracted Riemann tensors 
*G,;; and Gi; of these two spaces in conformal 
correspondence are connected by the equations 


*Gis= Gij— $0, 5 — ww; 
+9: —}4.w+tAw], 
Aww => GD, ab 


(22) 


where 
A,w = 9 war, 


and W,,; is the second co-variant derivative of 
the scalar function w in the associated space. 








(15) 
(16) 


(18) 
(19) 


(20) 


d that 
entical 


n (14) 
static 


(21) 


es. In 
’ func- 
»y the 


in (5) 
€ con- 
ensors 
formal 
ations 


(22) 


ive of 
ce. 








The relation between the contracted Riemann 
tensors Gis (¢,7=1, 2,3) and *G;; of the four- 
space and three-space in (5) is 


Gij=*Gist }wist Pwwj—F* (tj, a} We 
=*G,,+3*w,;+iwm;. (t,j7=1,2,3) (23) 
Further we have the relations* 
*{ij, k} = [aj, Rj 
+3[ — wb! —wib;§ + WoG ig] 

*{ij, a}wa=~ (47, a}Wa—wwjt 39 Aw. (24) 
From Eqs. (22)—(24) we have the relations 
Gig=Gist}wwj—PGisAw. (t,7=1, 2,3). (25) 
Equation (18) can also be written in the form 
~8xE,j=}wwj—}gijAw. (i, 7=1,2, 3) (26) 


By the use of the expressions (25) and (26) for 
Gi; and E;;, the field equation (7) can now be 
written as 


Gi;—439:([4ew—}Aw)]=0. (7, 7=1,2,3) (27) 
The expression (17) for G4, can also be written 
in a very compact form 
1 0 
Gu=-- —(e?”g>w) + de?"Aw 
2 dx* 
‘ 
re) 
— fe” g"u— log(—g)! 
Ox? 


1 age 
om al gwar +——W 
2 ax* 


0 
Pia tog(a)'|. (28) 
ox? 


By use of the relations 
agi 
dx" 





= — fak, i} gr— fak, j} 9%, 


0 
— log(g)*= © {1a, a}, 
ox' 


* Editor’s note: for typographical reasons the Christoffel 
symbols associated with the metric coefficients *g;; and 


9ij are indicated with the superposed signs * and “. 
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Eq. (28) can be written as 
Gus =— he?"A.”. (29) 


Because of (19) and (29) the field equation (8) 
now takes the form 


—fe’A.w=—te”Aw, or Asw=4A,w. (30) 


From (27) and (30) follows the interesting 
result 


G,;;=0. (4, 7=1, 2, 3). (31) 


Since in a three-dimensional space the vanishing 
of the contracted Riemann tensor implies the 
vanishing of the Riemann tensor itself, we arrive 
at the conclusion that the associated three-space 
is Euclidean. We are, therefore, justified in 
adopting (21) as the line element of the most 
general static field in which the relation (14) 
holds. According to (30) the field equations will 
then reduce to the single equation 


Vw = }(wr+w.?+w;’). (32) 
The substitution, 
w= —2 log(1+2), (33) 
transforms it into Laplace’s equation 
Vv=0 (34) 


for the function v. 
Also according to (14) and (33) the relation 
between v and ¢ is 


v= —$¢/(¢+V2). (35) 


5. NATURE OF THE SOLUTION OBTAINED 


The Newtonian approximation of the above 
solutions is obtained by neglecting all powers of 
v higher than the first. We then have 


2a =e” =1—2v=1—29, 


where © is the gravitational potential in gravi- 
tational units. Also 


¢/(8x)'= +0/2(x)!. 


The potential ¢’ in electrostatic units is, there- 
fore, 


whence 


Q=+¢’. (36) 
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From this we see that if the field be due to 
several charged bodies separated by empty 
space, then to the Newtonian approximation 
they must remain in equilibrium under their 
mutual interaction. In fact it is because of this 
interesting property of the solutions that a 
reduction of the field equations to Laplace’s 
equation has been possible. 


6. VANISHING OF THE MATERIAL STRESSES FOR 
A PROPER CHOICE OF THE INTERNAL 
FIELD 


In the Newtonian theory if the densities of 
charge and matter bear the same constant pro- 
portion everywhere, and if in the above-men- 
tioned units the constant of proportionality is 
unity, then owing to the exact balancing of the 
gravitational and electric forces on every piece 
of matter, the stresses in matter will vanish. It is 
remarkable that such a result exactly holds in 
the relativistic theory also. Of course, given an 
external field in empty space the distribution of 
matter and charge producing that field is far 
from being uniquely determined. We shall show 
below that, corresponding to the external solu- 
tion obtained above, a particular internal solu- 
tion can always be constructed so that the 
stresses in matter may vanish at every point. 
We choose the internal solution in the following 
manner. 

In the interior of matter: (i) there exists a 
line element of the form (21), (ii) v is a solution 
of Poisson’s equation, (iii) gas and @ are con- 
nected by the relation (14), (iv) the first, second 
(and whatever derivative may be required), of w 
are continuous at the boundary surface, sepa- 
rating matter from empty space. 

The internal field should satisfy the equation 


Gi—4}giG = —8x(Mi+E,), (37) 


where M,/ is the material and £;/ the electro- 
magnetic energy tensor. Since Eqs. (25) and (29) 
hold for any static field, and Eqs. (19) and (26) 
follow from the assumption (iii) alone, they 
continue to hold in this case also. On account of 


the assumption (i) we can put 
Gij=0, 9j=85;, (i, 7=1, 2, 3) 


in (25), (26), (29), and (19) and obtain the ex- 
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pressions 
Gi;= 4ww;— }6,;V*7w _ 
— 8x Ej;= }ww;— 36,2," ; Stee 
Gu= —}e’V*w, 
G= —hevDw,?+e"Vw, 
Gi—}g/G= —}e°wwj+}6,e°2w,?, 
(7, j=1, 2, 3) (38) 
—8rEi = — he"ww;+ }e"s;,2w,?. 
(7, j=1, 2, 3) (39) 
Combining Eqs. (37)-(39) we arrive at the 
interesting result 
M,=0, (4, 7=1, 2, 3), 


that is, the material stresses all vanish at every 
point. But the situation is quite different in the 
case of the component M,‘ which represents the 


density of matter. 
By (19), (29) 


Gat — }gs'G = —e"V*w-+ je"=w,’, 
—8rE,‘= —je"=w,’. 
Because of Eq. (37) we, therefore, have 


—8rM = —e"[V*w—42w,7 |, 
whence 
V2v = —4n(1+0)*§M,'. 


For small v this reduces to Poisson’s equation. 
In empty space M,‘=0, and we get back 
Laplace’s Eq. (34). 


7. EXAMINATION OF WEYL’S AXIALLY 
SYMMETRIC SOLUTIONS 
(A) The Case of a Pure Gravitational Field 


Weyl has shown that the most general static 
axially symmetric gravitational field in empty 
space can be represented by the line element 


(ds)? = —e“[ (dz)?+ (dr)? ] 
—re—”(d0)?+e"(dt)?. (40) 
Weyl chooses the internal line element also in 


this form. We then have the following expres- 
sions for the components of the material energy 
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tensor : 
8x7 '= —82T?? 
= }e-"(2u2/r+2we/r+wi?—wy*), (41) 
89733 = feo" (Mu tun +Wu + Wet pwr + w,’), 
—8aT f= — 37" (Wir t+ Wee — 11 — Mee 
+2w-/r—4w,?—}w-"), 
—8x7}2= —8xT2! = — he" (wiw2—wi/r —u;/r). 


(42) 


Weyl makes the assumption that the stress 
component 7;* vanishes everywhere. This leads 
to the two differential equations, 


Aw=Wi1t+Wo2tWe/r = 82T ye", (43) 


V=011 +022 = —3 (wi? +w2”), (v=utw) (44) 


for determining u and w. In empty space T,‘=0, 
and Eq. (43) reduces to Laplace’s equation for 
axial symmetry. This equation is first solved for 
w and the solution substituted in the right-hand 
side of (44). The latter equation can then be 
treated as a two-dimensional Poisson’s equation 
for the function v. The solution for v is made 
unique by the condition that it vanishes on the 
axis and at infinity. Thus we see that the two 
assumptions, 7;*=0 and 7,‘=0, made by Weyl 
enable us to find any desired solutions for u and 
w which, when substituted in (41) and (42), will 
determine the stress components 7;', T;*, and 
T? uniquely. In order that the solutions thus 
obtained may be the correct solutions for empty 
space, it is essential that 7;', 7;*, and 7,2 should 
vanish in addition to T;* and 7,4. Weyl’s 
method of solution, in many cases, does not lead 
to the vanishing of these components. In order, 
therefore, to arrive at the correct solutions for 
empty space, we set the left-hand sides of Eqs. 
(41) to (42) equal to zero. The functions w and v 
then satisfy the equations 


Aw =0, (45) 
V1 =rWiWe, (46) 
V2 = $r(w2?—w,’). (47) 


As before, we first integrate Eq. (45) by intro- 
ducing an arbitrary axially symmetric density 
distribution ¢. This solution for w is then sub- 
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stituted in Eqs. (46) and (47). If a solution for 
v can now be found satisfying Eqs. (46) and 
(47), then Weyl’s Eq. (44) will be automatically 
satisfied in the empty regions of space. But the 
converse is not true. A necessary and sufficient 
condition that a function v may exist satisfying 
(46) and (47) is that the line integral 


I= f (v\dz+vedr). (48) 


(i) along any closed curve, A, which lies entirely 
in empty space should vanish, and (ii) along any 
curve, B, which starts from the axis and ter- 
minates on the axis and lies entirely in empty 
space, should vanish. 

First we consider curves of the type A which 
may surround matter distributed in the shape of 
a ring. By Stokes’ theorem, 


hes f f u—eddS, 


where dS is an element of plane area in the 
(r—z) half-plane enclosed by the curve, and 


Vi2= dv; /dr, V1 => Ov2/dz. 


1, V2 are looked upon as two different functions 
and not as the derivatives of a single function v. 
By (46) and (47) 


Vi2—V21 = 7W\° Aw. 


If the enclosed area is empty then Aw=0 and, 
therefore, J vanishes. If, however, the area. 
contains matter, then 


Aw = —4ro, 


ey se 
= -2f {ff resedsirao 
a 


the last integral having been taken over the 
entire ring. Now w, is the component of force in 
the z direction acting on the mass element edr, 
and w; consists of two parts, 


and 


W1 = Wy, itW. ly 
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where w;, is the internal force attributable to 
the remaining part of the ring and w,, is the 
resultant of external forces acting on the element 
odr. The contribution to the integral J due to 
w;,, therefore vanishes, and we have 


r=-2f f fr sodr. (49) 


Thus we obtain as the necessary and sufficient 
condition for the existence of v that the “‘re- 
sultant force acting on the ring according to 
classical calculations must vanish.” 

Next we consider a curve of the type B which 
starts from a point P on the axis, terminates at 
a point Q on the axis, and lies entirely in empty 
space. After going from P to Q along the curve 
if we go back to P along the axis we complete a 
circuit, and the above proof can be easily ex- 
tended to this case also. 

From this we conclude that solutions corre- 
sponding to two spheres separated by empty 
space do not exist, whereas solutions correspond- 
ing to a sphere at the center of a circular ring do 
exist. 


(B) The Case of an Electrostatic Field 


Weyl] has shown that the most general axially 
symmetric electrostatic field in empty space can 
be represented by the line element (40). He 
takes the internal line element also in this form 
and obtains a class of solutions of the type (13). 
As his solutions may not lead to the vanishing of 
all the components of stress in empty space, we 
proceed, as before, to investigate the conditions 
for existence of solutions for strictly empty space. 
The field equations in this case are 


8rE|! = —82rE,? = fe” ($1? — $2”) 
= }e-“(2v2/r+w1?—w-’), 
—8nE;* =e" ($17 + $2”) 
= pe-“(V0+ wi? +}w,"), 
—8rEy= —}e~“” (o1 + 2”) 
= —}e-"(2Aw —V*v— }w,? —3w-?), 


—8rE,?*= —8rE,' = —e-"“did2 


= —4e—"(wywe—v1/r). 
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The equations to be satisfied by w and ¢ are 


therefore, 
Aw =e~"($1?+¢,"), (50) 
Ad = W1d1+ Wedr. (51) 


We shall consider Weyl’s solutions only, which 
are of the type (13). The two equations, (50) and 
(51), then become identical. Weyl introduces an 
auxiliary function x defined by the equation 


x= [1+ B6+4)-Ws. (52 


We then have 
Ax =e—“ [Ag — wid1 — Wo | 
= (B+¢)~'[Aw—e~"o:?—e-"g,"]._ (53) 
Hence (50) and (51) reduce to Laplace’s equation 
Ax =0. 


Inside matter we take the auxiliary function x 
to be a solution of Poisson’s equation. The equa- 
tions to be satisfied by v are 


V1 =TWW2— 2e~“’ pide, (54) 
V2 = $r(wWe* —wy”) +re-"(g1"— 2"). (55) 
In this case, 


Vi2—V.. =rwiAw+re—[ wid? — wide” 


=r[w:(B+¢) —2¢1 JAx by (53) 
=r(B?—2)x1-Ax by (52). 


Therefore, the line integral 


r= -2(B*=2) f ff x..0dr. (56) 


We therefore arrive at the same condition for 
existence of solution as before. From (56) we 
draw the interesting conclusion that if B=+v?, 
then a solution always exists irrespective of the 
distribution of the masses. This is in agreement 
with what has been said in Section (5). 


8. BOUNDARY VALUE PROBLEMS 


If gas and ¢ are functionally related, then every 
level-surface of gas (thereby meaning a surface 
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on which gas is constant) is also a level-surface 
of ¢. We shall now prove the proposition that if 
there exists one surface S (either closed or 
extending to infinity) on which ga, and @ are 
both constant, and if one of the two domains into 
which this surface divides the entire space be 
completely free from matter and charge, then 
in this domain, D, every level-surface of ga, will 
also be a level-surface of ¢, and therefore g4, and 
¢ will be connected by an equation of the type 
(13). 

To prove this we shall require only Eqs. (6) 
and (8). Equation (8) can be written as 


a 
—{[f-"(—g)!e fe ]=f-"(—g) *g"* bags. (57) 
axe 


Multiplying Eq. (6) by ¢ and subtracting from 
Eq. (57), we have 


a a 
—{f-'(—g)!g™f. ]-—@—Lf-"( —g) #2 oe ] 
ax* axe 
— f-'(—g) igo.» =0 
or 


a 
—[f-(—8)'e"(f—2¢*)»)]=0. = (58) 
axe 


Equations (6) and (58) are of the form 


a ny +b'u a =0, (59) 
where 
ait= — f-1(—g) git, (60) 
. 0 » 
Pe mt fee"), (61) 
Ox* 


u=o—A, or f—A, or f—}¢?—Be-1, 


Uji =OUu/Ox*, Urpin) =u /dx'dx*. 


We now require the help of a theorem given 
by Courant and Hilbert* on the uniqueness of 
the solutions of linear elliptic differential equa- 
tions. The proof, as given by these authors, is for 
a finite domain only, but the following modi- 
fication holds for both finite and infinite domains. 

We make the substitution 


u=zv, z=1-—e"", uw>O0 


*R. Courant and D. Hilbert, Methoden der Mathe- 
oe Physik (Julius Springer, Berlin, 1937), pp. 274- 


EINSTEIN’S FIELD EQUATIONS 








with 
r=L(x!)?+ (x)? (x4)? = [8,yeter}, 


The origin of space-coordinates is chosen to be 
outside the domain D, so that r is everywhere 
positive in D. Equation (59) is then transformed 
into 


av (inj) +B) +ov=0, (62) 
where 
Bi =22-a‘*g 14) +5! 
= 22-"pe"(a'*5,;x7/r)+b', (63) 


c =2- (a*21 54) +b'2;;)) 





O"* 55 5 emX?X™ 
= -riem| po -— 
r? 
On 5 jb emX?x™ bb 3x? 
=o (=-— ain |} (64) 
r r8 r 


Now we choose the space coordinates in such 
a way that the g,,’s approach Galilean values at 
infinity. We then notice that: (1) the a‘*’s tend 
to 6“ at infinity, and (2) 5‘, B*, c all tend to zero 
at infinity. Concerning the sign of c, which is very 
important, we see first of all that the a*’s form 
the matrix of a positive definite quadratic form. 
Therefore, the quantity a**6;;d:nx’x"/r? in the 
first term of (64) is everywhere positive in D and 
approaches the value unity at infinity, and 
hence possesses a positive lower bound. Similarly 
the coefficient of u in the bracketed expression 
in (64) tends to zero at infinity, and, therefore, 
its absolute value possesses a positive upper 
bound. Hence the quantity » can be chosen so 
large as to make this expression positive through- 
out D. Since z is positive throughout D, it follows 
that c is negative throughout D. 

We have thus established the two main 
requirements for the proof, namely, that a* 
should form the matrix of a positive definite 
quadratic form and that c should be negative. 
The proof now follows the same course as that 
given by Courant and Hilbert,* and we arrive 
at the conclusion that if v (and therefore u) 
vanishes on the boundary of D (and at infinity 
if the domain be infinite) then v (and therefore 1) 
vanishes throughout D. 
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DEDUCTIONS FROM THE ABOVE THEOREM 


(1) First we apply the theorem to the function 
u=f—}¢?—Bo-1. 

If there exists a surface S on which f and ¢ are 
both constant, then we can always choose the 
constant, B, so as to make u vanish on S. Hence 
it follows that if one of the two domains into 
which the entire space is divided by this surface 
be empty, then in this domain u vanishes and, 
therefore, f and ¢ are functionally related. That 
is, f and ¢ have the same family of level-surfaces 
in D. This is the theorem stated at the beginning 
of this section. 


(2) Applying the theorem to the function 
u=@—A we draw similar conclusions. 

(3) In the absence of an electric field we apply 
the theorem to the function u=f—A and draw 
similar conclusions. In this case it can be easily 
shown that space-time is flat in D. 

My thanks are due to Professor M. N. Saha 
D.Sc., F.R.S., for his constant interest and 
encouragement during the progress of the 
work. 

Thanks are also due to the authorities of the 
Council of Scientific and Industrial Research, 
India, for granting me a scholarship to Carry on 
this work. 
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A detailed discussion of the energy loss of negative mesotrons in matter is presented. The 
energy range considered is from +2000 ev to the lowest quantized orbit of the mesotron. The 
most important mechanism for energy loss is that of electron collisions except very near the 


nucleus, where radiation losses are important. 


The time for the over-all process is of the order of 10-" sec. in condensed matter and 10~° sec. 
in normal air. In chemical compounds the probability of capture near the various atoms is 


roughly proportional to their atomic numbers. 





1. INTRODUCTION 


ECENTLY the significance of experimental 

results! on the capture of negative mesotrons 
in matter has been discussed from the point of 
view of the information that it gives concerning 
the interaction between mesotrons and nucleons. 
The interaction of slow negative mesotrons with 
matter has been described as consisting of two 
steps:’ first, one in which the mesotron is cap- 
tured in the Bohr orbit with a radius of the order 
of 10-" cm near the nucleus; second, a step in 
which more typically nuclear interactions play a 
role during which the mesotron is destroyed by 
its collisions with the nearby nucleons. The 
present paper will be primarily concerned with 
the detailed description of the first step. 

The chief purpose of a detailed description of 
the capture process is to make sure that the 
time required for it is short compared with the 
natural decay time of the mesotron (~2X10-* 
sec.). We propose to discuss in particular how 
the physical and chemical state of matter in- 
fluences the capture process. In this connection 
we shall investigate also the relative probabilities 
of the capture of mesotrons near various types 
of nuclei in case the slowing down material is 
not a pure element. 

Throughout the greater part of the capture 
process the wave-length of the mesotron is short as 
compared to the geometric dimensions of the field 
in which this particle is moving. It is therefore 
permissible in most of our arguments to consider 


the motion of the mesotron as purely classical. 


1M. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 
71, 209 (1947). 

*E. Fermi, E. Teller, and V. Weisskopf, Phys. Rev. 71, 
314 (1947). 


2. ENERGY LOSS OF ELECTRONS OVER 2000 EV 


As long as the energy of the mesotron is more 
than 2000 ev, the velocity of the mesotron is 
greater than the velocity of the valence electrons. 
The slowing down of the mesotrons can then be 
treated according to the conventional methods 
applicable to fast heavy particles. In the slowing 
down of the mesotron the longest time is spent 
in the state when the mesotron moves with 
relativistic velocities. The consequences of the 
decay of the mesotron during this phase of cosmic 
radiation are well known and will not be dis- 
cussed here. The time needed to slow a mesotron 
from the relativistic 10° volts to 2000 ev is 
about 10~* to 10-?° second in condensed matter, 
or 1000 times as long in air. This part of the 
slowing process is again not our primary concern. 
It will be found that the time involved is con- 
siderably longer than the time needed for the 
later parts of the capture process. The proba- 
bility of spontaneous decay during this phase, 
which corresponds to a range in condensed 
matter of a few centimeters, is only of the order 
of 10-*. Consequently a negligible fraction of the 
decays observed with the ordinary experimental 
arrangements can be attributed to it.* Most of 
this time of 10-!° second is spent in the phases 
when the mesotrons still have velocities close 
to the light velocity. Actually the formula for 

dW 4ne*NZ 
Here W is the energy of the mesotron, V its 
* Decays during this phase are completely eliminated if 


energy loss per unit time is 
Onc 
=——— In ( ; (1) 
dt mV bmin 
the observations are carried out by delayed coincidence. 
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velocity, m the electron mass, N is the number of 
atoms with atomic number Z per cubic centi- 
meter; Dmax and bwin are the extreme values of 
the collision parameters. The logarithmic factor 
decreases during this interval of time from a 
value of the order 10 to zero. The latter value is 
reached when the mesotron velocity becomes 
equal to the velocity of an electron. Thus the 
contributions of successive electrons vanish. 
Finally, the velocity of the mesotron drops 
below that of the valence electrons. From then 
on we shall discuss the process of further energy 
loss in detail. 


3. LOSS OF ENERGY TO A DEGENERATE 
ELECTRON GAS 


When the energy of the mesotron has dropped 
below 2000 ev, and its velocity is therefore less 
than the velocity of the valence electrons, for- 
mula (1) no longer represents a useful approxi- 
mation, and the loss of energy to electrons can 
better be approximated in the following way. 

We consider the mesotron moving inside a 
degenerate electron gas with a velocity V much 
smaller than the maximum velocity vo of the 
electrons. In this case we can estimate the 
energy loss as follows: in an individual collision 
between an electron and the mesotron, the 
change in speed of the electron will be of the 
order of magnitude V. Indeed for a head-on 
collision it would be 2V. Since the electrons 
belong to a degenerate gas, it is clear that all 
the collisions for which the final velocity of the 
electron lies inside the occupied zone of the 
velocity space will be forbidden on account of 
the Pauli principle. Only electrons with speeds 
close to v9 by amounts of the order of V will, 
therefore, be capable of colliding. Their number 
per unit volume is of the order of magnitude 


n= mp? V /h'. (2) 
The collision cross section o for collisions in 


which the deflection is appreciable is, on the 
other hand, of the order of magnitude: 


o = (e?/mv,?)?. (3) 
The energy transferred in collisions of this type 


not forbidden by the Pauli principle will always 
be positive and of the order of magnitude 


W=mvoV. : (4) 
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From (2, 3, and 4) we can calculate the order of 
magnitude of the energy loss per unit time: 


—(dW/dt) = Wonvy = me V2/hs 
= m*eT / (uh*) a T/to, (5) 


where T is the kinetic energy of the mesotron 
to = ph® /m*e = 4.84 X 10-5 sec., and y is the main 
of the mesotron. We have set 1. = 200m. 

The difference in the velocity dependence of 
the energy loss according to (1) and (5) which 
hold, respectively, for high and low velocities of 
the mesotron should be noticed. For high veloci- 
ties the energy loss per unit time is inversely 
proportional to the velocity of the mesotron: 
for low velocities it is directly proportional to the 
square of the velocity. There is, therefore, a 
maximum in the energy loss which is found near 
the boundary of validity of the two formulae; 
namely, for mesotron velocities of the order of 
the electron velocity vp. One might wonder why 
the energy loss (5) is independent of the density 
of the degenerate electron gas through which 
the mesotron moves. Actually the collisions occur 
between the mesotrons and the fastest electrons, 
and the collision cross section decreases as 1/yp°, 
This strong dependence on vo just suffices to 
cancel the effect of great electron density, great 
energy loss per collision, and great relative 
velocity of the colliding particles. Naturally, if 
(5) were taken strictly, one would obtain the 
absurd result that the energy loss remains un- 
changed even when the density of the electrons 
becomes extremely small. 

Actually there are two reasons that limit the 
validity of (5) for low electron density. One is 
expressed by the condition 


V<Kv (6) 


that the mesotron should move slowly with re- 
spect to the electrons. The second is due to the 
fact that when a negative mesotron moves 
through an electron gas, the density of electrons 
near it is reduced by the electrostatic repulsion. 
This rarefaction of the electrons near the meso- 
tron effectively neutralizes its charge at distances 
of the order: 


(ah/mvo)*, (7) 
where a is the Bohr radius 
h?/me?. (8) 
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The quantitative treatment of the electron- 
mesotron collision can be carried out according 
to the Born approximation method. This is 
justified if the formula 


e?/hvo <1 (9) 


holds. Now if the Born approximation is appli- 
cable, the region in which the relevant collisions 
take place has dimensions equal to the deBroglie 
wave-length of the scattered particle. In order 
that expression (5) be applicable we must de- 
mand that this wave-length be less than the 
length given in (7). It follows that 


h/mvo<(ah/mv)* or mvo/h>1/a. (10) 


The left-hand side of the last inequality is 
approximately equal to the cube root of the 
density of the electron gas. Condition (10) 
means, therefore, that the density of the electron 
gas must be such that, on the average, more 
than one electron is found in a cube having a 
side equal to the Bohr radius. Inequalities (9) 
and (10) are identical in content. It is easy to 
show that if they are not fulfilled (5) does not 
apply, and the mesotron loses energy at a much 
slower rate than the one given by (5). Conditions 
(6) and (10) are independent, and the more 
restrictive of the two will apply. Condition (10) 
is usually fulfilled approximately in condensed 
matter of not too small density. In the case of 
gases, however, it will be satisfied only within 
the atoms, and therefore the energy loss will be 
confined to these regions. 

Quantitatively, the energy loss of a slow 
mesotron in an ideal degenerate electron gas is 
expressed by the following integral: 


me’ V2 ® © 2r 
a | J J 
(cosé —cos6’) sinédé sin@’dé’dy 
x , 


sin*(2y) 


The meaning of the integration variables is 
the following: @ and @’ are the angles between 
the directions of initial and final electron velocity 
and the direction of the mesotron velocity ; y is 
the angle of deflection for the electron; and ¢ is 
given by 





cosy = cos@ cos6’+sin@ sin@’ cose. (12) 


The integral diverges logarithmically near 
¥=0. The divergence can be removed by taking 
into account the fact that collisions involving a 
small value of ¥ occur at large distances where 
the mesotron charge is screened. 

Collisions at a distance greater than (7) will 
not contribute, and one needs to consider only 
momentum changes greater than 


(Amvo/a)}, (13) 


which is h divided by the length (7). From this 
one finds that the integration need be carried 
down only to a value 


Vmin = (Amvo/a)*(mvo)— 
= (e?/hvo)'=(c/137v9)*. (14) 


One should notice that if (10) is fulfilled, Yuin is 
small compared with unity, as it should be. 
Eliminating the divergence with this prescrip- 
tion, one can evaluate (11) and find finally the 
following estimate for the rate of energy loss: 


qw 4 meV? 1 








ancunasiy aaie doaeamniaMitiiaeis 
dt 3r h* Wmin 
2 m*e*V? 13706 
=— In » (15) 
3n Ohh c 
or finally, 
dW 4 me*T 137 
—-— =— In——. (16) 


dt 3x ph’ c 


The logarithm will be of a small numerical order 
of magnitude and will have somewhat larger 
values in the deep portions of the atom. For 
estimates of the order of magnitude it will be 
permissible to use formula (5). 

The previous theory will be applicable to 
those cases in which the electrons in the vicinity 
of the mesotron can properly be described as a 
degenerate gas. In particular, the theory will 
break down within the K shells of the atoms, 
since the electron density there is less than the 
value corresponding to the degenerate gas. We 
shall return to this question in the next section. 
We shall have to give special consideration to 
the case of insulators and gases where electrons 
may be excited only by discontinuous amounts. 
The case that most nearly approaches the ideal 
conditions is that of the metals, which we shall 
now discuss. 








gE. 





4. METALS 


In order to calculate the energy loss as a 
function of the mesotron energy, we have to 
calculate the average value of the kinetic energy 
T which is to be substituted in formula (5). In 
calculating average values of JT, one makes use 
of the fact that the probability of finding the 
mesotron in a given volume element is weighted 
by the square root of the kinetic energy that 
the mesotron has at this position. The reason 
for this weighting is that the volume in the 
momentum space available to the mesotron is 
proportional to 





T'dT =Tidw. (17) 
The average kinetic energy is therefore : 
f (W— U)'dr 
T= ’ (18) 
f (W—U)*dr 


where U is the potential energy and dr is the 
volume element. 

For high mesotron energy U will be negligible 
compared with W, and T will be equal to W. 
As the energy approaches zero, namely, that 
value for which the mesotron can no longer 
freely pass from one atom to another, the average 
kinetic energy becomes appreciably larger than 
W because U is negative. For W negative, the 
mesotron is bound to a definite atom. For nega- 
tive values of energy the kinetic energy is of the 
order of magnitude of the absolute value of W, 
as will be discussed later. 

For positive W values one obtains a low limit 
for T and the energy loss replacing (W—U)! by 
the smaller expression W!+(— U)! in the numer- 
ator of the integral, and replacing (W— U)! by 
the larger expression W!+(—U)! in the de- 
nominator. The error caused by these substitu- 
tions is a maximum when W is equal to the 
absolute value of U, and is then a factor 2. We 
write thus 


[wart {(- U)'dr 
T ~ 


f Widr+ f oo oi 


(18a) 
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We shall use the potential obtained from the 
statistical model 
Z3/3e@2 (x) 


Ua — sain 
b x | (19) 





where we have set for the distance from the 
nucleus r=xbZ~— and the length d is 


b= (9x?/128)*X Bohr radius 
=0.47X10-§ cm. (20) 


The function ¢ has been tabulated.* In some of 
the following calculations we use the crude 
approximation** 


g=0.4/x. (21) 


Approximating a lattice cell by a sphere one 
obtains 7 by integrating (18a) over the cell, 
The second integral in the numerator can be 
performed by partial integration using the differ- 
ential equation’ for g. One finds 


je N-W1+4ne*b!Z[1 — 9(x0) +x0¢' (x0) ] 





T 
z0 
N-"W'!+42eb'/2Z-4 f gixidx 
0 


Here N is the number of atoms per cubic centi- 
meter, and xo is the value of x at the edge of 
the cell. Using (21) we get 


N-'W!1+4nZe*bi(1 —0.8/x0) 
T _ ’ (23) 
N-'W'!+3.96eb5/2Z—1x,? 





where xo is given by the relation: 
1/N =4rb'x°/3Z. (24) 


Values of 7 as obtained from (23) are given in 
Table I for graphite and iron. One finds that 7 
has a flat minimum at 7 ev and 20 ev, respec- 
tively. This will be, therefore, the value of the 
energy at which energy is lost at the slowest 
rate. At higher energies Table I shows that 7 
becomes less than W. This is because of the 
approximation which we made in substituting 
(18) by (18a). Actually 7 > W always holds, and 
we underestimate 7 if we use 7 =W for high 
values of W. From Table I and formula (5), 
one can calculate, for the two cases in question, 


2 E. Fermi, Zeits. f. Physik 48, 73 (1928). ; 
**From x=0.5 up to almost x=8 the quantity x¢ 
remains between the limits 0.3 and 0.5. 
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‘the time needed for the mesotron to lose energy 
from 2000 ev to zero. One finds 2.6 X 10-" sec. in 
graphite and 2.2X10-* sec. in iron. Somewhat 

longer times would be found in condensed matter 
of lower density. As a practical average time for 
crossing the interval from 2000 ev to zero we 
take for all types of condensed matter about 
3x10-" sec. 

We proceed now to the question of energy loss 
when W is negative, when the mesotron can be 
considered bound to a special atom. If the 
energy is negative and its absolute value is 
sufficiently large, it is a sufficient approximation 
to set the kinetic energy 


T=a\|WI, (25) 


where a is a number of the order of unity. 
(Actually a=1 for a Coulomb field and a>1 for 
the statistical potential.) Expression (25) will 
lead to a very small loss of energy for small 
absolute values of W. Actually at W=0 the 
kinetic energy does not vanish, and may be 
obtained from (23): 
2622 4/3 
T (0) =————_(1 —0.8/x»). (26) 
bxo? 

Since at W=0, T increases with decreasing W we 
may use (26) as a lower limit of T at negative 
energies. We shall use for the energy-loss ex- 
pression (5), and substitute for 7 expression (25) 
or (26), whichever is the greater. Expression (26) 
will be relevant from W=0 to —W=50 ev. The 
time required to cross this energy region is of 
the order of tp =4.84X10-" sec. In the range 
where T is estimated by (25), W as a function of 
time is given by 


—W=T(0)er!, (27) 


This formula is valid as long as the statistical 
model is permissible, i.e., as long as the mesotron 
moves outside of the radius of the K-shell. At 
distances smaller than this radius the actual 
electron density is less than the density obtained 
from the statistical model. Nevertheless the 
energy loss of the mesotron continues to proceed 
according to the formulas (16), (5), and (27) 
even when the mesotron is somewhat closer to 
the nucleus than the radius of the K-shell. The 
reason for this is that the energy loss of the 
mesotron does not depend on the electron den- 
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TaBLeE I. Average kinetic energy as a function of the 
total energy. 











T(ev) 
Wiev) Graphite Iron 
0 36 86 
5 23 62 
10 24 56 
20 27 55 
30 32 57 
50 45 65 
100 83 95 








sity, and thus the failure of the statistical model 
to predict the correct density within the K-shell 
has no direct effect on the behavior of the 
mesotron. 

When the mesotron moves inside the K-shell 
it is not permissible, however, to apply the 
method of deriving formulas (16) and (5) which 
we have given above. In a collision the energy 
mvoV is exchanged, and a time not shorter than 
h/mvoV is needed to describe such a collision. 
During this time the mesotron moves through a 
distance h/mvo. Our discussion so far has assumed 
that during a collision the mesotron moves in a 
straight line. When the mesotron moves along 
the K-orbit of the atom its path has a radius of 
curvature equal to h/mvo. For orbits on or 
within the K-shell it will be a better approxi- 
mation to consider the time-dependent dipole 
consisting of the mesotron and a unit positive 
charge located on the nucleus. One may then 
calculate, by applying perturbation theory to 
the statistical model, the energy exchange be- 
tween the mesotron and the electrons. 

If one assumes that the mesotron moves on a 
circular orbit one may evaluate the rate of energy 
loss dW/dt. The result of the calculation is 
identical with (16) provided that the minimum 
momentum change of the electrons (Amv /a)! 
obtained from (7) is greater than the minimum 
momentum change hw/vp compatible with the 
requirement that a mesotron with angular ve- 
locity w exchanges energy with the electrons in 
quanta fw. At radii somewhat smaller than 
the radius of the K-shell the inequality hw/vo 
> (Amvo/a)* holds for not too light atoms. Then 
(16) must be replaced by 





dW 8e'm?T mv? 
-—= ( ). (28) 


hw 


dt 3aph*® 


















Expression (5) remains a satisfactory approxima- 
tion as long as the logarithmic factor does not 
become small compared to unity, i.e., as long as 
the energy fw transferred to the electrons is 
smaller than the maximum kinetic energy mvo?/2 
of the electrons. For mvo?<hw the perturbation 
calculation gives 


(29) 





dW 25!2axm'noetr? 7 Ze?\? 
dt 3hi =) 
where mp is the electron density near the nucleus 
and r is the radius of the circular orbit of the 
mesotron. Expression (29) agrees with the results 
of calculations on internal conversion, for the 
non-relativistic case, provided that the energy of 
the ejeced electron is great compared with the 
energy of the K-electrons.‘ It is of interest to note 
that (28) holds if the radius r of the mesotron 
orbit is greater than 


ry =(m/p)'*rx =0.171rk. (30) 


Here rx is the radius of the K-shell. For r<r; 
expression (29) is valid. 

For very small values of r, energy loss by 
radiation is faster than energy loss due to 
interaction with electrons. For circular mesotron 
orbits one finds by comparing (29) with the 
radiated energy that radiation becomes predomi- 
nant for r<r2, where 


= 4-U1(m /p)*( hc /e?)3!7 Z-5/2y 
= (4.5/Z5/2) x10" cm, (31) 


where r-=h/mc is the Compton wave-length 
divided by 27. 
From (29) one can calculate the time needed to 


TABLE II. Summary of slowing-down times in graphite 
and iron. 











Time 
Energy range (ev) Graphite Iron 
0— — Ze?/r; 5.3to 7. Ato 
—Ze*/1r; +> —Ze*/re 4.6to 0.3to 
— Ze? /re—>—Ze*/r9 3.6to 0.3to 
2000—> — Ze?/ro 18.9% 12.5to 
(total ‘‘slow’’ range) 9.2X10- sec. 6.1107 sec. 








"4 See, e.g., Rasetti, Elements of Nuclear Physics (Prentice- 
Hall, New York, 1936), pp. 134-139. 
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While the classical expression for the radiation by 
an accelerated charge gives the time spent 
between rz and the lowest quantum orbit 
ro=h?/Zye’*, 








2¢3 
tre rg = (723 —_ ro*) 
4e’ 
Mu re\3 ro\3 
sal (-) -(5) 
4mZL \r, re 
78 16 
-| ag (33) 
Zr Z} 


For heavy nuclei 7; becomes smaller than rq, 
and there is no region in which (29) determines 
the energy loss. Further consequences of high 
Z-values are that ro becomes less than the nuclear 
radius and that pair production begins to play a 
role in the last steps of the slowing-down process, 
All these effects shorten the time the mesotron 
needs to get close to the nucleus. 

For carbon and iron, Table II summarizes the 
times (measured in units tp = 4.84 X 10—"* sec.) the 
mesotron needs to cross the energy regions indi- 
cated in the first column of the table. 

One sees that the slowing down time is less than 
10-* sec. This is very short compared with the 
lifetime of the mesotron, 2 X10-® sec. 


5. INSULATORS 


The case of insulators differs from that of 
metals because the amount of energy that may be 
delivered to electrons in a metal can be arbitrarily 
small, whereas in an insulator it must be at least 
as large as the gap between two Brillouin zones. 
This usually amounts to several volts. The loss of 
energy to electrons will be thereby reduced in 
those cases in which energy is transferred in small 
individual amounts. 

In a collision the amount of energy transferred 
to the electron is of the order of magnitude mvoV, 
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and we can expect that Eq. (5) will hold only if 
mvoV >G, (34) 


where G is the “Brillouin gap,” i.e., the minimum 
energy that the electrons can accept. In case (34) 
is not fulfilled, the rate of energy loss will be 
smaller. This limitation will necessitate a change 
in formula (18) where the integral in the numer- 
ator will no longer be extended to all portions of 
space accessible to the mesotron but only to those 
for which in addition (34) is fulfilled. For positive 
W, the effect on dW/dt is most important at and 
near W=0. For Z 2 6 actual calculation based on 
the statistical atom model shows that the rate of 
energy loss is not changed in this critical region 
by more than a factor 3. The corresponding in- 
crease in the slowing-down time is less than 
10- sec. 

For positive energies, regions close to nuclei 
where (34) is fulfilled will occasionally be visited 
by the mesotron. For negative energies regions 
close to the nucleus may be avoided if the 
mesotron is captured into an almost circular 
orbit. Therefore the possibility arises of the 
mesotron spending a long time in such a circular 
orbit. We shall show that actually the mesotron 
will spend in a circular orbit a time which is not 
long compared to the other times we have ob- 
tained in the slowing-down process. 

Circular orbits of negative energy exist only at 
a considerable depth within the atom. This can 
be proved using the statistical potential (19). 
Circular orbits of negative energy exist only 
where x¢(x) is an increasing function of x.*** The 
value of x g(x) is zero at x =0, increases to a maxi- 
mum reached at x=2.25, and decreases beyond 
this point. Consequently, circular orbits of nega- 
tive energy are to be found only within a distance 
from the nucleus corresponding to x=2.25; 
namely, 

r=2.25b/Z}. (35) 


For a circular orbit at exactly this radius the 
energy is equal to zero, and the kinetic energy of 
the mesotron is equal to 


0.0924? /b = 2.7248 ev. (36) 


*** This statement does not refer to circular orbits out- 
side the core of an atom. In condensed matter there is 
usually no room available for orbits of such very large 
radius. 


If a mesotron is captured into an orbit which is 
not circular, its closest approach to the nucleus 
will be even smaller than the distance given by 
formula (35), and the maximum kinetic energy 
will be larger than that given by (36). Applying 
condition (34), we find that energy loss to the 
electrons will be possible from at least one part of 
the orbit if the condition 


0.09e°Z*3/b ev > }(u/m)G (37) 


is fulfilled. If we set for G the fairly large value of 
7 volts, we see that relation (37) is fulfilled if Z is 
9 or greater. Actually relations (34) and (37) are 
not to be taken in a quite strict sense because 
head-on collisions between electrons and meso- 
trons will give an energy exchange twice as large 
as assumed in (34). If this is taken into account 
an additional factor of } appears in the right-hand 
side of (34). If we again assume 7 volts for G the 
limiting value of Z drops from 9 to 6. 

Stable circular orbits of positive energy exist 
for x values greater than 2.25. The condition of 
stability for a circular orbit is: 


g—x¢’—x9"'>0. (38) 
This condition is fulfilled up to x =3.3, or 
r=3.3bZ-3. (39) 


A mesotron moving on a circular orbit of this 
radius has a smaller velocity than one moving on 
the radius given by (35), and one may expect in- 
creased difficulties in the energy exchange be- 
tween mesotron and electron. Actually, in light 
atoms the greatest stable circular radius (39) 
differs from (35) by less than the uncertainty due 
to the spread of the wave function of the mesotron. 
At the same time the difference in angular 
momenta between the greatest circular orbit and 
(35) is small. For Z = 6 this difference is less than 
h. Thus a mesotron moving in a stable orbit of 
positive energy will not lose energy at a much 
smaller rate than will a mesotron whose wave 
function has its maximum at the radius given 
by (35). 

We can conclude that the circular orbits will 
hardly ever be too stable. Even in case they 
should be formed around an element like lithium 
or beryllium the total time of energy loss will 
probably be less than the lifetime of a mesotron. 

Furthermore, the actual size of the Brillouin 
gap is affected by the localization of the mesotron 
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on one lattice atom. Since, at least in the critical 
cases, the mesotron is captured fairly far inside 
the atom, the atom is effectively turned into an 
element with atomic number (Z—1). If we were 
dealing with an isolated atom, this would lower 
the ionization energy of the atom and turn it into 
the much smaller value which usually is called the 
electron affinity of the atom of charge (Z—1). 
Actually this electron affinity may even be zero. 
In the special case of mesotron capture by the 
hydrogen atom, it is found that when the 
mesotron approaches the nucleus to a distance of 
0.639 Bohr radii, the binding energy of the 
electron bcomes zero. In the closed shell struc- 
tures usually found in insulators, an electron 
affinity of two or three volts is likely to remain. 
As a consequence of this the mesotron, after it is 
captured, may lose as a first step an energy 
smaller than the Brillouin gap. After this loss the 
atom in which the mesotron is now localized does 
not have a closed shell, and further excitation 
may still require less energy than the width of the 
Brillouin gap. Of course, further ionizations 
would tend to raise the ionization potential, but 
the local electron deficiency will be promptly 
filled by capturing electrons from neighbors. 

A special situation arises when the mesotron is 
captured on a hydrogen atom, as may happen for 
instance in paraffin. In that case the hydrogen 
and the mesotron circulating around it form a 
small neutral system which will move along and 
will readily permeate to any part of the lattice. 
As a result one will expect that the mesotron will 
eventually be caught in the field of a more highly 
charged nucleus. 

After the mesotron has attained negative 
energies of an absolute value greater than 100 ev, 
the further energy loss proceeds in insulators as it 
does in metals. 

In conclusion we see that the total time needed 
for energy loss in insulators is apt to be a little 
longer than in metals, because of the difficulty in 
bridging the Brillouin gap. There will be, how- 
ever, no change in the order of magnitude of the 
total time which the mesotron needs to reach its 
lowest orbit. 


6. GASES 


In gases, as in insulators, electrons cannot 
accept from the mesotron arbitrarily small 
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amounts of energy. The lowest electronic excita. 
tion energy, J, of a gas molecule plays the same 
role as the Brillouin gap does in insulators, 

For positive values of W the ene 
proceeds according to (5) and (18). The upper 
limit of the integrand in the numerator jg de- 
termined by the condition 


MVo V=/f, (40) 


The integral in the denominator is extended over 
the whole space. As a result the average value of 
the energy loss —dW/dt is reduced by the ratio of 
gas density to insulator density. Under conditions 
of normal temperature and density the energy 
loss in a gas is about a thousand times smaller 
than in a solid insulator. The time needed to slow 
down the mesotrons from W = 2000 ev to W=0is 
approximately 3X10" seconds. 

For negative values of W the mesotron js 
localized on a specific molecule. As the energy 
transfer to electrons proceeds, it is likely that 
enough energy will be given to nuclear motion to 
cause dissociation of the molecule on which the 
mesotron is found. Thus we may confine our 
attention to the atom which carries the mesotron 
along. The energy loss of the mesotron causes 
progressive ionization of this atom, and, as a 
result, the minimum excitation energy J is in- 
creased. This may lead to an effective stoppage of 
the energy loss if the mesotron happens to be 
moving on a nearly circular orbit. The energy 
loss may stop when the relation 


hu = I (41) 


is fulfilled. Here w is the frequency of the 
mesotron in a circular orbit. 

Further energy loss of the mesotron may thus 
be delayed until the ion carrying the mesotron 
makes a collision with another molecule. Then a 
substantial part of the ionic charge will be 
neutralized, and the energy transfer from the 
mesotron to the electron cloud will be resumed. 
We assume that the excitation energy J of the 
atom which carries the mesotron may be written 
in the form 


Il=K?r, (42) 


where i is the degree of ionization and K is a 
constant. The total energy E needed to raise the 
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jonization from zero to 7 is 
E=tK#. (43) 
We assume that between two collisions of the 


mesotron-carrying atom the energy of the meso- 
tron W will change by E/2. We have then 


dW \ 
~—=}NE=}NKi', 
dt 


(44) 


where NV is the number of collisions per unit time 
~10! sec. 

The energy loss of the mesotron will be smallest 
if the mesotron continues to move on circular 
orbits. For the radius r of this orbit we write as 


previously 
r=xbZ-. 


Using (19), (21), and (41) we get 


heZ d /e\} 
[= |-—(<)] =2.06Zx-*ev. (46) 
widixiL dx\x 


From (42), (43), and (44) we now obtain 


(45) 





dw 
a V)IZix-8, (47) 
t 
If W is calculated for a circular orbit from the 
statistical model one obtains 
Z 4/8 e2 
W = ————d (x¢) /dx. (48) 
2bx 
Here the approximation (21) is not sufficient. In 
the region where circular orbits of negative 
energy are possible we may use for the expres- 
sion xe 
xg=0.489 —0.025(x — 2.25)? 
+0.015(*—2.25)*. (49) 


For x <0.5, expression (49) is not a good approxi- 
mation. This region, however, contributes little 
to the slowing-down time and (49) suffices, there- 
fore, for our purpose. 

We calculate the total time by integrating 


° dt 2.25 rdw dw 
ft fp 
2.25 OX 0 dx dt 
L288 Zs%0%(26)-*[ (0.34/24) —0.045] 
x 
0.493Z'(eV)IK-x-3 


Nv 
Kv\izZz- 
=155(—) ’ 
eV. N 








(50) 
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A reasonable value of K is 5 ev. This gives 


z-16 
t=35 ; 


N 





(51) 


The total time needed for the energy-loss 
process at negative energies will therefore be of 
the order of 10-*® second even if the mesotron 
continues to move on circular orbits. 


7. CHEMICAL COMPOUNDS 


It is of interest in experiments on the disap- 
pearance of negative mesotrons in chemical com- 
pounds to find out with what relative probability 
the mesotron is captured by the different kinds 
of atoms. We are led by crude estimates to the 
conclusion that the capture probability is pro- 
portional to the nuclear charge Z. This may be 
seen as follows. 

It is simplest to set the capture probability 
proportional to the energy loss of the mesotrons 
near the various atomic species. At low positive 
mesotron energy, which is the relevant region for 
our argument, the energy loss is given by an ex- 
pression whose numerator contains the numerator 
of (23). For W=0 this is proportional to Z. The 
denominator is a constant integral for all atomic 
species and does not enter in the evaluation of the 
ratio of energy losses. 

The capture probability will actually be pro- 
portional to the rate of energy loss only if the 
ratio of these rates does not change too rapidly 
near W =0. In particular it is necessary to demand 
that this ratio should not be strongly altered by 
the energy change due to a single passage of the 
mesotron through an atom. For this energy loss 
5W we may write 


w= f [hu(W—U)}ds/t. (52) 
Near W=0 we get 
5W =2.3eVZ8 In(x0/Xmin) (53) 


where Xmin Corresponds to the distance of mini- 
mum approach. The energy given in (53) is small 
enough so that it does not yet affect significantly 
the numerator in (23). Thus the ratios of energy 
losses are hardly affected by 6W, and we may 
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conclude that the ratios of capture probabilities 
are proportional to Z. 


8. CONCLUSION 


The over-all conclusions can be summarized as 
follows. In condensed substances, both con- 
ductors and insulators, a negative mesotron is 
captured in its orbit nearest to the nucleus in 
about 10-" second. In a gas the corresponding 
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time is a little longer than is indicated by the 
ratio of densities. In particular, in normal air it is 
of the order of 10~* sec. In both cases this time is 
very short compared with the mesotron Natural 
lifetime of 2X10~® sec. so that the mode of 
ultimate disappearance of the negative Mesotron 
is governed by the balance between natural decay 
and typically nuclear phenomena leading to 
mesotron disappearance. 
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Some information has been obtained on the spin dependence of scattering of slow neutrons 
by Be, Al, and Bi by measuring the scattering cross section for filtered neutrons. The result 
is that in none of these three cases does the sign of the scattering length change when the spin 
orientation is changed. But in the case of Be and Bi the magnitude of the scattering length 
for one spin orientation may be up to twice as great as that for the other spin orientation, and 
in the case of Al the variation may be by a factor of three. 


OME information on the spin dependence of 
the scattering of slow neutrons can be ob- 
tained by measuring the cross sections of 
some microcrystalline substances for filtered 
neutrons.!? 

When a slow neutron is scattered by an atom 
having nuclear spin J, two values for the scat- 
tering length? can be expected according to 
whether the spin of the neutron is parallel or 
antiparallel to J; these will be indicated by a, and 
a_. If these two values are equal, there is no spin 
dependence of the scattering. In this case inter- 
ference phenomena are not influenced by the 
spin, and the neutron waves scattered by the 
atoms behave as fully coherent. When a, and a_ 
are different, the coherent scattering of the atom 
is determined by an average scattering length: 
(see reference 2, formula (6)) 


I [+1 


= a_+ a... 
2I+1 2I+1 


1H. L. Anderson, E. Fermi, and L. Marshall, Phys. Rev. 
70, 815 (1946). 
2 E. Fermi and L. Marshall, Phys. Rev. 71, 666 (1947). 





a 





The remaining scattering behaves as incoherent 
for interference phenomena. 

In order to discuss the significance of coherent 
and incoherent scattering it is necessary to dis- 
tinguish between collisions in which the spin 
orientation is not changed and those in which it 
changes. The first type of collision is responsible 
for coherent scattering, the second for incoherent. 
The reason is that interference takes place only 
when the scattering is due to the cooperative 
action of all atoms. This is the case when there is 
no spin change to indicate which atom has been 
responsible for the scattering. If there is a spin 
change, however, the scattering is attributed to 
the individual action of that atom whose spin has 
changed. 

One can prove by elementary quantum me- 
chanics that the coherent scattering cross section 


IS 
I I+1 : 
a_ a, 
27+1 2I+1 
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SCATTERING OF SLOW 





NEUTRONS 


TABLE I. Scattering cross sections for slow neutrons. 














—EEE 
<x 10™ 2 x< 10% 2 f pa hite —- 
7 cm o cm? for or i ering 
for unfiltered BeO filtered filtered P aa 
neutrons neutrons neutrons Theoretical section “scattering Limits for value 
Element Spin A~1.8A A>4.4A A>6.7A edge in A < 10% cm? <10™ cm? of a,/a. 
Be 3/2 6.0 0.49 0.64 3.95 0.47 7.5 0.6 to 1.8 
al 3/2 1.7 0.93 1.23 4.67 0.21 1.5 0.25 to 3.3 
Bi 9/2 7.2 1.67 to 1.4 0.93 8.00 0.85 8.9 0.5 to 2.1 
Pb a 9.8 4.6 to 1.9 1.61 5.7 0.85 ~ Ol —_ 








and the incoherent scattering cross section is 


I(I+1) . . 
"a1 (2) 


7 incoherent > 


This last vanishes when a, is equal to a_. 
The scattering cross section ¢ is the sum of (1) 


and (2). 
I+1 I 
o=4r( a4?+ a). (3) 
2I+1 2I+1 








The scattering of neutrons by micro-crystalline 
substance is largely due to the Bragg reflections 
on micro crystals which happen to be oriented in 
the proper way. Such Bragg reflections are 
possible only if the wave-length ) is less than 2d, 
where d is the maximum lattice spacing of the 
crystal. Correspondingly, a large drop in the 
total scattering cross section of micro-crystalline 
materials is usually observed at this critical wave- 
length. The cross section usually does not drop to 
zero, however. The residual scattering is partly 
caused by the effect of irregularities of the 
crystal, both permanent and those due to thermal 
agitation, and partly by the incoherent scat- 
tering of Eq. (2), and to the effect of isotopes if 
they exist. 

If the residual cross section could be measured 
for a substance containing only one isotope and 
made of perfect crystals at absolute zero, a 
measurement of the incoherent cross section of 
Eq. (2) could be obtained. From it one could 
calculate |(a,—a_)|. Practically only an upper 
limit to this quantity can be obtained on account 
of the unavoidable irregularities of the sample. 

Measurements have been performed on the 
residual cross section of micro-crystalline Be, Al, 
Bi, and Pb. A beam of neutrons from the thermal 
column of the Argonne heavy water pile was 
passed through a BeO filter (100 g/cm?) or a 





graphite filter (57 g/cm?). In this beam of filtered 
neutrons the transmissions of the various ele- 
ments were measured for samples of increasing 
thickness. In some cases the cross section de- 
creased appreciably with increasing thickness of 
the sample. This is due to a further filtering of the 
neutron beam by the sample, as will be discussed 
later. 

The results of the measurement are summa- 
rized in Table |. The cross sections observed for 
neutrons filtered through BeO and through 
graphite are given in columns 4 and 5, respect- 
ively. For comparison the cross section observed 
for unfiltered neutrons is given in column 3. 
These last values, in particular those of Bi and 
Pb, are appreciably decreased by the fact that 
the samples had a coarse crystalline structure. 
Consequently there is extinction of the Bragg 
reflection within the individual single crystals 
causing an apparent lowering of the cross section. 

The cross sections of columns 4 and 5 are much 
lower than those of column 3 because the filters 
have removed a large fraction of the neutrons 
capable of being Bragg-scattered. 

For Be, the theoretical edge (2d) is 3.95A. 
Consequently both filters pass only wave-lengths 
longer than 2d and which cannot be Bragg-scat- 
tered. The increase in cross section from column 4 
to column 5 is probably due to the Doppler effect 
resulting from thermal agitation of the Be nuclei. 
The residual cross section 0.47 given in column 7 
has been chosen as the lower of these two values 
corrected for the small absorption of Be. 

In the case of Al, the BeO filter does not com- 
pletely eliminate all neutrons that may have 
Bragg reflections. Most of the rise in cross section 
from column 4 to column 5 is due to increase in 
the 1/v absorption Actually the minimum scat- 
tering cross section is observed with the graphite 
filter. This value of 1.23 is decreased by 1.02 





because of absorption of the Al. The residual 
scattering cross section is therefore taken as 0.21. 

For Bi, the theoretical edge 2d is 8.0A. There- 
fore neither filter is able to remove completely 
neutrons that can be Bragg reflected. For BeO 
filtered neutrons there is a decrease in cross 
section with increasing thickness of Bi pre- 
sumably due to a partial filtering of the neutrons 
by the sample itself. The graphite filter, because 
of its larger critical wave-length (2d), gives a 
lower cross section. This last has been decreased 
by 0.08 because of absorption, and therefore 0.85 
has been taken as the upper limit to the residual 
scattering. 

The theoretical edge for Pb is at 5.7A. The 
BeO filter is therefore inadequate as is apparent 
from the fact that the cross section varies from 
4.6 to 1.9 with increasing thickness of the ab- 
sorber. As residual scattering cross section we 
have taken the one measured with the graphite 
filter corrected by 0.76 for absorption. 

The values of the residual scattering cross 
section so obtained are collected in column 7 of 
the table. In column 8 are listed the scattering 
cross sections of the bound atoms obtained by 
correcting the cross sections of the free atoms by 
the factor (A+1/A)? for the reduced mass. 

For the elements Be, Al, and Bi, which have 
only one isotope, the residual scattering cross 
section is due only to spin dependence of the 
scattering length and to crystalline irregularities. 
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The assumption that the residual cross section is 
due entirely to spin dependence leads to an 
overestimate of the difference between a, and a 
The limits for the ratio a,/a_ given in column 9 
are calculated with this assumption and using the 
data of columns 7 and 8 together with formulae 
(2) and (3). 

In the three cases investigated it is found that 
a, and a_ have the same sign. This is in accord. 
ance with the fact? that almost all elements have 
a positive scattering length. If one is allowed to 
generalize from these few cases, it would appear 
that the same is true for the scattering lengths for 
individual spin orientations. This is quite plav- 
sible also on theoretical grounds. On the other 
hand, differences of as much as a factor of 2 in the 
two scattering lengths produce a relatively small 
amount of incoherent scattering. The observation 
of the residual scattering becomes a sensitive 
method for investigation of spin dependence only 
when the difference between the two scattering 
lengths is quite large. It is quite possible that this 
difference is much less than the maximum values 
given, and that the observed residual scattering 
is due largely to crystal irregularities. An analysis 
of the data for Pb has not been given because 
this case is complicated by the presence of several 
isotopes. 

This document is based on work performed 
under Manhattan Project sponsorship at the 
Argonne National Laboratory. 
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The °S State of Carbon 


A. G. SHENSTONE 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 


(Received May 19, 1947) 


The °S state of carbon is important in thermo-chemical theory, and its position has been 
the subject of much speculation. Evidence is here given that it is 33735.2 cm™' above the 
ground state *P and combines with it to emit two lines at \2965 and 2967. 





N the theory of the chemical binding energies 
of carbon compounds, the sp* §S state of the 
carbon atom plays a large part. Because the 
tetravalent condition of carbon is ascribed to that 
state, the excitation energy of S enters intimately 
into the calculations of the carbon-bond energies 
derived from observed reaction energies and from 
band spectra and, perhaps more important, it has 
been thought to be involved in the derivation of 
the heat of sublimation of carbon. This whole 
field has recently been discussed by Long and 
Norrish! in a paper on The Thermochemistry of 
Carbon, from which it is clear that there is serious 
disagreement on the value of 5S between various 
workers in the field. 

It is recognized that the well-developed theory 
of iso-electronic sequences has given quite accu- 
rate predictions of the positions of such terms as 
this °S. Using the means so provided, Edlen? 
predicted that °S should be about 32900 cm 
above the ground state s*p? *P of carbon. Bacher 
and Goudsmit’ calculated 35,000 from a wave- 
mechanical theory, and Ufford,‘ also from the 
wave-mechanics, predicted 22,000 to 26,000. The 
value favored by Long and Norrish from the 
chemical data is in the range given by Ufford’s 
theory and so falls considerably below the theo- 








retical value of Bacher and Goudsmit and the 
extrapolation of Edlen. There is no doubt that 
Edlen’s value would appeal to most spectro- 
scopists because no radical exceptions have ever 
been found to the rules of iso-electronic sequences. 
It appears to have been assumed by spectro- 
scopists, however, that the observations of the 
spectrum of the carbon arc made by Paschen and 
Kruger® were so complete that there was no hope 
of finding the absolute value of °S by actually 
observing its'spectroscopic combinations. 

Before the war when I was working on spectra 
which extend into the fluorite region, I found that 
there is in the literature a great lack of data on arc 
spectra in that region. This is due to the fact that 
most of the observations have been made with 
vacuum sources, which are notoriously poor 
emitters of arc lines. For that reason I had 
adopted the type of arc source first used by 
Selwyn,* and I had photographed with it the 
arcs of Cu, Ag, Au, Fe, Ni, Pt, and several other 
metals in nitrogen for myself and for Dr. Meggers 
at the Bureau of Standards. In those spectra the 
nitrogen lines always appeared, and, when carbon 
electrodes were used, the carbon arc and first 


TABLE II. Evidence for correctness of assignment in 


























Table I. 
TABLE I. 
Regular doublet law (» —~RN%5 /36) X107* 
MA) Int. »(cm~!) Ay Assignment Ay (av)? Diff. Diff. 
1431.597 20 69852.1 2s2p? *S.—2s2p73s88P; CTI 45.5 2.695 54.61 
24.8 .759 43.36 
1432.105 15 69827.3 5S. — ‘iP, NII 126.6 ssel 97.97 
20.7 .757 43.40 
1432.530 10 69806.6 5S.— 5p, O III 285.6 4.111 { 141.37 
—" — 752 > 43.55 
FIV 559.4 4.863? 184.92/ 
'L.H. Long and R. G. W. Norrish, Proc. Roy. Soc. 187A, = == 
337 (1946). Se 
*B. Edlen, Nova Acta Regiae Soc. Sc. Upsala IV, 9, No. ‘F. Paschen and G. Kruger, Ann. d. Physik 5, 7, 1 
6 (1934). (1930). 


*R.F. Bacher and S. Goudsmit, Phys. Rev. 46, 948 (1934). 
‘C. W. Ufford, Phys. Rev. 53, 569 (1938). 


6 E. W. H. Selwyn, Proc. Phys. Soc. 41, 392 (1929); A. G. 
Shenstone, Trans. Roy. Soc. A237, 453 (1938). 
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TABLE III. Wave-lengths of multiplet at \1656. 
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TaBe IV. Wave-lengths and assignment of lines 











— 
(A) Int. »v(cm~) av XA) Int. »(cm™~!) Assignment 

$$ 
1656.255 5 60377.18 2964.846 2 33718.74 2 adi 
1656.918 3 60353.02 +27.08 Av=27.02 PP 2p, 
1656.998 10 60350.10 2967.224 5 33691.72 2p'P.—2y/ 
1657.368 2 60336.63 2582.901 5 38704.58 2p!S,— gap" 
1657.891 4 60317.60 >27.11 a" 
1658.113 5 60309.52 








spark lines were observed. In the latter case there 
was always in addition a group of three lines at 
1432 which definitely associated themselves in 
character with the carbon lines and yet did not 
appear in any published wave-length lists. Their 
separations and position were incompatible with 
any unknown but predictable transition in C J or 
C II, except that from 2s2p?3s5P down to 2s2p? 5S. 
Because | had many other spectra in hand and | 
knew of Professor B. Edlen’s interest in the light 
atoms, I sent him these data in 1939 to use as he 
wished. He recently published the material given 
in Table I in Nature’ and gave the following 
evidence for the correctness of my assignment of 
origin in Table II. The data in Tables I and II 
differ slightly from those originally given to 
Edlen because more recent measures have been 
included. 

From the comparison with other spectra in the 
same iso-electronic sequence and from general 
familiarity with spectral relations, he estimated 
the position of the sp* 5S above the ground state 
at about 32,700+2 percent cm~. This prediction 
was questioned by Long and Norrish on the 





‘. 2. Hg. 








grounds that the evidence was insufficient and 
that the value was incompatible with the known 
thermochemical energy relations of carbon, jt 
therefore seemed important to fix the position of 
the ®S term definitely, by finding, if pogg. 
ble, the lines corresponding to the transition 
s*p? *P —sp* 5S. There should be two such lines 
with a separation corresponding to the separation 
s*p?(®P;—*P2), those being the components of #p 
with which 5S2 can combine without violation of 
the AJ rule. This separation is given in the litera. 
ture as 27.5 cm™ from measurements of yp. 
resolved multiplets in the Schumann region, 
I have recently obtained photographs of the 
s*p? §P — s*ps*P multiplet at 41656 completely re. 
solved, and I have made many measurements of 
the differences. Because of the scarcity of accu- 
rate standards, the wave-lengths may be in error 
by as much as +0.02A but the differences should 
be much better, perhaps only +0.002A. From 
these measurements I find the intervals to be: 


27.1 and 16.4+0.1 cm, 
40.6 and 19.0+0.1 cm—. 


sp? 3P. 10 
sp*s®Po, 1,0 


The wave-lengths are given in Table III. The 








Fic. 1. Carbon arc in argon. 
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7B. Edlen, Nature 159, 129 (January 1947). 


I 2479 2p'Sp—3s'Pi. 


2964-7 2p*P;, »—2p! §So. 
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transition from the structure 2s2p* to 25?2p? is of 

course an allowed one, being from odd to even, 

but the intensity of the inter-combination 1 ee &S 

lines is zero in strict Russell-Saunders coupling. 

However, with a triplet width as large as 43.5 

there is a sufficient departure from Russell- 

Saunders coupling to make one expect lines of 

intensity about 10~* or 10~* of the allowed lines. 

Believing that the past failure to observe the 
desired transition was probably due to the very 
dense band spectrum which appears in the ordi- 
nary graphite arc, I tried the arc in nitrogen and 
in argon. Neither of these eliminates the band 
spectrum at all completely, but they both have 
the effect of lowering the temperature of the 
electrodes. The most obvious difference in the 
spectra is the almost complete disappearance of 
the ordinary impurities such as Ca and Mg when 
N- or A is used. This may be associated with the 
fact that the graphite burns away rapidly in air, 
more slowly in nitrogen, and scarcely at all in 
argon. The experimental arrangement was a 
simple enclosure of the arc in a glass bulb with an 
open neck toward the spectrograph. The gas was 
purified of oxygen by passage over hot copper 
turnings and was allowed to escape from the arc 
enclosure into the open air. The argon un- 
doubtedly contained a percent or two of nitrogen. 
The arc was operated at about four amperes 
without any ballast resistance from a constant 
current source described elsewhere and was 
focussed on the spectrograph slit by a quartz- 
water achromatic lens. The potential across the 
electrodes was about 40-45 V in nitrogen, 55—-60V 
in argon, and 60—62V in air, when the gap was 
about 3mm. The rods were spectroscopic graphite 
} in. in diameter. 

Figure 1 isa 6X enlargement of two sections of 
an exposure taken on a medium Hilger quartz 
spectrograph using an Eastman spectrum-analysis 
No. 1 plate. The shorter wave-length section 
exhibits the two known C J lines at \2479 and 
\2583, which stretch from one electrode to the 
other, and the C JI lines at \2509-11 which are 
very strong at the negative pole. The other line 
which stretches across the arc appears only 
occasionally. It is mercury (2536. In the longer 
wave-length section of Fig. 1 the two marked 
lines have exactly the characteristics of the other 
C I lines and are the only ones observable on the 
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Fic. 2. Plot of (J)' and of 2p°P—295S for iso-electronic 
sequence C I to F IV. 


whole plate which do. They can be photographed 
faintly in less than a second and then appear 
practically alone in that region of the spectrum. 
These two lines and \2583 have been measured on 
several plates against copper and iron standards. 
The best values, obtained from a thirty-five 
minute exposure on a 21-ft. 30000 line grating, 
are given in Table IV. The lines at a dispersion of 
1.25A per mm appear only slightly wider than the 


TABLE V. Revision of term values. 











Electron Name 
configuration (Edlen's notation) Value 
282 2pP» 90880.5 
16.4 
2p*P, 90864.1 
27.1 
2p°P. 90837.0 
2992p" 2p'So 69232.1 
2s2p° 2p’ 8S, 57145.3 
2s*2p3s 383P 5 30546.5 
19.0 
388P; 30527.5 
40.6 
3s°P 2 30486.9 
2s2p%3s 3s 5P, — 12661.3 
20.7 
3s’ 5P, — 12682.0 
24.8 
3s’ ®Ps — 12706.8 
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Fic. 3. Carbon arc in (1) argon, (2) nitrogen, and (3) air. 


iron comparison lines. The consistency of the 
measurements indicates an accuracy of +0.002A. 

My assignment of the lines at \2964 and 2967 
as 2p*P —2p’ 5S is substantiated by the following 
positive evidence and a lack of any contrary 
evidence. 


1. The two lines do not correspond to known lines of any 
element, and they are not band lines. 

2. They are very close to the position predicted by 
Edlen. 

3. They have the same character as other C J lines and 
are the only newly observed lines with such a character. 

4. They have a wave-number separation which agrees 
within the small experimental error with the one found 
from Schumann region multiplets. 

5. They have approximately the relative intensities to 
be expected in such a transition. 

6. They have an absolute intensity of a reasonable 
order of magnitude, roughly 10-* or 10~* of \2479.322 
2p'So— 3s'P; and of the same order as A2583.65 2p'S—3s*P;, 
which is the only other known inter-system line. 

7. The plot of (7) for §S for the iso-electronic sequence is 
very nearly a straight line, asisalso the plot of 2p°P —2p’ 5S, 
which illustrates the irregular doublet law. Both plots are 
given in Fig. 2. 


Figure 3 exhibits the difference between arcs in 
air, nitrogen, and argon, and shows that the 


identification would have been possible, though 
with much more difficulty, in an arc in air, The 
three spectra are of about the same exposure 
time but the intensity of the air arc is strikingly 
greater corresponding to its greater visual bril- 
liance. It should be noted that the polar C J/ 
lines are practically absent in the air arc and that 
the known C J lines are relatively weaker than in 
the other gases. 

The use of the measurements here reported 
allows the following (Table V) revision of old 
term values and the addition of the new terms 
2p’ §S and 3s’ 5P. This places °S at 33735.2 cm™ 
above the ground state. Its excitation requires 
4.16 ev or 96.4 kcal. per mole. 

If one defines the resonance lines as those due 
to the transition to the lowest state from the first 
excited state of opposite parity, then the new 
lines here discussed are the resonance lines of 
carbon. It is odd that they are almost the last 
important lines of C I to be discovered. Their 
existence implies that °S can hardly be con- 
sidered as a metastable state, and should not, 
therefore, enter in that role into thermochemical 
theory. 
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The ,-Instability of Mesons 


R. J. FINKELSTEIN* 
University of California, Berkeley, California 
(Received May 10, 1947) 


The following examples of y-instability are investigated: (a) Mp—yi+~72, (b) Mi-yi+72 
+3, (c) Mi>Mot+y, (d) Mi*—M,*++7, where My and M, mean a pseudoscalar and a 
(heavier) vector meson, respectively, and the superscript + indicates the charge, if there is 
one. If the coupling constant and meson masses appropriate to the Schwinger mixture are used, 
the lifetimes are found to be: rz=1X10~* sec., m=2X10-™ sec., re-= 1107" sec., and rg=4 
10-8 sec. In all cases except (b) the calculation leads to logarithmically divergent integrals. 





1. INTRODUCTION 


CCORDING to current meson theory both 
the nucleon and the electron-neutrino carry 
a virtual meson cloud; conversely, a free meson 
carries a nucleon cloud and can emit an electron- 
neutrino pair. The possibility of emitting an 
electron-neutrino pair permits the B-decay of the 
free meson, and the nucleon cloud, coupling even 
a neutral meson to the electromagnetic field, 
permits y-decay. Unfortunately, the whole status 
of meson theory is at present quite uncertain, and 
our numerical results probably have at most an 
order of magnitude significance, not only because 
of uncertainties about the form and numerical 
constants of the theory, but also because of the 
dangers of working with a divergent formalism at 
its limits of validity. Nevertheless the investiga- 
tion given here may be of interest, since it illus- 
trates a kind of instability common to a class of 
theories. The following examples of y-decay, 
which are probably typical of the more interest- 
ing cases, will be considered. 


(A) Mo-r+72 Ta =1X10-"* sec. 


(B) Mimmityetys te=2X10-" (1) 
(C) M:->Mot+y tTe=1X10-* 
(D) Mi*->Mo*+y tp =4X10-'* 


where My, and M, mean pseudoscalar and vector 
meson, respectively, and the superscript + indi- 
cates a charge, if there is one. Processes of this 
nature have recently been proposed to explain 
the soft component of the Auger showers.! 
Process (B) is the simplest mode of decay of a 


- On leave from Argonne National Laboratory, Chicago, 
inois. 

'J. R. Oppenheimer, New York Meeting of American 
Physical Society, January, 1947. 
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vector meson, since the binary fission is forbidden 
by Furry’s selection rule.? Processes like (C) 
and (D) are logical possibilities in a meson theory 
such as the Schwinger mixture in which two 
kinds of mesons of different rest mass are postu- 
lated. Another possibility, not described by (C) 
or (D), is the decay of a heavy meson into a 
lighter one with a change of charge, e.g., M,+— 
M,+e*. In fact, Hulthén has proposed a mixture 
of charged pseudoscalar and neutral scalar fields 
with different rest masses, and one may ask 
whether the heavy (neutral) particle decomposes 
into the lighter one. However, the decay of one 
meson into a lighter one with the ejection of an 
electron should not occur; instead the heavy 
meson should undergo simple §-decay. On the 
other hand, the processes (A)—(D) are all much 
faster than 6-decay, since they are caused by the 
coupling of the meson field to the nucleons whose 
mesonic charge is of the order of e, whereas 
B-decay arises from coupling of the meson field 
to the electron and neutrino sources whose meson 
charge is only of the order of 10~%e. Whenever 
there is competition between 8- and y-decay, the 
latter will win. For example, consider M)y—et 
+e-; this reaction, whose lifetime is of the order 
of 10-* sec., should not be expected since (A), 
with a lifetime of 10~** sec., is much faster. 

The lifetimes are given at the right of the 
corresponding reaction. According to these re- 
sults a charged vector meson does not disappear 
through direct 6-decay but only as the result of 
(D) followed by the B-decay of Myo*. It would 
also follow from the above scheme that the 
eventual fate of a neutral vector meson is its 


2S. Sakata and Y. Tanikawa, Phys. Rev. 57, 548 (1940). 
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transmutation into three photons, either directly 
by (B) or indirectly by (C) and (A). Since, 
however, (B) and (C) are competing decompo- 
sitions and (C) is much the faster, (B) should 
not be expected. 

2. HAMILTONIAN AND MATRIX ELEMENTS 


The interaction terms in the Hamiltonian 
responsible for these processes are to the first 
order in the coupling constants 


(4x) [err Tal gol UaBys | 


+(fo/Ko)[vs00oVatoD.Ua|}bdx, (2a) 


(4m)! f 8°. rales Use Ue] 


+(f:1/Ki)[Voey+Vahe}}bdx, (2b) 
where 


a and ra=isotopic spin index and matrix, 
(Uo, U) = wave function of M,, 
U and @=wave function of M, and nucleon, 


respectively, 
V=DxU and V,= —0,.U—DU), 
7? —iBa, o=%75a, Y¥5s= — 10110203. 


f and g have dimensions of charge, and K =yc/h, 
D,=V —t(e./c)A, where A is the electromagnetic 
vector potential; D=¥V except in case (D). For 
simplicity we put g=0 and consider only the 
second invariants; in addition, the analysis is 
limited to transverse mesons. The non-relativistic 
terms, oDU and V6se, play no role at all in 
decompositions (A) and (B), since these terms 
are proportional to the velocity of the meson, 
and one can assume that the meson is initially at 
rest. In (C) and (D), however, either the initial 
or final meson must be moving. 

These interaction terms give rise to transitions 
which have matrix elements with the following 
absolute values: 


Mo pd 
(2x/G)*(fo/Ko) Lu n* 7 stm 
+(cp/€)(un*eum) ], (3a) 
M, (transverse) : 
(29 /G) *f hice (€/p10?) (Un*~OmUm) 
+[(p/p1c) X€m](UnBoum) ], (3b) 


FINKELSTEIN 


where (p, «) are the momentum and energy of 
the meson, while its polarization js €n. When 
the meson is charged, there are additional terms 
depending on A; these new terms are not needed 
before Eq. (28). Matrix elements for the emission 
and absorption of a photon by a proton are also 
required : 


(2%/G)ehce“*(u,,* eeu). (3c) 


The functions u represent plane wave solutions 
of the Dirac equation. These matrix elements 
correspond to processes in which a nucleon jp 
one state is created while a nucleon in another 
state is annihilated, and at the same time a meson 
or photon is created or annihilated. If both 
nucleon states have positive (or negative) energy, 
the corresponding process can be described as the 
emission or absorption of a photon or meson by 
a nucleon (or antinucleon). If one energy state 
is positive while the other is negative, the process 
described is the creation or annihilation of a 
nucleon pair with the simultaneous emission or 
absorption of a photon or meson. 

Terms of exactly the same form as (2) couple 
the meson field to the electron-neutrino. The 
8-decay of a free meson is an illustration of pair 
creation associated with the matrix elements 
(3a) or (3b), if wu is interpreted as belonging to 
an electron-neutrino instead of to a proton- 
neutron. In this process a charged meson dis- 
appears, while simultaneously a neutrino in a 
negative energy state is annihilated, and an 
electron in a positive energy state is created. 
Hence the initial charged meson is replaced by 
an electron and an antineutrino, which is equiva- 
lent to a neutrino. One should expect an analo- 
gous process with nucleons; instead of Mt+t—e«+ 
+n, one has M+-—+P++N, or if the meson is 
neutral, M°—P++P-. Since the mass of the 
meson is less than that of the nucleon, these 
transitions can take place only virtually; but 
because the coupling constant is 10* times greater 
for the nucleon, the decomposition into heavy 
particles is correspondingly more frequent than 
the decay into light particles. Although such a 
decomposition cannot lead directly to a final 
state in which the nucleons separate, it can play 
an intermediate role in a higher order process 
which is still very fast, since the virtual decompo- 
sition is so frequent. 
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DECOMPOSITION OF A NEUTRAL PSEUDO- 
¥ SCALAR MESON INTO PHOTONS' 


The lifetime, 7, is given by the usual expression 
r= (2x/el VI" (4 


where p is the density of final states and V is the 
matrix element for the process. Both photons 
get the same energies and opposite momenta so 
that 
p = G4mpo’c/16h', (5) 
where G is the fundamental volume and yp is 
the mass of the meson. 

The decomposition of a neutral meson into 
two photons requires three steps, e.g., Mo—P* 
4+P-, Pt—Pit+71, and P,;++P-—y:, where 
P,+ means that the momentum of P+ has been 
altered by the emission of y:. These three steps 
can take place in any one of six orders; further- 
more, since the theory is symmetrical with 
respect to nucleons and antinucleons, each of 
these six possibilities gives rise to another if the 
roles of the nucleon and antinucleon are inter- 
changed. Hence there are 12 terms in V, the 
probability amplitude. Further, since the con- 
servation laws do not fix the momentum of the 
proton pair, there is an infinity of intermediate 
states, and the probability for the transition 
from the initial to the final state must be 
integrated over the momentum space of an 
intermediate pair particle. 


v=(G/m) f f(oND-)PaPas, (6) 


where the N; are products of matrix elements, 
the D; are energy denominators, + runs over the 
12 possibilities mentioned, and w is solid angle. 
In order that (6) be convergent it is necessary 
that N vanish faster than 1/E. 

It is convenient to expand the ratios N/D in 
powers of x=yo/M: 


N,Dj* = (°N;+'Nix)(°Di—'D x). 
Writing °D;=1/4E*, where E is the energy of a 


pair particle, and anticipating the result that 
°N;=0, one has to the first order in x 


LND- =x(¥ Ni) /4E (7) 


* This instability has also been considered in unpublished 
calculations by Oppenheimer, Christy, and Lewis. 
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OF MESONS 


By Eqs. (4)-—(7) 
t= S124*(h/poc*) (hic*/etfe*?) I~, (8) 


where 


s=f (> wi) (y? — 1) "dy/y), (8a) 


ly = Jnde, (8b) 
in; ="N;/[2e®foh*(2x/Gus)"), (8c) 
y=E/Me, (8d) 


n;, y, and J are dimensionless. 

We now calculate the dimensionless numera- 
tors nm; These may be distinguished by three 
subscripts having the order in which the three 
intermediate steps occur. For example, by (3), 


Mois= Lo (u-|y5|u)(u| we;| u;)(u;| ae;|u-), 


where the sum is over the spins of the inter- 
mediate states. The quantities u and u>~ are 
plane wave functions having momentum P and 
energies +E and —E, respectively. Similarly u,; 
corresponds to momentum P; and energy £,. 
To satisfy the conservation of momentum, one 
must have P;=P—p,, where p; is the momentum 
of the ith photon. Here e; and e; are unit 
polarizatior vectors of the corresponding pho- 
tons. Then 


Noij= Spr~-yer(ae;)A;(ae;), (9) 


where the X’s are the projection operators defined 
by 
A* = +}(+1+e(V/c)+8/y), (9a) 


A=A(Vi, Fi). (9b) 
Evaluation of mo; gives 


noij=tLe; Xe; ]{ Ly ys —y* (V/c) 
—y*(V;—V)/c}. 


Expansion of E; and V; to the first order in x 
gives 

noij= (tx /2y*) [es Xe; Jk, (10) 
where k; is the unit propagation vector of the ith 
photon. Since k; = —ke, it follows from (10) that 


No12 = No21. 


In addition nijo = Sp[A~(ae,)As(ae;)A(ys) ]. Take 
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the complex conjugate : 


nijo* = SplA~ysd(ae;s)Ai(ae;) ] = — Mois, 


Ni jo = Noij- 


Hence 012 = M021 = "120 = M210. All of these numer- 
ators lead to convergent integrals. The situation 
is quite different if the meson is annihilated in 
the second step; the corresponding integral is 
then indeterminate. One rather natural way of 
evaluating this integral leads to a cancellation 
and a finite result again equal to m2. We shall 
tentatively assume this result, but will after- 
wards return to it. The six ”; values obtained by 
permutation of 0, 1, and 2 are then equal, and 
by Furry’s theorem are also equal to six others, 
obtained by interchanging the roles of nucleon 
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Fic. 1. Transition schemes (on left) and momenta 
corresponding to mjoj* and mjo;~. On the left the minus 
superscript means a state of negative energy and on the 
right it refers to a negatron. P; means momentun P— pj, 
while P_.; and P_ means momenta —(P—p;) and —P, 
respectively. 
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and antinucleon. Hence 


Cwi=(12)(4n) ix/2y*)LerXertis. (1p 


After squaring and summing over both polar. 
izations of both photons, one has by (11) and (8) 


r= 4n*(hh/woct) (We /efe) yy 


The Elements nj; 


Let n* denote the numerator when P+ absorb. 
a meson, and let m~ designate the numerator 
when P~ absorbs a meson. The two transition 
schemes are shown in Fig. 1a and 1b. Then 


niojt = SplA~(ae;)Avysdi(ae,) |, (13) 
Ning = —Sp[A~(ae,)Ai(@e;)A~y5]. (13a) 


The negative sign is introduced in (13a) because 
two of the vacuum protons are interchanged jp 
the transition to which this equation refers. 
Evaluation of the spurs gives 


Nioj* = iLe;Xe;]{ _- [2 —yi?t+yiy1](V/c) 
—[1i+9y7 ](Vi-V)/c}, 


nioj =iLe, Xe; ]{[2—y?*+yy,7 ](V/c) 
+[1-—y?}(Vi—V)/c}, 


and 


niojt +niox =tLe; Xe; ]{ (yi? —y) (V/c) 
+(-—y?-y797)(Vi-V)/c}. 


Both n+ and m~ contain the logarithmically 
divergent term (V;—V)/c which cancels in the 
addition. To the first order in x this sum is 2m. 
But according to Furry’s theorem m* and n- are 
of the same absolute value so that they either 
double or cancel and never give the partial 
cancellation just described. The reason for this 
difference is that Furry calculates the m~ by 
reversing all f/nomenta, energies, and directions 
of transition. This procedure is equivalent to 
choosing the momenta as shown in Fig. Ic. 
Since there is an ultimate integration over solid 
angle, the relative directions of the protonic 
momentum in n+ and m- should not affect the 
integral (n+-+-n-), if this integral is well defined. 
But since it is not well defined, its value can be 
changed by reordering the terms in the integrand. 
Hence there is an uncertainty in the correspond- 
ing lifetimes. The terms in which the meson is 
annihilated in the second step are similar to 
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those appearing in the self energy, since one 
hoton is used to create the pair and the other 
to collapse it. As usual, there is no satisfactory 
procedure for dealing with this situation. 


4, DECAY OF A NEUTRAL VECTOR MESON INTO 
: THREE PHOTONS 


According to Furry’s theorem a vector meson 
cannot decompose into two y-quanta because 
the nt and the m~ cancel. Hence we consider 
the decomposition into three photons, suggested 
by Sakata and Tanikawa.? 

The momenta of the three photons close a 
triangle. For given pi the ergodic surfaces are 
ellipsoids with foci at the two ends of p;. The 
density of final states is determined by the 
number of states common to the ellipsoidal shell 
(E, E+dE) and the spherical shell (p2, p2+dp2). 
The volume of this ring is 


x,y 
po, 7 





aR=2ryJ( )épsdn, 


where (x,y) are Cartesian coordinates in the 
plane of the moments, and »=2+ 3. Hence 


a pine Sais 14 
r*ody pido, =2r sch® Pidpodw;. (14) 


The factor 6 appears in the denominator of p 
because the photons are indistinguishable. The 
product pip2ps; favors a symmetrical decomposi- 
tion of the meson, but this initial bias is lost 
when p is combined with V?. The lifetime now 
is given by 


t= (1/96) (urc?/h) (e®f2/cth*) 
4 a 
x J*dfidfe, (15 
J _ Tiidiy (15) 


‘W. H. Furry, Phys. Rev. 51, 125 (1937). The original 
statement of this theorem refers primarily to transitory 
electron-positron pairs interacting with the electromagnetic 
field. The more general statement of Furry’s theorem 
needed here is the following. Let 0 be any Dirac matrix 
such that (kr|0|k’r’) = +(—k’—1'|0| —k—1), where k is 
the momentum and r is the sign of the energy. The theorem 
now states that a process is forbidden if it involves an 
odd number of even matrix elements, where “even” means 
that the positive sign is correct in the equation just written 
for 0. The even matrix elements are contained in the 
vector and antisymmetric tensor interactions; those con- 
tained in the scalar, pseudoscalar, and pseudovector 
combinations are odd. 
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where 

8J= f “(Zede-)(y*—1)hy*dy, (15a) 
1 

_ f niu, (15b) 

D,;=8E*d;, (15¢) 


n; = N;/((ehc)*(fihc) (€r€2€suic®) "(24 /G)* |. (15d) 


Here €;(=¢*f;) is the energy of the ith photon. 
The n; and d; are the dimensionless numerators 
and energy denominators. It is again convenient 
to expand in powers of y,c?/E(=x/y). 


d;=1+'d,(x/y), 


The ; will now be considered in detail. 

The decomposition takes place in four steps of 
which one is needed to annihilate the meson and 
three to create the three photons. These four 
steps may be permuted in 24 ways. Again new 
possibilities may be realized because of symmetry 
between nucleons and antinucleons; it is possible 
to choose three intermediate states of positive 
(negative) energy and one with negative (posi- 
tive) energy or two with positive and two with 
negative energy. Hence there are 24X3=72 
chains in all. A typical numerator is 


Nmijk = Sp[A~(Baem)A(ae;)A:( ae ;)A;;(@ex) J, (16) 


n= n,+'n(x/y). 


where 


A=A(v,), Ps=P—p;, Pij=P—pi—p,. 


Equation (16) may be written 
Nmijk = A = LdsSp(G'a'). 


Aw=0 unless :=4 or 6. A. involves x only 
through the energy of the virtual proton. But 
E/E;=1+Ax/y, where A;=fVk;/c and where 
k; is the unit propagation vector. Hence terms 
linear in x, being also linear in V, vanish in the 
angular integration. A similar argument is ap- 
plicable to the linear portion of A,., insofar as x 
appears in the energy. But A,» also has terms 
linear in x in virtue of 


Vi/c=V/c+fi —ki+ (Vki/c)(V/c) }(x/y). 


Here the factor of x is even in V, and so are all 
other factors in the spur except V itself. Hence 
the term is odd and may be ignored. In addition 
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one has }>°w;=0 (see Eq. (18)). It follows that 
Lads = (x/y)Dwite;, (17) 


where 'c; is the part of 'd; independent of V. It 
is therefore only necessary to determine the zero- 
order part of m; and one can put A;=Aij=A in 
all the spurs. The ; contain terms independent 
of V/c and in addition terms of the type 
Sp[(a@m)(aV/c)(ae;)(aex)(aV/c)(ae;)], which 
are quadratic in V/c and must be integrated over 
solid angle. Let 
@; = 2%,’ + %w,"", 


where w’ corresponds to those processes in which 

one intermediate state is positive (negative) and 

the other three states are negative (positive) ; 

and where w” refers to those processes in which 

there are two positive and two negative states. 
The results are 


Gmijk = {ijk ]+[mjik]}} F(y), 

imix = | —[mjki]+[mjik]} F(y), 
where 

F(y) =x(y*—-y™) 
and 
(mijk |= Sp[(aen)(we;)(ae;)(aex) J, 
[mijk ]=([mkji]=4{ (enes)(e;ex) 
+ (€m@z)(€:€;) — (Cm ;)(€:€x)}. 

The remaining elements can be found from the 
following relations, which can be established by 


taking the complex conjugate of (16) and the 
corresponding equation for wim ;x. 


(18a) 
(18b) 


@mijk = — @kjim, 
Dimjk = — Wk jmi- 
The ‘c; are also antisymmetric in the same 
indices, and 
‘Cmisk = 1+ (fr —fi)/2, 
‘Cimik = 3 +fk/2—fi. 


Hence 


72 
> wc; = 16FE, (19a) 
1 


where 


E=) (eme,)(ejex)[1+2(fe—f;)]. (19b) 


ijk 


By (15a), (17) and (19) 


J= —2exE [ (y?—1)!y-*dy = —29xE/5 
1 
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and by (15) 
tt = (1/6002*) (ui/M)?(e*f1?/c*h*) (urc?/h) 


4 pi 
x f E*(fu fe)dfrdfe. (29) 
0 }-fi 


An approximate integration of this equation 
followed by a sum over both directions of polar. 
ization of the three photons and by an average 
over the polarization of the vector meson yields 
the final result 


T= (2.4/6002*) (ui/M)*(ef1?/cth*) (uic*/h). (21) 


5. DECAY OF A NEUTRAL VECTOR MESON INTo 
A NEUTRAL PSEUDOSCALAR MESON 
AND A PHOTON 


The density of states is now 
p=G(4r/h*)(u;*c) F,', (22) 


where the photon carries off the fraction F, of the 
rest energy of the meson. The reciprocal lifetime 
is now 


th = (1/1282) (uic?/h) (Cfo°f1?/Ch*) 


X (M/po)*(F,4/1—F,)J*, (23) 
where 


Sua f (Cwd-)[y*-1]}4ydy. (23a) 
1 


The total process requires three steps; for 
example, M,;>P++P-, Pt-P++y7:, P#+P- 
—+M>. Again there are six permutations, each of 
which appears twice in the way previously men- 
tioned. This case is slightly different from those 
already considered, because the pseudoscalar 
particle is created in a state of motion, and both 
non-relativistic and relativistic terms in the 
Hamiltonian must be considered. The relativistic 
numerator, except for a factor of 8, is the same 
as the numerator appearing in case (A). For 
example, 


N10 = Sp[A~(Baes)A(ae,)A, (5) ]. 
One finds 


m1,0=1[e: Xe, ly"(V—V,)/c, 
N07 =ife; Xe, ly" (V+Vo)/c, 
Nory = —iLe1 Xe, (vy +90") (Vo/c), 
nyo1 =ife, Xe, ly (V+ V,)/c, 
Moi =i: Xe, ly (V —Vo)/c, 
Nyi0= — tLe: Xe, |(y +977) (V,/c). 
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The leading term of }-nd;' caused by the 
relativistic part of m; is of order (u/M)*. Hence we 
consider only the non-relativistic term which is 
of the order (u/M)°.. A typical non-relativistic 


numerator is 
niyo = Sp[A~(Baes)A(ae,)A(ysako) J. 


The subscripts on the projection operators have 
heen omitted because only zero-order terms need 
be calculated. The non-vanishing terms of the 
spur contain @ as a factor three or five times. 
Those containing a five times are quadratic 
functions of V and must be averaged over all 
directions of V. The results are 


W170 = Wo1y = (82i/3)(y* —y")[e, Xe, |Ko, 
wry = (8xi/3)(y*+2y") [ei Xe, }ko. 


Again 
W170 =Woy1, Woly=Wy10, W10y = Wy01- 
Hence 
Ew (32i/y*) Le: Xe, |Ko. (24) 


In this case also, logarithmically divergent terms 
cancel. After squaring J, summing over the two 
possible polarizations of e,, and averaging over 
€;, one finds by (23) and (24) that the reciprocal 
lifetime is 


r= (16/274?) (uic?/h) (fers? /ch*) 
X(M/po)?Fy4(1— Fy), (25) 


and by the conservation laws 
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6. DECAY OF A CHARGED VECTOR MESON 


The decomposition (D) is very similar to (C). 
Although the electric charge present in (D) 
introduces some new complexity, it alters the 
y-lifetime only very slightly because the processes 
permitted by the charge are strongly forbidden 
by selection rules. Before describing the effects 
of the charge, however, we consider the decom- 
positions common to both neutral and charged 
mesons. For example, one such possibility is: 
M;t—Pt+(N), PtoPit+y7i, Pi#+(N)-Met, 
where (N) means antineutron. There are six 
permutations and the usual doubling. In the 
neutral case all of these twelve combinatorial 
possibilities are realizable, but in the charged 
case only six of them are permitted because (a) 
a neutron cannot emit a photon if its magnetic 
moment is neglected, and (b) reactions of the 
type P++M+t—P* are forbidden. The result is 
that the matrix element due to these processes 
for the charged case is only half its value for the 
neutral case, and the lifetime for either charged 
meson would be 


to =4re, (26) 


if the charged and neutral cases were otherwise 
the same. 

However, the charged meson, by interacting 
directly with the electromagnetic field, can 
decompose in other ways. It is convenient to 
divide these new processes which lead to decay 
into two classes according to whether they 
require two or three steps. The three-step pro- 


F, = (1/2)(1 — (wo/m)?]J. (25a) cesses are 
(a) M,t-M,t+y M,t->Pt+(N) Pt++(N)—-M,* 
(b) M,t-M,t+y 0-P-+N+M,;t M,*+P-+N-0 
(c) 0-P-+N+Myt Myt-Myt+y7 P-+M,t+N-0 (27) 
(d) 0-P-+N+M,t M,t-M,*+y7 P-+M,t*+N-—-0 
(e) 0-P-+N+Mot P-+Mit+N-0 Mot—Mott+y7 
(f) M,t+-—>P+(N) P+(N)—-M,* M,t—Mo* +7. 


It will now be shown that all processes (27) are 
forbidden. The nucleons are involved in only 
two of the three steps. In (c), (e), and (f) M, 
disappears and M, is created in a state of rest. 
The corresponding spur is Sp[A~(Baem) As]. This 
must be multiplied by the matrix element de- 
scribing the emission of a photon by Mo*; but 
since this spur vanishes, the corresponding 








process may be ignored. In (a), (b), and (d) 
both M, and M, move and the relevant spur is 


Sp{d-[Aovs+Bovsa [Aba +B,8ys0)} 
= Sp{d-[A oys\:(ArBa) + (Boysa)A,(Ai8a) 
+Aoysr:(BiBysa) + (Boyse)A:(BiByse) J}, 


where the A and B coefficients are given in Eq. 
(3). The first term vanishes identically; the 
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fourth term vanishes since By 1 B,. The second 
term reduces to a linear function of V,=[A; 
XB, ]-V; and since the energy denominator is 
an even function of V,, this term vanishes after 
the angular integration. The third term vanishes 
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for a similar reason, and hence all proc 

are forbidden. ae 
Finally the meson may interact simultaneoug| 

with the electromagnetic and mesonic fields in 

the following ways. 


(a) O0-P-+N+7+Mo+ P-+N+M1t-0 


(b) Mit+—P*++(N) 


(c) Mit-Pt+(N)+7 


(d) 0-M,t+P-+N 


In each case one step involves a pair, a photon, 
and a meson. The corresponding matrix elements 
follow from (2), where it is no longer permissible 
to replace D by V; they are (except for constants 
which do not concern us) 


(28a) 
(28b) 


(Un*G€,Um), 
Le, Xx en | (u n*BY5QUm) . 


The probabilities for the processes (a) and (b) 
vanish, as one can verify by direct computation 
or by the following observation. The interchange 
of the two members of the pair is equivalent to 
a reflection with respect to the origin of momen- 
tum space, since they have equal and opposite 
momenta. But this interchange reverses the sign 
of the probability amplitude ;* hence the proba- 
bility is an odd function in momentum space 
and vanishes when integrated over all directions. 

In (c) and (d) M,j is at rest so that it is only 
necessary to find 


nyo = Sp[A~(CiBysa)Ay(A ovst+Boyse) J, 


where C; is given in (28b). The symbol mo 
means that M;, is absorbed before Mo is emitted. 
It is also necessary to find m9. The results are 


nio=AoCi{[y"(V,/c) +9,1(V/c) ] 
+ (cpo/€o)(y*—9y")}, 


noi = AoCi{y—"(Vo/c) +901 (V/c) J 
+ (cpo/€o)(y — yo) }. 


On expanding mDio"'+m01Dor"' in powers of 


Pt++(N)—Mot+y 
P++(N)—-M,* (27) 
P-+N+My;-. 





y™ one finds that the leading term contains yt. 
hence the integral of }-”,;D,;~ over momentum 
space is convergent even though there js only 
one energy denominator for these processes. 
Furthermore, the leading term, which is of order 
(u/M)?, can be neglected in agreement with our 
previous neglect of the relativistic terms, also of 
the second order in (u/M), which appear jn 
case (C). Hence the lifetime in case (D) js 
given by Eq. (26). 


7. NUMERICAL RESULTS 


The significance of our numerical results js 
quite uncertain because of doubts concerning the 
correctness of meson theory and the applicability 
of quantum theory. For consistency, however, 
the following parameters were used in all cases: 
1/0 = 1.6, fe /he =0.05, fo/mo=fi/m. By Eqs. 
(12), (21), (25), and (26), the lifetimes given in 
Eq. (1) were computed. These parameters are 
given by Jauch and Hu® for the Schwinger 
mixture. The two constants ui/o and fo are 
chosen so as to give the correct *S binding energy 
and the correct 4S virtual level of the deuteron, 
while the mass yo is assumed to have the (low) 
value of 177 electron-masses determined by 
cosmic ray measurements. 

I would like to express my very sincere thanks 
to Professor J. R. Oppenheimer for having 
suggested this problem and discussed it with me. 


5 J. M. Jauch and N. Hu, Phys. Rev. 65, 289 (1944). 
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In an attempt to verify the theory of the long spark of Loeb and Meek, d.c. sparking poten- 
tials in a plane parallel gap in air have been measured as a function of plate separation and 


pressure at and below atmospheric pressure. 


Below atmospheric pressure, the sparking- 


potential curve as a function of gap length at constant pressure, displays a peculiar levelling-off 
at large gap lengths. This has been shown to be due to field distortion of a pressure-dependent 
type. This field distortion is, in large measure, caused by the walls of the large chamber used 
in this study. This field distortion prevented the verification of the Loeb-Meek theory. How- 
ever, there is considerable evidence for believing the results at atmospheric pressure to be 
valid. Meek’s equation represents the data accurately at atmospheric pressure. Water vapor 
causes a small but definite scattering in sparking potentials. 


INTRODUCTION 


HE streamer theory of sparking developed 

by Loeb and Meek! invokes processes 
known to occur in discharges at large values of 
the product of pressure and plate distance. From 
the theory sparking potentials may be calcu- 
lated. Other equations are either empirical repre- 
sentations of data such as the one given by 
Peek‘ or else they are based on mechanisms 
known to be inoperative at large values of the 
product of pressure and gap length. Loeb and 
Meek?* have extended the theory to include the 
long spark. They postulated that, at a given 
pressure, above a certain plate distance, 40, 
Meek’s equation does not apply. According to 
Loeb and Meek the slope of the sparking-poten- 
tial vs. gap-length curve should decrease with 
increasing gap length until 59 is reached. At this 
point the slope should increase and remain con- 
stant thereafter. The straight line above 4o, if 
extended, should pass through the origin and 
should lie below the values of the sparking poten- 
tial for gap lengths less than do. It was also 
postulated that this inflection in the sparking- 
potential curve would occur at shorter gap 
lengths the lower the pressure. On the basis of 


*Now at New York University, University Heights, 
New York. 
1]. M. Meek, Phys. Rev. 57, 722 (1940). 
usa” Loeb and J. M. Meek, J. App. Phys. 11, 438, 459 
+L. B. Loeb and J. M. Meek, The Mechanism of the 
Electric Spark (Stanford University Press, Stanford, 1941). 
‘F. W. Peek, Jr., Dielectric Phenomena in High-Voltage 
Engineering (McGraw-Hill Book Company, Inc., New 
York, 1929), p. 73. 
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the only data then or now available, Loeb and 
Meek calculated 9 for air to be 88, 15, and 3.3 cm 
at 1520, 760, and 380 mm of Hg, respectively. 
The theory is thus conveniently open to experi- 
mental test at pressures of the order of one-half 
an atmosphere. 

The theory of Loeb and Meek as applied to 
the long spark is admittedly semi-empirical, and 
the constants are not too well known. While the 
underlying principles may apply, the quantita- 
tive bases of the predictions may be off con- 
siderably. The validity of the choice of constants 
has been discussed recently.® 

Irrespective of the details of the calculation, 
which are somewhat in doubt, the very definite 
predictions require experimental verification, and 
toward this end, at the suggestion of Professor 
Loeb, the writer undertook to measure sparking 
potentials in air at half an atmosphere pressure. 
Unfortunately, it was just in the region of pres- 
sures and gap lengths which should prove of 
interest that the apparatus failed to function 
properly because of field distortion. In view, 
however, of the fact that very consistent spark- 
ing-potential data at atmospheric pressure were 
obtained, and that the instrumental difficulties 
encountered are of considerable interest, it is 
believed worth while to report the findings. 


APPARATUS AND EXPERIMENTAL PROCEDURE 
The chamber used was the one employed by 
Sanders® and Posin’ for measurements of the 


5 L. H. Fisher, Phys. Rev. 69, 530 (1946). 
* F. H. Sanders, Phys. Rev. 41, 667 (1932). 
7D. Q. Posin, Phys. Rev. 50, 650 (1936). 








first Townsend coefficient of air and nitrogen. 
This ionization chamber was constructed with 
the aid of a grant in support of Sanders’ re- 
search by the National Research Council. 
A diagram of the chamber as originally used by 
Sanders (and also for some time by the writer) 
is given in Fig. 1. For a detailed description of 
the chamber one may refer to Sanders’ original 
paper. For the purpose of the present discussion 
it will suffice to recall just a few features of the 
chamber. The shell consisted of a cylindrical 
brass casting 45 cm in height and 50 cm in 
diameter. Both electrodes were 20 cm in diam- 
eter. The lower plate was the high voltage 
electrode and was profiled as shown. This elec- 
trode rested on a porcelain insulator, J,; the 
insulator in turn rested on a circular metal 
plate P. The lower insulator, J2, was made of 
Pyrex glass and was filled with oil. No mercury 
was used in the chamber in the present work. 
With the arrangement shown in Fig. 1, at 
half an atmosphere corona from the high voltage 
electrode occurred at voltages too low to test 
the Loeb-Meek theory. Figure 2 shows the 
arrangement finally used. The high voltage elec- 
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Fic. 1. Ionization chamber as originally used by Sanders. 


424 LEON H. 





FISHER 


trode of Fig. 1 was replaced by another of equal 
diameter the edges of which had a radius of two 
centimeters. In addition, it was necessary to 
remove the supporting brass plate, P, to the 
bottom of the chamber. A new hollow porcelain 
insulator, J,’, was used which reached from the 
brass plate, P, to the high voltage electrode. 
This increased the distance from P to the high 
voltage electrode from 10 to 20 cm. A new ojj. 
filled insulator, J:’, entered the chamber as 
before, but reached all the way to the lower 
electrode. With this arrangement, it was possible 
to reach 60 kv at 380 mm of Hg before corona 
appeared. 

The plates were levelled to within 0.002 inch 
outside the chamber by means of a telescopic 
plug gauge. Whether the plates remained level 
after the electrode system was lowered is yp. 
known since mechanical stresses are set up. The 
electrodes were polished with emery paper and 
rouge, and the chamber surfaces were carefully 
washed with ethyl alcohol before being closed 
for a series of measurements. Fiducial marks 
were placed on both electrodes and gap distances 
were measured by means of these marks with a 
cathetometer. Gap separations were read to a 
few thousandths of a centimeter. 

The d.c. potential source was of the same type 
as the one used by Sanders; it was stabilized by 
an electronic stabilizer of the saturable-core 
reactor type. In parallel with the spark gap 
were wire-wound Taylor resistors of 100 meg- 
ohms. These resistors are manufactured by Shall- 
cross and carry a one percent accuracy rating 
when used under the present operating condi- 
tions.* The voltage was read by potentiometer 
across 1000 ohms set on a Leeds and Northrup 
dial box in series with the Taylor resistors. The 
accuracy of the voltage measurements was there- 
fore determined by the Taylor resistors and is 
of the order of one percent. 

A current-limiting resistor of 0.1 megohm was 
placed in series with the chamber in order to 
prevent pitting of the plates and also to avoid 
excessive production of oxides of nitrogen. The 
Taylor resistors were placed directly in parallel 
with the chamber. Leakage across the current- 
limiting resistor and across the oil-filled insulator 


*L. S. Taylor, Bur. of Stand. J. Research 5, 609 (1930). 
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‘troduced no uncertainty in the measurements. 
The chamber was grounded during the meas- 
urements. 

The central portion of the cathode was illumi- 
nated by ultraviolet light in order to prevent 
appreciable time lags, and also to induce the 
spark to pass in the uniform part of the field. 
During the entire study, however, ultraviolet 
light had no effect on the sparking potentials, 
and in no way influenced the position of the 
spark. This situation must have been due to the 
presence of adequate initiating electrons in the 
chamber at all times. 

Air from outside the laboratory was admitted 
through a capillary, the air passing over two 
traps of aged liquid air, and one of phosphorus 
pentoxide. At no time was liquid oxygen observed 
in the trap. The pressure of the gas was read on 
a mercury manometer to one millimeter. All 
pressures have been corrected to 22°C. 

To measure a sparking potential, the voltage 
was raised in steps of 100 volts, the electrodes 
being in contact with the power supply at all 
times. After each setting, a 30-second interval 
elapsed before the voltage was increased again. 
The time lag was almost always so short that 
the 30-second wait was unnecessary. The voltage 
taken as the sparking potential was the reading 
at the passage of a spark. For a given reading, 
two measurements were taken with ultraviolet 
light incident on the cathode, and one was taken 
with no irradiation of the gap. 

Previous observers*!® working with air have 
reported the lowering of sparking potentials 
because of the formation of oxides of nitrogen. 
In this study, no such effect was observed, the 
chamber being large enough to prevent impurities 
from building up to an appreciable concentra- 
tion. The current-limiting resistor also helped to 
reduce the amount of chemical action during a 
spark. 

Sparking potentials in the closed chamber in 
dry air for any one filling were reproducible to 
100 volts regardless of the values of the sparking 
potential (10 to 60 kv). From one filling to 
another, the results scattered by about 300 volts. 
No large spread in individual measurements, as 
is often reported, was observed. This fact is 


*K. S. Fitzsimmons, Phys. Rev. 61, 175 (1942). 
”W. R. Haseltine, Phys. Rev. 58, 188 (1940). 
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Fic. 2. Ionization chamber as used in the present work. 


probably due to the exclusion of dust particles 
from the chamber, as well as to the presence of 
adequate initiating electrons. The exclusion of 
water vapor is believed to have increased the 
reproducibility of the measurements. On closing 
the chamber, or even after a new filling, sparks 
occurred at low voltages, but after five or ten 
sparks passed, the sparking potential became 
constant. 


EXPERIMEN1aL RESULTS AND DISCUSSION 


Table | represents the results of a run, series |, 
at atmospheric pressure with varying gap length. 
The results are plotted in Fig. 3. Series I yields 
30.1 kv as the sparking potential for a one- 
centimeter gap at atmospheric pressure at 22°C. 

Table II gives the values quoted by Schumann" 
as a compendium of the results of a large number 
of workers, and also the careful observations of 
Spath.” 

Values for corresponding values of pé from 
series I are also given in Table I]. Spath’s values 

" W. O. Schumann, Elektrische Durchbruchfeldstdirke von 


Gasen (Julius Springer, Berlin, 1923), p. 26. 
# Reference 11, p. 173. 
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Fic. 3. Sparking potential vs. pé for series | 
(near atmospheric pressure). 


fall consistently below those given by Schumann, 
and this difference has been ascribed to Spath’s 
use of ultraviolet light. The present results above 
0.60 cm fall below those of Spath, but the agree- 
ment below 0.60 cm is excellent. Measurements 
made by Haseltine’ give V, as 15.4 kv at 
225 = 324 mm Xcm as compared to 15.5 kv from 
the present data. 

In Meek’s equation, a constant, K, appears 
which is the ratio of the field due to the positive 
ion space charge to that of the applied field. 
It is assumed that at the sparking voltage K is 
a constant independent of pressure and gap 
length. That K is a constant independent of 
pressure has been demonstrated for streamer 
formation in hydrogen in inhomogeneous fields." 
It is of interest to apply Meek’s equation to 
consistent data in order to evaluate K. As a 
matter of fact, the reverse procedure may be an 
extremely valuable method of estimating the 
consistency of a set of data. As is commonly 


TABLE I. Series I measurements of pm potential at 
constant pressure [p(22°C) = 756 mm Hg ] and varying gap 











length. 

6 p25 Vv; In.K ad 
1.985 cm 1501 mm Xcm 54.7 kv —7.6 $3.4 
1.801 1362 50.1 —7.7 12.1 
1.606 1214 45.4 —7.2 12.4 
1.425 1077 40.7 —7.9 11.6 
1.207 912 35.4 —7.2 12.0 
1.022 773 30.5 —7.5 11.5 
0.831 628 25.5 —7.0 11.5 
0.599 453 19.3 —7.0 11.5 
0.401 303 13.8 —6.8 11.0 








%L. H. Fisher and G. L. Weissler, Phys. Rev. 66, 95 
(1944). ‘ 
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realized, V, is an extremely insensitive function 
of K; therefore, K is an extremely sensitive 
function of the value taken for V,. For example 
if it is assumed that a one-centimeter gap in 
air at atmospheric pressure has the sparking 
potentials given in Table III, the values of K 
given there result. Thus, for these conditions a 
range of two kilovolts introduces a change in K 
of a factor of 1000; a one-kilovolt difference 
introduces a factor of 30, and a difference of 
0.1 kv introduces a change of 20 percent. It ig 
necessary, therefore, to be cautious in interpret. 
ing calculated values of K. In Tables I and Jj 
are listed values of the natural logarithm of x 
calculated from series I and from the measure. 
ments of Spath. It is seen that the present 
measurements give values of K constant to 
within a factor of two, while those of Spath give 
values of K which vary by a factor of 60. Thus 
one value of K can be used to represent the data 
of series | within the limit of experimental error; 
this is not true for the measurements of Spath. 
In Tables I and II are also listed calculated 
values of aé at sparking, and it is seen that the 
condition aé=constant represents the data of 
series I in a satisfactory way. Again this is not 
true for Spath’s measurements. 

The small absolute values of K calculated 
from Spath’s and the present data are disturbing, 
since it is difficult to see how photo-electrons 
formed in the gas would have much tendency to 
feed into the positive space charge. Since, how- 
ever, no great accuracy is claimed for the con- 
stants in Meek’s equation, errors in this might 
account for the low value observed. 

The fact that series I gives lower values of 
the sparking potential than are usually reported 
is in itself not too disturbing. Most of the former 


TABLE II. Comparison of present results with those quoted 
by Schumann** and with those of Spath.*** 

















Ve V. V; (This In.K ab 
é p25 (Schumann) (Spath) study) (Spath) (Spath) 
0.40 cm 306 mm Xcm 14.4 kv 13.8 kv 14.1 kv —73 106 
0.50 382 17.4 16.8 16.9 —6.5 11.6 
0.60 459 20.3 19.8 19.6 —5.9 124 
0.70 536 23.2 22.75 22.3 —5.2 182 
0.80 612 26.1 25.6 24.9 —48 138 
0.90 688 28.9 28.6 27.5 —3.8 149 
1.00 765 31.7 31.55 30.1 —3.2 15.6 
2.00 1530 59.6 55.7 
** See reference 11. 


*** See reference 12. 
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Tasie III. Values of the constant K of Meek’s equation for 
s various sparking potentials. 





! 





Sr 
V 33.0kv 32.1 32.0 31.0 
peestealated) 20 0.6 0.5 0.02 








work used a.c. potentials. Loeb,'* among others, 
has pointed out that a.c. sparking-potential 
measurements may be expected to be slightly 
higher than d.c. measurements. The peak value 
of an a.c. voltage is applied only instantaneously, 
and because of a statistical time lag, the gap may 
not break down until a value of the peak voltage 
is impressed which is somewhat greater than the 
true sparking potential. Reukema’® found that 
the use of ultraviolet light reduced the average 
observed a.c. sparking potential by 3.5 percent. 

Measurements of sparking potentials at half 
an atmosphere and with varying plate separation, 
series IV, are plotted in Fig. 4. It is seen there 
that the sparking-potential curve at half an 
atmosphere levels off at large values of the plate 
separation. In order to investigate this levelling- 
off, runs were made at pressures of 50 and 
25 cm. Series V was made at 500 mm, and series 
VI was made at 261 mm. The results, including 
those of series I and IV, are plotted in Fig. 5. 
It is evident from Fig. 5, that for p5>800 mm 
Xcm, Paschen’s law is not obeyed, and some 
kind of field distortion comes into play, at least 
for pressures below one atmosphere. The data 
for series 1V, V, and VI above p5= 800 therefore 
have no significance. 

In order to understand the nature of the field 
distortion, a cylindrical piece of sheet metal was 
placed in the chamber and was grounded. This 
cylinder reduced the effective diameter of the 
chamber by two inches, and measurements were 
made again at pressures of one atmosphere, 50 
cm, 38 cm, and 25 cm. The measurements at 
atmospheric pressure and at 50-cm pressure did 
not differ from the values obtained with the 
original chamber. Therefore, for these high pres- 
sures, the effect of reducing the size of the 
chamber by two inches was unimportant. This, 
however, was not the case for the lower pressures. 


4L. B. Loeb, Fundamental Processes of Electrical Dis- 
charge in Gases (John Wiley and Sons, Inc., New York, 
1939), p. 550. 

*L. E. Reukema, Trans. A.I.E.E. 47, 38 (1928). 
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Fic. 4. Sparking potential vs. 6 for series 1V 
(near half an atmosphere pressure). 


The data for the 38-cm (series VI1) and 25-cm 
(series VIII) pressure runs with the reduced 
size of the chamber are plotted. Series IV and 
VII are plotted in Fig. 6, and series VI and VIII 
are plotted in Fig. 7. The levelling-off of the 
sparking potential vs. plate-separation curves is 
more pronounced with the smaller chamber than 
with the original chamber at pressures in the 
neighborhood of 380 mm and below. It is clear, 
therefore, that there is at low pressures and 
large gap lengths a field distortion imposed by 
the walls of the chamber, and that the question 
of the behavior of the sparking potential vs. gap- 
length curve at half an atmosphere must wait 
until a larger chamber is available. Clearly, a 
chamber of much greater diameter will be 
needed in order to study sparking potentials at 
low pressures and large gap distances. 

The question now arises as to the reliability 
of series 1. Below p5= 800 mm Xcm, there seems 
to be little doubt that the values are reliable. 
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Fic. 5. Sparking potential vs. pé at four different pressures. 
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Fic. 6. Sparking potential vs. pd at half an atmosphere 
pressure with original chamber (series IV) and with slightly 
smaller chamber (series VII). 


Above this value of ~6, one can only say that 
consistent values of K are obtained. There is 
thus some reason for believing series I to be 
valid in its entirety. 

Various other attempts were made to discover 
in more detail how the geometry of the chamber 
affected the measurements. The polarity of the 
electrodes was reversed, and a complete set of 
measurements was taken at various pressures. 
The results were, within the experimental error, 
identical with the values obtained with the 
usual polarity. 

To determine the effect of the small radius of 
curvature of the grounded electrode, the elec- 
trode that was used by Sanders as a cathode 
(the lower plate in Fig. 1) was inserted as the 
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Fic. 7. Sparking potential vs. 6 at one-quarter of an 
atmosphere pressure with original chamber (series VI) and 
with slightly smaller chamber (series VIII). 
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grounded electrode in the arrangement of Fig, 2 
In spite of the decreased curvature of the a: 
electrode, the results were again unchanged 
This result indicates that the contour of the low 
voltage electrode made no contribution to the 
field distortion. Whether the contour of the high 
voltage electrode added to the field distortion is 
at present unknown. A model study is now in 
progress. 

The one fact which seems to require the high 
voltage plate to have caused some distortion js 
that the sparking potentials with the reduced 
chamber at a pressure of 50 cm were the same 
as for the large chamber at the same pressure, 
and yet the sparking-potential curves at pres. 
sures of 50 cm and 76 cm do not coincide when 
plotted as a function of pé. 

The electrode system of Fig. 2 was lifted into 
the open room (the surrounding metal walls were 
effectively removed) and sparking potentials 
were determined ; the electrode system was then 
lowered into position so that the electrodes were 
now surrounded by the metal walls, and measure. 
ments were made again. The two sets of observa- 
tions (at atmospheric pressure) were identical, 
The complete removal of the walls, just as 
decreasing the diameter of the chamber, had no 
influence on the sparking potential at atmos- 
pheric pressure over the range of plate separa- 
tions studied. 

In the experiment described in the preceding 
paragraph it was found that the sparking poten- 
tials both in the open and closed system scattered 
considerably ; this probably was caused by the 
presence of dust or water vapor. The scattering 
persisted even after considerable sparking and 
was not removed by ultraviolet light. When 
clean dry air was admitted in the normal manner, 
the sparking potentials were again reproducible. 
Water vapor was then allowed to diffuse into 
the chamber, and it was found that a small but 
definite spread of 300 to 500 volts is introduced 
in the value of the sparking potential. 

The author wishes to express his gratitude to 
Professor Loeb who suggested the problem and 
under whose direction it was carried out. 
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Fig. 2. The y-Radiation of Na*‘ and the Energy-Level Scheme of Mg** 
he new M. L. WiEDENBECK 
langed, Department of Physics, University of Michigan, Ann Arbor, Michigan 
he low (Received July 23, 1947) 
to the 
~ high INCE Na* can be produced easily and be- beta-spectrograph study of the secondary elec- 
ton be cause it emits high energy y-radiation it has trons ejected from a body irradiated with the 
now in long been of interest from the applied standpoint -rays in question. Such work has been carried 

. as well as from the standpoint of the structure of out by Itoh,‘ Elliott, Deutsch and Roberts,' 
7 high the levels giving rise to the radiation. Over the Mandeville,* and most recently by Siegbahn.' 
—s past several years the radiation has been in- Their results all indicate that only two y-rays are 
apt vestigated by many workers by various methods emitted by the Mg* nucleus and these are of 
— and in many cases the results have been incon- ‘“‘nearly’’ equal intensity. The values given are 
esure, | sistent. 1.38 and 2.76 Mev. 
: —_ The sodium-8-spectrum has been most accu- The most simple arrangements of the levels of 
— rately studied by Siegbahn' and by Lawson? and Mg™ would occur if the y-rays were in cascade. 
d into has been found to be simple. Absorption experi- However Guthrie and Sachs” have pointed out 
a walk ments by Langer, Mitchell, and McDaniel’ in that such an arrangement would lead to a mass of 
ential 1939 have shown that the 8 —y-coincidence rate Na* which is ~1.5 Mev higher than that calcu- 
pn is independent of the energy of the B-radiation, lated by Barkas* from a study of the fine 
2 wale thus showing that a single 8-transition takes structure of the mass defect curve. They further 
oni place. The upper limit of the 8-spectrum is pointed out that this difficulty could be elimi- 
ari 1.390 Mev. nated by assuming the level scheme shown in 
ntical. The history of the y-radiation following the Fig. 1b with a branching ratio of 2:1. gt} 
te 8-decay of Na**is somewhat more complex. This Absorption measurements of the y—y-coinci- 
ad no has, however, been described in considerable dences were carried out by Cook, Jurney, and 
‘tmos- detail by Siegbahn.! It would appear that the Langer.* These indicate that the 2.8-Mev radia- 
para- most accurate results are those obtained by a_ tion is involved in the coincidence and hence 

supports level scheme of Fig. 1a. 
seding _ Recently Barker'® has measured the average 
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de to e b ‘ |. Itoh, Proc. Phys. Math. Soc. Japan 23, 605 (1941). 
1 and 5 L. G. Elliott, M. Deutsch, and A. Roberts, Phys. Rev. 


Fic. 1. ‘‘a” shows the two possible cascade energy-level 
diagrams of the decay of Na*™. “‘b”’ indicates the parallel 
decay proposed by Sachs. 


1K. Siegbahn, Phys. Rev. 70, 127 (1946). 

?]. L. Lawson, Phys. Rev. 56, 131 (1939). 

*L. M. Langer, A. C. G. Mitchell, and P. W. McDaniel, 
Phys. Rev. 56, 962 (1939). 
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* W. H. Barkas, Phys. Rev. 55, 691 (1939). 

*C. S. Cook, E. Jurney, and L. M. Langer, Phys. Rev. 
70, 985 (1946). 
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Fic. 3. Counting arrangement used in determining N7,7, 
and N7,8. The source, S, is fixed rigidly to counter 71. 


y-energy per 6-particle for Na* by means of 
ionization methods. His results give 4.17 Mev 
(+10 percent) per 8-particle. This would also 
substantiate the level scheme indicated in Fig. 1a. 

Tocheck the simplicity of the Na* 8-spectrum, 
a determination of the 8—-y-coincidence rate was 
made as a function of B-absorber thickness. The 
minimum thickness of absorber was equivalent to 
approximately 4 mg/cm? of aluminum. The curve, 
Fig. 2, is seen to be independent of the 6-energy 
and again attests to the simplicity of the B- 
spectrum. 

An additional method which can be used to 
distinguish between the level schemes of Figs. 1a 
and 1b lies in the y—y- and 8—y7-coincidence 
rates. If we have two counters designated by 7; 
and 72 (Fig. 3), each having an efficiency ¢1.4 for 
the 1.4-Mev 7-radiation and an efficiency ¢2.3 for 
the 2.8-Mev y-ray, then 


Ny, 7, = (207, €1.4€2.8/€1.4+ €2.8) N7,, (1) 


for level scheme of Fig. 1a. Nv,7, being the y—7- 
coincidence rate, oy, the fraction of the total 
radiation intercepted by counter 71, and Ny, the 
counting rate of counter y2. Now if the counter 72 
is moved farther from the source and converted 
to a #-counter, the 8—vy-coincidence rate for 
scheme of Fig. 1a is: 


Ney, =o7,(€1.4+ €2.8) Ne. (2) 
Dividing Eq. (1) by Eq. (2), and remembering 
that the source, S, is fixed rigidly to the counter 
v1, we have: 
(Nv,7,/N7,)/(Nv,8/Ne) 
_ 2€1.4€2.8/ (€1.4+ €2.8)°. (3) 
Similarly for level scheme of Fig. 1b: we obtain 


TABLE I. Comparison of y—~y and 8—vy coincidence rates. 











N7172/N 72 Total relative 
a efficiencies 
Radiator Nv18/Ng (brass =1) 
Brass 0.525 1 
Cu 0.525 1 
Al 0.51 1.2 
Fe 0.502 1.06 
Pb 0.517 0.95 
Lucite 0.51 1.3 
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Fic. 4. Ratio, R, of Nri2/N,/N7,8/Ne asa function of 
the ratio of the counter efficiencies for the 2.8- and 1.4-Mey 
y-rays. 


the ratio 
(Nyv,7,/Ny,)/(Nv,6/ Ne) 
= 3(€1.4)?/2(€1.4+¢€2.8)". (4) 

It should be noted that these equations are 
derived on the assumption that the efficiency of 
the counters is low, actually they are ~one 
percent efficient. 

Equations (3) and (4) are plotted in Fig. 4. It 
is clear that scheme la gives rise to a ratio 
(Nyv,7,/Nv,)/(N7,8/Ne) nearly equal to 0.5 over 
a very wide range of efficiencies while scheme 1b 
gives a rapid variation in the region of interest. 

The counters were of brass construction witha 
4-mil Cellophane end window. To convert the 
counter for y-counting, a small cylindrical plug 
was inserted in front of the window. 

Table I gives the results of six determinations 
made with various radiators in front of the thin- 
windowed counter ($-mil Cellophane). The total 
relative y-efficiencies of the counters are also 
included in the table. It is seen that over this 
wide range of radiators the ratio remains es 
sentially constant near the value 0.5. The experi- 
ment with the lead radiator was repeated with a 
greater thickness of lead to shift artificially the 
efficiency ratio (€2.s/€1.4) to a higher value. The 
counting rate with such an arrangement is 
low but within the statistical error the ratio 
(Ny,7,/N7,)/(N7,8/Ne) remains in the region of 
0.5. It must, therefore, be concluded that the 
cascade level diagram of Fig. 1a is the correct one. 

We wish to express our thanks to Professor 
Uhlenbeck for his discussions of this problem. 








i 


2 Se |] 


pre 


qu 


usi 
19: 
ant 
in 1 
ene 








20 


INCtion of 
| 1.4-Mey 


» (4) 
ONS are 
ency of 
 ~one 


ig. 4. It 
a ratio 
.5 Over 
eme Ib 
rest. 

witha 
ert the 
al plug 


nations 
e thin- 
e total 
re also 
er this 
ins es 
experi- 
with a 
lly the 
e. The 
ent is 
» ratio 
rion of 
at the 
ct one. 
fessor 


n. 








pHYSICAL REVIEW 


Letters to the Editor 
UBLICATION of brief reports of important discoveries 
in physics may be secured by addressing them to this 

department. The closing date for this department is five weeks 
prior to the date of issue. No proof will be sent to the authors. 
The Board of Editors does not hold itself responsible for the 


opinions expressed by the correspondents. Communications 
should not exceed 600 words in length. 











The Yield of Nuclei Formed by Fission 


A. J. DEMPSTER 
Argonne National Laboratory, Chicago, Illinots 
July 9, 1947 


HE recently published yield curve' for the nuclei 

formed in fission emphasizes the lack of any explana- 
tion for this striking feature of the fission phenomenon. 
Calculations of the potential energy associated with de- 
formations of the liquid-drop model of the nucleus** have 
met with difficulties in attempting to predict any of the 
details of the fission process. The purpose of this note is to 
point out that a physical condition is imposed on the fission 
process by the consideration that, since the pair of fission 
products are structures of known energy content, they 
must be formed with a definite potential energy. This 
condition is automatically satisfied in the liquid-drop model 
of the nucleus, at least to the extent that this model 
accounts for the masses of the nuclei, but its formulation as 
a quantitative physical condition independent of any 
theoretical interpretation of mass defects is of interest as it 
suggests a reason why some types of fission are more 
probable than others. 

The curve marked E in Fig. 1 is the kinetic energy ac- 
quired by the different pairs of fission products, whose 
masses are plotted as abscissae. It is computed from the 
mass changes between uranium and the fission nuclei, 
using the packing-fraction curve‘ given by the writer in 
1938, and allowing 22 Mev for the energy of excitation,® 
and the equivalent energy of eight extra electrons that are 
in the fission products when formed. E is also the potential 
energy at the instant of separation of the two fission nuclei. 
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Fic. 1. Curve E gives the kinetic energy of fission products, or the 
potential energy at separation. Curve d gives the separation of centers 
of spherical fission products that would account for E. Ri+Re gives 
the theoretical sum of the radii of the fission products. 
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The most probable type of fission gives the masses 95 and 
139, indicated by the arrows, and not two equal masses, 
even though this latter process would release a greater 
amount of energy. The curve d is the computed distance 
between the centers of the fission products, supposed 
spherical, at which the potential energy, Z(92 —Z)e*/d, due 
to Coulomb repulsion, is equal to E. The value of d comes 
out to be larger than the sum of the radii Ri +R: given by 
the liquid-drop theory of atomic nuclei, using the formula :* 
R=r,A}, with ro= 1.48 X 10- cm. Even smaller values for 
ro have been proposed.? The sum of the radii of the two 
fission products, calculated from this formula, is given by 
the broken line at the bottom of Fig. 1. There appears to be 
a difficulty, analogous to the paradox met with in alpha-ray 
emission, in that the fission particles must have no kinetic 
energy at a distance between centers (d), greater than the 
sum of their radii. It may, however, be possible to increase 
the theoretical values of the radii enough to make the lower 
curve cut across the curve d. 

The curves, however, suggest possible reasons for the 
different fission yields. For the cast of large differences in 
the masses, d must be large and E small, and only an 
unusual manner of separation would give the small po- 
tential energy permitted. With the fission masses 95 and 
139, a normal separation can be considered as leaving the 
products with the required potential energy. In the case of 
even division, on the other hand, only unusual separation 
processes would be able to provide the large potential 
energy in the fission products required by the curve E. 

1 The Plutonium Project, Rev. Mod. Phys. 18, 539 (1946). 

?R. D. Present, F. Reines, and J. K. Knipp, Phys. Rev. 70, 557 
(ON. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 

4A. J. Dempster, Phys. Rev. 53, 870 (1938). 

5K. Way and E. Wigner, Phys. Rev. 70, 116 (1946) 


*H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8. 97 (1936). 
7 E, Feenberg, Phys. Rev. 55, 504 (1937). 





Geiger-Mueller Counter Detection of Light Radi- 
ation from the Paths of High 
Energy Particles 


Paut B,. Weisz AND B. L. ANDERSON 


Socony-Vacuum Research and Development Department, 
Paulsboro, New Jersey 


July 14, 1947 


ETECTION of Cerenkov radiation has been at- 
tempted by use of photographic methods! and photo- 
multipliers.* Furry* recently pointed out the possible value 
of Geiger-Mueller counter tubes of high light sensitivity to 
cosmic-ray studies. 

An experiment was designed in an attempt to detect 
Cerenkov radiation accompanying the passage of cosmic- 
ray particles and of radiation from a natural radioactive 
source by means of light-sensitive Geiger- Mueller counters. 

The counters used had been developed in connection 
with other work‘ and had their spectral response in the 
ultraviolet, between 2000A and 3000A. They consisted of a 
cathode cylinder of wire screen of 50 percent coverage, 
covered with a 0.005-in. thick layer of electroplated gold, 
and contained in a Corning 9741 glass envelope capable of 
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Fic. 1, a, Arrangement of the cosmic-ray experiment. b. Arrangement 
using thorium-nitrate salt with its decay products. 


transmission in that spectral range. The quantum efficiency 
of the surface can be affected within wide limits by the 
nature of the degassing process and baking atmosphere 
before filling. 

The experiment is based on an arrangement in which the 
counter is immersed into a liquid medium which itself is 
transparent to the light radiation in the 2000-3000A range 
to which the counter is sensitive. The counting rate N, of 
the counter should then be due to high energy radiation 
penetrating the counter itself, and in addition to photons 
emitted by high energy particles in the liquid medium 
outside of the counter. In the experiment a second counting 
rate No is obtained after adding a small amount of strongly 
ultraviolet absorbing chemical, which will neither disturb 
the geometry nor the density of the liquid by an appreciable 
amount, but should eliminate those counts obtained from 
light radiation emitted in the liquid. Effectively, the liquid 
medium acts to increase the ‘effective counting volume”’ of 
the Geiger counter while the liquid is transparent to the 
light to which the tube is sensitive. 

Figure la shows the arrangement of the cosmic-ray ex- 
periment (A); Fig. 1b for the experiment using thorium- 
nitrate salt (with its decay products) as emitter (B) and 
(C). In the latter experiment it was attempted to stop the 
beta-radiation from the thorium within the pure water in 
order to keep the count from direct beta-radiation as low as 
possible. The liquid medium was always pure water (its 
transmission was checked to give an optical density of 
D=0.005 at 2500A increasing to D = 0.034 at 2200A and no 
measurable absorption between 2650-3000A). The added 
chemical to make it opaque was 1 percent (by weight) of 
hydroquinone in (A ), and a variable amount (as indicated ) 
of potassium-chromaite in (B) and (C). 

Table I shows the results of the experiments. 

Percentages of chemical are given as weight-percents of 


TABLE I. Results of experiments on the detection of Cerenko ates 
from cosmic rays and natural radioactive decay products ation 








— 
(A) Cosmic rays—counter K-5 


N:(counts/min.) No(counts/mi 
Run pure H:O HO+1% Codon. 
1 23225 21045 
2 254245 22446 
3 263 +7 22445 
Average 249 219 





(B) Thorium nitrate—counter K-9 





N:(counts/sec.) No(counts/sec.) 
pure H:0 H20+0.1% KeCr0, 
1 24721 246+1 
2 250+1 247+1 
3 24931 24341 
4 247+1 243 +1 
5 248+1 245+1 
Average 248.2 244.8 
(C) Thorium nitrate—counter K-9 
Liquid N(counts/sec.) 
pers H:0 148.7 +0.8 148.7 
20 +0.5% KeCrO, 146.7 +0.8 ) 
H20+1.0% K2CrOx 146.2 +0.8 146.3 (Ay. 
H:0+1.5% KsCrOx 146.3 +0.8 3 (Average) 
H20+2.0% KeCrOx 146.0 +0.8 








chemical added to the pure water. In experiment (C), 
which was done with a smaller quantity of thorium salt, 
additional amounts of K:CrO, were added to see whether 
the counting-rate reduction would continue. This was obvi- 
ously not the case, as had to be expected if light radiation 
was the source of eliminated counts. 

The counter K-5 was unusually sensitive, while K-9 had 
the normal sensitivity which we had been accustomed to 
obtain from our surface preparations, and which has been 
estimated from their response to blackbody radiation to be 
of the order of 2 to 8X 10~* counts/quantum. 

Calculation shows that the cosmic-ray counting rate 
experiment could be expected to show a 10 percent effect 
with a Geiger-Mueller counter quantum efficiency of about 
6010-4 counts/quantum. Easy detection of Cerenkov- 
type light radiation accompanying cosmic-ray particles in 
an experiment as described is indicated by these calcula- 
tions at quantum-efficiencies of the order of 5X 10-* or 
better. In the thorium éxperiment a larger effect might at 
first be expected due to the large number of primary 
particles available. However, the very small effective path 
in emitting Cerenkov radiation before their velocity be- 
comes so low that Bn >1 is not satisfied any more, and the 
smaller sensitivity of counter K-9 makes these results 
appear reasonable. 

While such experiments appear relatively simple, it is 
felt that some attention should be diverted to the possible 
contribution of electronically excited molecular states in 
the stopping medium induced by the passage of the high 
energy particles. 

Work using more promising photoelectric surfaces such 
as zinc is continuing. 

1P, A. Cerenkov, Phys. Rev. 54, 499 (1939); G. B. Collins and 
V. G. Reiling, ibid. 499 (1939). 

21. A. Getting, Phys. Rev. 71, 123 (1947); R. H. Dicke, Phys. Rev. 
71, 737 (1947). 

3W. H. Furry, Phys. Rev. 72, 171 (A) (1947). 

4P. B. Weisz, Electronics, July 1946. 
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On Explicit Solutions of the Equations for Steady 
Compressible Flow 
Bruce L. Hicks 


Ballistic Research Laboratories, Aberdeen Proving Ground, Maryland 
July 21, 1947 


N classical aerodynamics the steady flow of incom- 
I pressible, inviscid fluids can often be described ex- 
plicitly in terms of functions of a complex variable. There 
are, however, few accurate solutions of the equations for 
compressible flow that can be obtained in closed form. In 
view of the importance of obtaining more exact solutions, 
especially in the transonic region, the following outline of 
a new method of approach is submitted (see related work 
of Munk'). The flows considered are adiabatic but not 
necessarily isentropic and iso-energetic. 

For uniplanar irrotational flow, the compressible flow 
equations can be written 


aN/dy =A1(N)dd/d¢, (1) 
aN/de = —[A2x(N) J '08/dy, (2) 


where the velocity vector V=[g(N)RT}N; |N|=N 
=|V|/Ce(N)RT}; 8 is the angle of inclination of the 
vector V to the x axis; and gand y are potential and stream 
functions for the vector N, defined by the equations 
N=Ve=VXky. To a given choice of the function g(N), 
there correspond flows with a given physical characteristic, 
namely, a relation between N and the pressure p on curves 
normal to the streamlines. Choice of g(N) also fixes the 
form of the functions A;, A2 and all other functions of N 
or XN which will be introduced subsequently. 

Except when g(N) «(1—N?)~, irrotationality of N re- 
quires,’ if the equation of motion is to be integrable, that 
(i) the stagnation pressure ~; vary according to the 


equations 
logp: = ¥(¥) = TU(N) — ¥(¢) (3) 


(where ® and W are arbitrary), and is constant (between 
shocks) on streamlines; (ii) V?—« £(N)®’(¢). If g(N) 
«(1—N*)", the integrability condition yields, instead 
of (3), the simple result Vp, =0. The second of Eqs. (3) can, 
however, be assumed here also, in order that explicit in- 
tegration of the flow equations can be achieved for this 
important case. It is possible, with other choices of g(N), 
that the detailed analysis may be considerably simplified 
as it is known to be for diabatic flow.” 

We now assume that (1) and (2) are integrable; then if 
N(¢, ¥) can be found by some means, # is calculable from 


€¥ ON ¢. vo _aN 
j= ~_ A2(N Bat fey [Ai(N)] Free (4) 


Equation (3) together with assumed integrability of (1) 
and (2), demands that the functional-differential equation, 
involving the functions @(¢) and ¥(y), must be satisfied. 
p2( SU) F’ (©) + 2yo’ (TN) F(S) 
+ui(SU)P'(¥)+2yi'(T)P(¥)=0, (5) 
where 
F(#)=(d@/dg)+; P(¥)=(dv/dy). 


After repeated integration with respect to @ and ¥, 
using (3), Eq. (5) can be reduced to a functional equation 
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in ®(¢) and ¥(y), which we shall schematically represent by 
F[®, W, ci(e), ki(e)]=0, (6) 


in which ¢;(¢) and k;(¢) are functions introduced in the 
partial integrations. 

If any functions (¢) and ¥(y) can be found that reduce 
(6) to an identity, then (3) and (4) immediately yield 
explicit, closed expressions for N and 3. Further integration 
would yield x and y as functions of ¢ and ¥. Simple radial 
and vortex flows correspond to choosing W and @, in turn, 
to be constant. No other exact solutions of the functional 
Eqs. (5) and (6) are known at present, although all pre- 
viously known compressible flows have not been tested 
for their possible conformity to these equations. 

The utility of the method, as compared to others, will, 
of course, depend on whether exact or approximate solu- 
tions of (5) and (6) can be found that yield flow patterns 
of some physical interest. There are indications that such 
solutions can be found in the transonic region. It is noted 
that the description N=V¢ applies to fields other than 
those of constant entropy and stagnation pressure, and 
may, therefore, make possible investigations of compressible 
flows down-stream from regions of heating and viscous 
dissipation. 

This investigation was begun at the Cleveland laboratory 
of the N.A.C.A. in 1945. The assistance of the author's 
co-workers there, Dr. P. E. Guenther and Mr. R. H. 
Wasserman, is gratefully acknowledged. 

1M. M. Munk, Phys. Rev. 72, 176 (A) (1947). 

2 B. L. Hicks, Phys. Rev. 71, 476 (A) (1947); B.R.L. Report No. 633, 


1947. Also Hicks, Guenther, and Wasserman, Quart. App. Math. 
(October, 1947). 





The Radial Density Variation of Gases and 
Vapors in a Centrifugal Field* 
J. W. Beams 
University of Virginia, Charlottesville, Virginia 
July 21, 1947 
F the equation of state of a gas or vapor is expressed by 
the well-known relation 
d @ 
re ae 
M + M 
where p is the pressure, T the absolute temperature, d the 
density, M the molecular weight, R the gas constant and 
where B, C, etc., are functions of T and are characteristic 
of the gas or vapor, then it can be shown theoretically that, 
under equilibrium conditions in a hollow centrifuge rotor 
with inside radius r spinning N r.p.s., 


p=RT( +05 +---) 


2eMN%?  d, 2B 3C 
————- slog —- +—(d,— ——(df—dg)+--+, (1 
RE 8G, + 3p Orde) +s (dt — det) + (1) 


where d, and do are the densities at the radius r and the 
axis, respectively. Since this latter equation may be used for 
determining the molecular weight of a gas or vapor, a 
precise experimental test of the relation has been under- 
taken. The apparatus consists of an air driven vacuum-type 
centrifuge which is similar to those previously described." 
The rotor is a short accurately machined hollow cylinder 
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(2 cm high and 4.4 cm radius) with flat ends and is sup- 
ported in the vacuum chamber by a thin flexible hollow 
shaft. The hollow shaft allows the pressure at the axis of 
the hollow spinning rotor to be measured at various rotor 
speeds. Measurements of this axis pressure were made at 
intervals of 100 r.p.s. up to 1200 r.p.s. with nitrogen, oxy- 
gen, and CO; gases, respectively. With these data, together 
with an accurate knowledge of the dimensions of the hollow 
rotor and auxiliary pressure measuring apparatus, the 
temperature, and the molecular weights of the respective 
gases, Eq. (1), could be tested. The results obtained show 
that this equation holds within the limit of experimental 


error. When Eq. (1) was used to determine the molecular 


weights of the nitrogen, oxygen, and CO: gases respectively, 
the values found agreed with the accepted values within 
considerably less than one percent. 

Experiments also were carried out on pure ether vapor in 
equilibrium with the liquid phase in the periphery of the 
rotor. Ether was distilled into the evacuated rotor while it 
was spinning until about a cubic centimeter of the liquid 
condensed at the periphery. The pressure of the vapor at 
the axis was then measured as a function of the rotor speed 
and Eq. (1) was found to hold within the experimental error 
of roughly one percent. Consequently the molecular weight 
of a gas or vapor can be measured by the centrifuge method 
with a precision of about one percent. In the case of a 
radioactive gas or vapor, the molecular weight may be 
determined by measuring the amount of radioactivity at 
the axis of the rotor as a function of the rotor speed. The 
advantage of this centrifuge method for determining the 
molecular weights of radioactive gases or vapors is that only 
an extremely minute amount of the radioactive material is 
required, especially when the material can be mixed with a 
gas such as COs, nitrogen or helium. 


* This work was supported by the Bureau of Ordnance, U. S. Navy 
under Contract NOrd-7873. 
1J. W. Beams, Rev. Mod. Phys. 10, 249 (1938). 





The Motion of a Particle in an Electromagnetic 
and Gravitational Field 


H. C. CorBEN AND MULAIKA CORBEN 
Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
July 21, 1947 


ITTING gravitational theory into the framework of 

special relativity theory is a problem that has proved 
practically insoluble. At the same time, the incorporation 
of electromagnetic phenomena into general relativity 
theory has proved extremely complex. We want to point 
out that the equations of motion for a particle in an 
electromagnetic and gravitational field may be derived 
from a simple Lagrangian which is based on an extended 
principle of special relativity theory. 

We consider the position and velocity space of the 
particle g*=(x*, ilv“), (a=1---8, w=1---4) where / isa 
constant of the dimensions of length, and v“=2*, a dot 
denoting differentiation with respect to ds = (dx“dx,)!. For 
any four-vector B‘= B,, Bi = —B;(j=1, 2, 3). We now de- 
fine an 8-invariant Lagrangian 


L=Agdg*/dr (a=1---8) 
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where dr = (dgqdq*)t and Ag is an 8-vector describing the 
interaction between the particle and the field. We write 
Aq=(exy, imle,) where x,y, €4 are two four-vectors and 
€, m are constants characteristic of the particle. Then 


& = Ldr/ds => ex,v" — me," 
Hence we define! 


sae) 
Pu dv" = ds\av*/" 


Upyss = —P, = dL /av". 


The variational principle 
6 { Lar=3 { £ds=0 


then yields p, =0,£, where 0,=0/dx". 

In general the Ag may depend upon both the position 
and velocity coordinates, but we consider here only the 
special case in which they are functions of the x# only, 
Then 


Pu=meytex,; Py=me,; 
Pu =ev’ Oyky au mv’) yey. 


The Hamiltonian may be defined by 


H = pa(dq*/dr) —L =(pa—Aa)(dq*/dr) 
so that 
dg*/dr=0H/dpa; dpa/dr = —dH/dq". 
Using d/ds =v’d, for the special case considered here, we 
may rewrite the equations of motion thus 


d/ds[v" (Auer +94) ] = (e/m)v" (Ouk» <—_ OpKy)+U'V" OOo. 
We now define 
Lur>= Oper t+ Oey ; Suv = Oyky —OyKy, 


so that gu», fy» are 4-tensors, respectively symmetrical and 
antisymmetrical, derivable from the 8-vector Aq. Thus 


Sur(dv”/ds) = (e/m)f yyv” + 4v'v" (OuSve — Pek uv — PvE ya), 


which suggests that we may interpret the f,, as the 
electromagnetic field strengths and the gy, as the gravi- 
tational potentials. For the general case in which the A, 
involves v* explicitly, other types of fields appear in these 
equations as well. 

Alternatively we have 


L£+ m(d/ds)(eyv") = exyv" + 4mgyyv"v" 


from which, clearly, the equations of motion also follow. 
We also have KH =H(dr/ds)=4mg,,v"v". We note that for 
€u = 4Xy, Zit = —1, 2c =1, Sux = O(u ~ v); then if Pu=PutP,, 
it follows that P,=mv,+exy, dP,/dt=ev’d,x,, the usual 
electromagnetic equations. 

The above results suggest that the equations of motion 
of a particle are covariant for arbitrary rotations of its 
position and velocity space. Such rotations are found to 
relate observers who are spinning or accelerated with re- 
spect to each other, and clearly mix the electromagnetic 
and gravitational components of Aq. Analysis of this group 
of rotations and of the field equations for Aq will be 
published soon. 


1E. T. Whittaker, Analytical Dynamics (Dover Press, New York, 
1944), p. 266. 
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X-Rays Associated with U**** 


R. L. MACKLIN AND G. B. KNIGHT 


tories, Carbide and Carbon Chemicals Corporation, 
Research a Oak Ridge, Tennessee 


July 16, 1947 


E should like to report the observation of x-rays 

associated with U***. L and M x-rays were found to 
occur spontaneously, the yield of the former, per alpha- 
particle, being between one-third and two. Detection was 
accomplished by means of a thin-walled Geiger counter. 
The efficiency of the tube was between 0.25 and 1 percent, 
making the quantitative determination of yield uncertain 
as indicated. Absorption curves in aluminum were meas- 
ured to determine the energy of the L radiation. Film and 
counter wall absorption made it impossible to obtain much 
information on the M x-rays. 

The source of the x-rays was a thin film of UsO0, contain- 
ing 0.245 mg of U** plated on a 22 cm? disk. Other ma- 
terials present contributed only 4 percent of the total 
alpha-activity. The U;Os film was prepared and counted 
within four hours after the removal of the natural decay 
products (UX;, UXe, UY, UZ). Figure 1, curve A, shows 
the results of an aluminum absorption experiment run on 
this source. In addition to the points shown, one millimeter 
of lead as absorber gave less than three counts per minute 
above background. The chief component found shows an 
absorption characteristic of 15.0+0.5 kev x-rays, resembling 
the L x-rays of thorium or uranium as closely as we were 





i 
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Fic. 1. Absorption*%of x-rays from U™ by aluminum. 
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able to determine. The M x-ray has not been identified so 
decisively but the presence of a soft component with the 
appropriate absorption coefficient has been shown in this 
work (Fig. 1, curve A) and in earlier experiments done in 
cooperation with Mr. L. E. Glendenin of the Monsanto 
Chemical Company. The absence of detectable K radiation 
was particularly notable in this work. 

A second set of aluminum absorption measurements was 
made using a solution of the same quantity of uranium as 
nitrate in 20 ml of solution contained in a shallow dish and 
of uniform depth, 0.78 cm. These measurements (Fig. 1, 
curve B) show no change in the absorption characteristics 
of the radiation. The intensity relative to curve A is 0.54. 
Calculating the effective intensity for 15-kev radiation 
from point sources distributed uniformly through the solu- 
tion gives 0.51. This experiment together with the high 
yield of x-rays leads us to conclude that we are not dealing 
with a delta-ray effect. 

Our interpretation is that the L and M x-rays observed 
arise from rearrangement of the electrons in the recoil 
atoms, after alpha-emission. The energy of these recoils is 
readily calculated to be about 82 kev by applying the law 
of conservation of momentum and the measured alpha- 
energy' of 4.76 Mev. This energy is insufficient to excite the 
K shell electrons (K limit = 100.4 kev),? but more than 
adequate to ionize the L and higher orbital electrons (Le: 
limit = 16.3 kev).? An adequate quantum-mechanical treat- 
ment of such excitation by recoil has not come to our 
attention. 

Dr. S. DeBenedetti has pointed out the failure of I. Curie 
and F. Joliot* to observe K shell x-rays from polonium, 
though they did observe L and M x-rays. Other natural 
alpha-emitters are accompanied by interfering beta- and 
gamma-radiation to an extent sufficient to make observa- 
tion of these recoil x-rays improbable. 

*“This document is based on work performed under Contract 
No. W-7405-eng-26 for the Manhattan Project and the information 
contained therein will appear in Division IV of the Manhattan Project 


Technical Series as part of the contribution of the Carbide and Carbon 
Chemicals Corporation.” 

1G. T. Seaborg, Rev. Mod. 16, 1 (1944). 

2 Int. Crit. Tab., Vol. VI, p. 

* 1. Curie and F. Joliot, J. - aoe et rad. 7, II, 20 (1931). 





The Nucleus as a Crystalline Solid 


J. G. Wrxans 
University of Wisconsin, Madison, Wisconsin 
June 26, 1947 


N a previous letter,' it was pointed out that the relation 
between mass and charge for stable nuclei was that to 

be expected on the basis of a crystalline-solid model for 
the nucleus with the protons on the surface. On the basis 
of this model, it was proposed that nuclear fission might 
result from the splitting of a nearly perfect crystal along a 
cleavage plane by the entrance of a neutron. Davison and 
Watson? have raised an objection to this model and mecha- 
nism of fission on the grounds that the unequal mass 
fragments after fission should be expected to contain a 
different charge to mass ratio from each other, and from 
the charge to mass ratio for other elements. They point 
out that experimental observations of fission fragments 





show that these fragments possess the same charge to 
mass ratio. They therefore conclude that the mechanism 
of splitting a crystal will not account for fission, and that 
the charge cannot be considered as residing on the surface 
of the nucleus. 

The objection raised by Davison and Watson on the 
grounds that fission fragments have equal charge to mass 
ratio seems not to be valid, since the fragments which are 
observed after fission are observed a long enough time 
after the splitting to permit a reorganization of fragments 
to obtain the charge to mass ratio corresponding to the 
resulting radioactive isotopes. The rearrangement of frag- 
ments may be considered to begin as the fragments start 
to separate. They can thus contribute their change in mass 
defect to the kinetic energy of the fragments. 

The crystalline-solid model offers a ready explanation 
for the fact that splitting occurs into fragments of unequal 
mass. If a nucleus of uranium (235) or plutonium (239) is 
considered as made up of nuclear particles in spherical close 
packing, then a crystal so formed can be considered as 
consisting of a central plane of symmetry with nuclear 
particles above and below in equal numbers, and with 
protons above and below the plane as well as in the outside 
rim of the plane. A model for uranium or plutonium would 
provide approximately 20 or 22 protons around the edge 
of the central plane of symmetry. 

If such a crystalline solid is split along a cleavage plane, 
it must be split above or below, but not through the middle 
of the plane of symmetry. This means that the fragments 
must have unequal rather than equal masses. The number 
of protons which would be required to be associated with 
the fragments is, for uranium (235), 36 and 56 for a clean 
split just above or below the plane of symmetry. This 
gives fission fragments krypton and barium. The number 
of neutrons in these two fragments would be such that the 
lighter fragment, krypton, would be left in a highly radio- 
active state and the heavier fragment, barium, would be 
left in a highly radioactive state after releasing about four 
neutrons. The various fission fragments observed could be 
associated with the splitting of the crystal with various 
parts of different planes going with the separate fragments. 

Another aspect or property of nuclei which also indicates 
that the nucleus may be considered as a crystalline solid 
is illustrated in Table I. In this table are listed the atomic 
numbers of the elements which have a given number of 
stable isotopes. Many elements consist of only a single 
isotope. All of these elements have an odd atomic number 
and an odd mass number. The elements of odd atomic 
number never have more than two stable isotopes, while 
the elements of even atomic number may have several 
isotopes. In Table I it is shown that atoms may be arranged 
in periods, where in some periods the multiplicity of stable 
isotopes alternates between one and several. In other 
periods it alternates between two and several. The known 
stable isotopes then form a periodic system shown in 
Table I. In this periodic system there are five periods in 
which, with successively increasing atomic number, the 
number of stable isotopes alternates between one and 
several, and these,periods are six atomic numbers long. 
There are four periods, also six atomic numbers long, in 
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which the number of stable isotopes alternates between 
two and several. There are three periods, each eight atomic 
numbers long, in which the number of stable isotopes 
alternates between two and several. All the elements {it 
into this periodic system except atomic numbers 4, 71, 73, 
and 83. According to this table it is to be expected that 
these elements, which are listed as having only one stable 
isotope,* may be expected to have two stable isotopes. 

The periodicity of six in atomic nuclei, as shown jn 
Table I, as well as the existence of elements of odd atomic 
number with only one or at the most two stable isotopes is 
readily accounted for on the basis of a crystalline-solid 
model for the nucleus. We may consider that the plane of 
symmetry of the nucleus contains an odd number of 
protons. Those elements which consist of only a single 
stable isotope may be considered as the ones with perfect 
symmetry of both mass and charge about the plane of 
symmetry containing an odd number of particles. Ele- 
ments which contain several stable isotopes may be con- 
sidered as those which are not perfectly symmetrical in 
both mass and charge, and can contain, for a given number 
of protons, a variable number of neutrons. 


1J. G. Winans, Phys. Rev. 71, 379 (1947). 

2 B. Davison and W. H. Watson, Phys. Rev. 71, 742 (1947). 

3 J. M. Cork, Radioactivity and Nuclear Physics (Edwards Brothers 
Inc., Ann Arbor, Michigan, 1946), Table of Isotopes. 





X-Ray Wave-Length Standards 


ELIZABETH ARMSTRONG Woop 


American Society for X-Ray and Electron Diffraction, Bell Telephone 
Laboratories, Murray Hill, New Jersey 


July 21, 1947 


-RAY wave-lengths have been expressed in x units. 
The x unit is defined in terms of the calcite spacing 

and is nearly 10-" cm, but is now known to differ from 
10-" cm by about 0.2 percent. During the last twenty-five 
years x-ray diffraction workers have expressed x-ray wave- 
lengths and crystal dimensions in terms of a unit which 
was 1000x units, but instead of calling it 1000x units have 
erroneously called it an Angstrém unit. In recent years, 
the x-ray diffraction groups have agreed to use the term 
kilo x unit (abbreviated kx) in place of the incorrectly 
used Angstrim unit, until agreement was reached on the 
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best conversion factor to use for converting from kx to 
Angstrom units. Agreement on the factor has now been 
_ of the American Society for x-ray and 
Electron Diffraction, I have been instructed to call the 
attention of American x-ray workers to the following 
announcement which appeared in the January, 1947, 
issue of the Journal of Scientific Instruments. Because of 
its importance it is here reproduced in its entirety. 


The Conversion Factor for kx Units to Angstrém Units 


At the annual conference of the X-ray Analysis Group 
of the Institute of Physics in July 1946 it was announced 
that agreement had been reached concerning the factor 
for converting measurements in kx units to Angstrém 
units. The factor agreed upon, after consultation with the 
American Society for X-ray and Electron Diffraction and 
Professor Siegbahn was 1.00202. This factor is probably 
correct to 0.003 percent. Since wave-lengths in X-units 
have been measured to an accuracy of 0.001 percent, the 
wave-lengths in Angstrém units can be taken as accurate 
to 0.004 percent in general. 

The following is a list of values of wave-lengths in 
Angstrém units of certain emission lines and absorption 
edges in common use. The column headed Ka gives the 
mean value of Ka; and Kaz, Ka; being allowed twice the 


weight of Kay. 
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Current values of the physical constants, such as those 
quoted by Birge in the 1941 volume of the Physical 
Society's Reports on Progress in Physics, should be used in 
conjunction with these wave-lengths. In particular density 
p is given by the equation 


p= 1.660202 4/V, 


where YA is the sum of the atomic weights of the atoms in 
the unit cell, and V is the volume of the unit cell in A*. 

















Absorption 

Ka Ka: Ka KA. edge 
Cr 2.28962 2.29352 2.2909 2.08479 2.0701 
Mn 2.10174 2.10570 2.1031 1.91016 1.8954 
Fe 1.93597 1.93991 1.9373 1.75654 1.7429 
Co 1.78890 1.79279 1.7902 ~ 1.62073 1.6072 
Ni 1.65783 1.66168 1.6591 1.50008 1.4869 
Cu 1.54050 1.54434 1.5418 1.39217 1.3802 
Zn 1.43510 1.43894 1.4364 1.29520 1.2831 
Mo 0.70926 0.71354 0.7107 0.63225 0.6197 
Rh 0.61326 0.61762 0.6147 0.54559 0.5341 
Pd 0.58545 0.58982 0.5869 0.52052 0.5090 
Ag 0.55941 0.56381 


0.5609 0.49701 0.4855 








It is recommended that in any published work the values 
of the wave-lengths used should be explicitly stated. 
W. L. Bracec 
Chairman, X-Ray Analysis Group of 
the Institute of Physics (England ) 





